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Journal of 
DRAU! ICS DIV ISION 


_ FLOOD- FREQUENCY IN THE PACIFIC NORTHWEST 


This p paper de describes, in general wang a method of of ¢ ‘determining th the mag- — 
nitude and frequency | of expected floods: that is applicable to any area in the 
‘Pacific Northwest of the United States. The development of a composite fre- 
quency curve applicable to— any site in a large homogeneous region | is ex-_ 
- amined, _ and the e effect of mean elevation and s ize of drainage area on 
. of the frequency curves in designated parts of the Pacific Northwest _ 


is demonstrated. A formula for determining the mean annual flood for any 
site on the west side of the Cascade Range in Washington is shown, and the 
drainage-basin characteristics used in its derivation are listed. flood- 
4 frequency curve can be drawn for any given site in the Pacific Northwest — 
Be Ps the results obtained from the mean annual flood formula and _ 
the composite frequency curve that are both applicable to that site. Com-— 
_ ments about limitations on the use of results obtained by this method are 


* 


_ “The purpose of this paper is to present a method of correlating stream-_ 
flow: data and basin characteristics to obtain information relative to a 


and frequency of flood discharges in the Pacific Northwest. The 7 


method reflects the latest developments ‘based on a continuing study of the 
subject | by engineers © of the United States Geological Survey, Dept. of the 


"ie Note.—Discussion open until April 1, 1961, - To extend the closing date one month, “a 
_ written request must be filed with the Executive Secretary, ASCE, This paper is part 
_ of the copyrighted Journal of the Hydraulics Division, Proceedings of 1 a aeatesil So- 
eiety of Civil Engineers, Vol. 86, No. HY 9, November, 1960, jal 
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November, 


Interior (UsGs). It is expected that the method w will be rev veviaed from time meto 
time as new ideas are conceived andtested. 
The discharge records collected by the ‘USGS and others are the data 
“on which flood- frequency studies are based. The records from about 800 
gaging stations in the Pacific Northwest were of sufficient length to furnish 
. data for deriving the flood- frequency curves. A curve based ona gaging- 
station record, however, applies only to that one particular site, whereas, 
generally, the information is desired at an ungaged point, frequently on an 
ungaged stream. The two parts to the problem are (1) computation of flood — 
frequency and (2) of the over abasin 
In general, of most engineering, studies are met by information 


‘their. drainage structures to. withstand floods that 
~ expected to occur, on the average, once in 50 yrs whereas their secondary 
highway structures are designed to withstand floods that: that are expected to. 
‘The Pacific Northwest is defined as the ‘Columbia River basin and the 
western slope of the Cascade Range in Washington and Oregon, ‘This area sal = 


designated as Parts 12, 13, and 14 by the 


‘FLOOD FREQUENCY ATA GAGING STATION 


‘The first step in a flood is to select the gaging stations to 
= There should be at least 10 yr of record at each station, although | : 


records as short as 5 yr are used when no others are available. - Adjust-— 


ments should be made for storage, or other artificial factors, that, tend to 


modify flood discharges significantly, or the record should not be used, 
Flood Series.—There are two kinds of flood series to be considered in > 4 


plotting frequency curves, the annual-flood series and the partial-duration — 


Series. The two methods give essentially the same results for recurrence 
intervals greater than 10 yr. A relationship between the two series hasbeen 
4 _ derived by Pegg een of the USGS, so that, when necessary, the annual- flood 


a recurrence ‘interval of 2. 0 yr the ‘annual- -flood series would have 
a recurrence interval of 1.45 yr in the partial-duration series. At a recurr- i 
ence interval of 10.5 yr the comparative value would be 10.0 yr and at 50.5 
50.0 yr. The -annual-flood method is used in these studies because of its 


‘Plotting Procedure. —Recurrence intervals are from the formula 


: . in which T is recurrence interval in years, N is the number of years of rec- 
ord, and M is the relative magnitude of the event beginning with the highest — 


Annual floods are plotted on a special form developed for analysis of flood 
; ae frequencies by the GumbelS method. Discharge is plotted as ordinate and time 


ican Geophysical Union Transactions, Vol, 30, 1949, pp. 879-881.” 

_ i... “Floods Estimated by the Probability Method, ” by E. J. Gumbel, Engineering News. 


“Annual Floods and the Partial-Duration Flood Series,” by W. B. ‘Langbein, Amer- 
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FLOOD- FREQUENCY 
as abcissa sa. Both logarithmic ordinate scales 2 are ‘used, depending 
on which plot more nearly approximates a ‘straight | line. An example of each oa 
kind of plot is shown in Fig. cs . The same data are plotted here to show the ~ 
aracteristic difference in the shape of the finalcurves. 
a: Fitting Frequency Graphs .— After the annual floods are plotted on a fre-_ 
quency diagram, a curve is fitted to the data. Considering the fact that coma 
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1.—MAGNITUDE AND FREQUENCY OF ANNUAL FLOODS. 


streamflow r records are less than 35 } yr inlength, and analytically fitted func- 
has little value. Therefore, treatment only is applied, and 


_ Flood- frequency compilations, prepared in accordance with 


_ described, | result in a curve of flood magnitude ° versus recurrence intervals 

for- each station. studied. The curve is not applicable directly to ungaged 
areas. In most instances, the flood records are short and time- -sampling © 
errors may be large. That is, the flows sampled by the period of record 
may not be representative of conditions for long periods. Also, the records — 


are, generally, for different periods of time. By combining the records for "4 4 


all stations ina . drainage basin or alarger region, t ime- -sampling and areal- 

sampling errors are reduced, and the resulting composite curves are ap- 

plicable tc to ungaged areas. However curves show 
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various recurrence intervals in terms” of ‘mean annual flood, so the latter 
EE be computed before the composite curve can be applied to an ompaee 


iin records, than five records each 20 yr long would provide a 100-yr re- 

cord. Unfortunately, this cannot be done because the data are not independent. 
In a homogeneous region, many of the gaging stations record the same flood © 
event so we simply have a20-yr record at each of the five stations. However : 
since stations are measuring the same events, they may be averaged to- 
gether to ‘provide a better measure of the frequency characteristics. ‘os a 


of the slopes of the individual frequency curves. . The slope usedis the ratio of 
a ie yr to 2.33-yr flood. According to Gumbel’s theory of extreme values, 7 
the mean of bhp number of annual floods should equal the flood with re- 


Base. | that floods for one station can- 
not be _ compared with floods for another station with a different length or 
eriod of record. In studies where mean annual floods are used as a correla- 


tion factor, _ the records must be comparable so they must be adjusted to the 


¥ 


_ The base time period may be chosen equal to any period desired for which — 
ecords are available, but, ingeneral, itis best to use the period of the longest _ 
ecords. A37-yr base, 1921-57, was used for Idaho, and a 46-yr base, 1912- a 
was used for Oregon, Washington, and western Montana. = 
Mean Annual Flood.—The mean used for each frequency distribution is the 
_ ‘Sraphical mean determined by the intersection of the visually best fitting fre 
quency > line with the vertical line corresponding tothe 2.33-yr recurrence a 
interval. The arithmetic mean for relatively short periods of record can be _ 
greatly infaenved by the chance inclusion of a major flood, The graphical © 
‘mean is preferred because it is more stable, with | rs ' weight given to the 
floods than to the extreme floods. 


a the purpose of computing an average frequency curve applicable to that region. 


| For each station, the flood discharges for several recurrence intervals are “eg 
e determined from the individual frequency curves. These discharges are divided bi ; 
js by the mean annual flood inorder toplace them ona dimensionless basis. The — 
: an value of the ratios is thendetermined for each of the selected recur- 
rence intervals. Each median flood ratio is plottedto its corresponding recur- 
‘ rence interval on a frequency chart, and a frequency curve is drawn, — 
curve, showing discharge in ratio to the mean annual flood versus recurrence 7 
_ interval in years, is based on all signifieant discharge records available and > 


4 may be considered as representing the most likely flood-frequency values for 


places in the region. Gee Fig.2). 


7 Adjustments ‘to Composite Curves. — Recent ‘Studies of streams in Idaho 

have shown that, in some areas, the shape of the curve changes with elevation. 
Also, some investigators have found wat the size of the <rainage area has an 
effect on the shape of the frequency | curve, with 
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Elevation Ac —In as the Clearwater River 
_ basin in Idaho, there is a definite relationship between flood discharge son 
Be “increases. For example, a 50-yr flood for a stream, with an average basin — _ 
elevation of 2,000 ft is 36% larger than for one of the same size area with an 
F _ average basin elevation of 4,000 ft, and is about twice as } large as for a stream | 
= this area at an elevation of 6,000 ft. Fig. 3isa family of curves valk 
_ the variation of flood ratio with recurrence interval and elevation,  — 
Area use of the dimensionless flood- frequency curve 
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- River region, the frequency curve is not affected by variation in size of — 
age areas up to about 200 sq miles, but for areas greater than 200 sq miles | 

5 the ratios of the larger floods tothe mean annual flood become smaller as the ‘ 
increases. ‘The effect is not as pronounced in the Spokane River—Pend 
Oreille River area as in some other regions. Even so, the 50-yr ratio of 
2.05 determined for an area of 200 m: miles is reduced to 1.62 for 20,000 sq 


subtracted the flood of the c composite 


to obtain flood ratios are applicable to size area being 


irequency curve, However, an eliect has been foun | 
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FLOOD- -FREQUENCY 


— Derivation of Formula for Mean Annual Flood.— A plot of drainage area 


against mean annual flood in ‘any region usually results in a graph in which 
the standard error of the plotted points about the regression line is large, 
sometimes 100% or more. Sucha large standarderror indicates that drainage- 

_ basin characteristics other ‘than drainage area may effect the size of the mean _ 

annual flood. Sometimes, it is possible to obtain a satisfactory correlation 
_ by dividing the region into hydrologic zones in which the runoff per square _ 

i mile is similar. If the standard error of the plotted points for a hydrologic — 7 
zone still seems excessive, other drainage basin characteristics n must be a 
In this study, the mantare error of drainage area plotted against mean — 
annual flood was excessive, about 90%, so other drainage-basin factors 4 
were considered. Mean elevation, mean annual runoff, area of lages and 
- ponds, precipitation, and channel slope were some of the additional factors 
that were found to improve the correlations. Graphical-linear-multiple cor- 

- relations were made to determine those parameters that were most signif- 
icant before entering into mathematical correlations. 
In the area designated as “the west side of the Cascade Range in Wash- fs 
ington” drainage area, mean annual runoff and area of lakes and ponds were ~ 

found to affect the mean annual flood most significantly and were used 
independent variables in a correlation with mean annual flood. Residual de- 
-viation of about mm. magnitude remaining in these correlations were found — 


“defined basin characteristic such as geology. A map was prepared on which ff 
; residual deviations were shown. The map was divided into areas by ——— ; 


Ge lines enveloping residuals of about equal size. The median residual for each am q 
area was used to represent a geographical factor for thatarea. = Asaie : 
‘The error “of the ‘correlation, before introducing geo- 


about 0. After assigning for geographical location and reducing — 
the degrees « of freedom by one for each zone used, the standard error was — 
reduced + about 24%, and the coefficient of correlation increased to 0.99. a 
A standard ecver 06 24% indicates that about two-thirds of the time a mean ” 
annual flood computed by the method described s should be within 24% of the the \ 

_ Formulas for other areas in the Pacific Northwest, such as western Ore- 
gon and the Snake River basin, were developed in a similar manner Sarr tee a 


Formula and —The formula derived for western Washington 


F 


in which Q refers to the mean annual flood in cubic feet per second, A denotes fr, 


the drainage area in square miles, R is the mean annual runoff in inches for 
= period 1930-57, L refers to the area of lakes and ponds in percent of Grainag 2 
area, and Gisthe geographical factor, 
_ Eq. 2 may be solved by use of logarithms or by slide rule, but it is also 
presented in nomograph form for ease of application. Two nomographs were 3 
_ designed for the above formula, one for drainage areas of 0.1 to 10 sq miles 
—" another for areas of 10 to 1,000 sq miles. The nomograph for small , 
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Derivation of Frequency Curve for. Any Site. _—The — »for deriving 
a flood- frequency curve for any site ina region is as follows: (1) Obtain the 
; _ drainage area, mean annual runoff, area of lakes and ponds, and other ono 
acteristics, as necessary, from appropriate maps (Mean annual runoff is 
- tained from a map showing isopleths based on gaging-station records.); (2) use 
the nomograph or formula to determine the mean annual flood; (3) select the | 
proper flood-ratio curve and determine ratios to mean annual flood for se- 
lected recurrence intervals up | to 50 y yr; and (4) multiply { these ratios by the © 
mean annual flood. computed in step 2 and rm against recurrence interv endl 
_to define the frequency curve. | 
aw The methods described may be used to predict the magnitude of floods nl ; 
given return intervals on streams in the Pacific Northwest. However, they | 
must be used with full knowledge of their limitations. = = ~~ eT. 
1. The term “return interval”, when considering annual floods, represents 
the average length of ‘time, in years, between the individual occurrences of - 
flood peaks of a given magnitude over a very long period. It does not indicate 
that the flood will occur on a specific time schedule, 50-yr flood is one that — 
a 2% chance of occurring in any year. fie aga - 
i 2. . The relationships shown in this report were developed from data 
‘ - gathered on the natural flows of representative streams, and should not be 
: applied to streams on which the flood-flow potential has been altered by the 
works of man, or to streams with | unique drainage-basin characteristics. 
/ Many~ of the smaller streams in the semiarid areas of eastern Washington | 
- Oregon are essentially dry except when they experience flash floods — 
= caused by local thunderstorms, This type of stream does not lend itself to the | 
methods of of analysis” presented herein | but will require special studies. 
tote Extrapolation beyond the areal limits of the basic data is not advised, 
- _ The mean annual flood formulas should not be applied to drainage areas ope 
- outside the regions for which they were derived. = 
bi ron * The length of existing gaging-station records does not justify extending 
a frequency curves beyond the 50-yr recurrence interval. Extrapolation of _ 
the curves beyond this limit may introduce considerable error. However, 
such extrapolation may provide the best available estimate Paw aanten. 
i The confidence that may be placed in the accuracy of the results obtained *) 
by using these data is related to the areal density of the basic data used. ‘a 
The denser the gaging sites, the more confidently the frequency curves “oni 
be defined, and the more accurately. the geographical and runoff factors for _ 
the flood formulas can be determined. 


structures ‘commonly built | in or pith stream channels. Where design re- 
quires consideration of recurrence intervals” greater than 50 yr, or the de- 


SUMMARY AND CONC LUSIONS. 


: 
A method has bee ignitude and 
+s; , of floods to be expected at any site in the Pacific Northwest, with exceptions [yi 


8s noted in ‘item 2 of the limitations. Flood-frequency curves for i 
stations are combined ‘several regional flood-frequency curves The 
vegions include stations for which Slopes of individual flood- frequency 
7 curves are similar. The regional curves show flood magnitudes in terms of 
| _ ‘Mean annual floods. To obtain the frequency curve for an ungaged site, it is 
: a only necessary to compute the mean annual flood for that site and convert 
i the appropriate regional frequency curve to flow in cubic feet per second by 
_ multiplying the ordinates of that curve by the mean annual flood. = | 
-.. Magnitudes of mean annual floods have significant correlations with drain- 3s) 
age area, mean annual runoff, area of lakes and ponds and in some places with Es a 
: z elevation. After all these factors are taken into account in multiple correla- 
7 [- tions, large residual standard errors still exist. These residuals tend to be of | 
about equal magnitude in given geographical areas, suggesting that a factor 
such as geology of the area has a significant effect. So-called geographic 
ors based on residual errors are 
gg ‘The formula for m¢ mean annual floods in western Washington is expressed in 
terms of drainage area, area of lakes and ponds, mean . annual runoff, and the 
"geographical factor. Similar equations have been developedfor other areas in 
i Pacific Northwest. Drainage area and area of lakes and ponds are easily _ 
Beene from topographic maps. - Mean- annual runoff is obtained from an 
4 _ isopleth map based on gaging-station records. Geographical factors are ob- 
g tained from a map that has been divided into areas having residual errors of 
ae similar magnitude. The geographical factor is based on the median residual _ 
3 - Because all of these factors are readily available, the 
The methods presented herein are subject to limitations. They cannot _ 
¥ “applied to streams in which the regimen has been so altered by the works of © i 
4 man that satisfactory adjustments for such cannot be applied. Nor can such | 
adjustments be applied to certain small streams in semi-arid eastern Wash- = 


_ ington and Oregon in which practically the sole sources of runoff are local . 
a "thunderstorms. Extensions | of frequency curves beyond the 50-yr recurrence a * 
interval are not justified, beyond that limit uld be 
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By Ber Bernard Le Méhauté! 


“i 


i 
_ general theory is ; presented to study the transmission and reflection of = 
_ gravity waves on a discontinuity to the first order of approximation ‘ahd wae 
i. wad Particular attention is given to the cases of an obstruction, a change of 


_ width, anda change of basa the gravity wave arriving a 


at an angle is studied. 


In harbors, the study of wave motion, caused by the action of incident waves _ 


from the open sea, is generally a difficult theoretical problem. The complexity 


ee. 


of wave patterns is increased by successive reflections within the harbor. 
that thes authors2,3,4 gave linear solutions for some simple cases, assuming 
that there were no energy losses in the harbor, that is no beach or wave-trap. a 
Thus, the solutions for circular,2 and sector harbors4 connected to 
the open sea by a small opening are known. The general method consisted in y 
finding a potential function @, satisfying the usual free-surface, bottom and = 
- limit conditions. wads the harbor, and combining this with the potential function | 


— 


Fn Note.—Discussion open until April 1, 1961. To extend the closing date one month, as. 
_ written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- ‘a 
ciety of Civil Engineers, Vol. 86, No. HY 9, November, 1960, 0 i 
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“Sur l’entretion des oscillations des eaux portuaires, ” by Mc Known, Publications | 74 
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“Agitation dans les ports a jetées by Bernard Le e Méhauté, 
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4 liga in the limits. . (The reader can see, for example, the very simple 
case, studied by ‘Ursell, 5 the wave transmission in infinite depth a 


_. _ Therefore, this kind of solution would not be expected to be of practical ao 


q importance in harbor engineering. Furthermore, the solutions of many prac- 
tical cases are still quantitatively unknown. In particular, the effects of the 
; harbor neighborhood shapes, and of the beaches inside the basins, are usually” 

< oe As far as the motion can be considered approximately two dimensional, _ 
there now exist linear computation methods: more simple than potential 

F. Biesel and the author used the number computation method to 
study®, 6,7 wave transmission and reflection on symmetrical obstacles and 
applied these results to study. agitation between a combination of these ‘ob- 

_ stacles, that is, in simple shape basins. In the Sogreah laboratory in France, | ; 

more complicated shape basins that had not been analyzed were 


The results of these experiments are available. 


a “nea- symmetrical discontinuities using a Similar method of computation. These 7 
Cee results can be applied to establish a general formula of the agita- ‘= 


~ ously published experimental results. Furthermore, this theory allows one 
_ to generalize the experimental results and to analyze, quantitatively, | the | effect i 
a _ Perhaps the calculus Lidl may seem too simple for a rigourous logi- — 
cal mind, but the relative > simplicity of established formulas allows engineers 
‘to use them | in practice. Above all, the assumptions are justified bye close = 
agreement between theoretical and experimental values. 
_ Notation.—The letter ‘symbols adopted for use in this paper are defined — 
5 "where: ‘they first appear, in the illustrations or in the text, and are arranged 
for convenience of reference, in Appendix II. 


second or power of the wave steepness are that is, the 


of the velocity (that is valid if the product of the steepness of the wave and the 

cube of the inverse of the relative depth is small— 


= exists athird dimensional component, it must have ‘secondary € effects, , and only 


5 “The Effect of a Vertical Barrier on V Javes | in Deep Water, : by | Ursell, Admiralty = 


— 
& 
4 
: 
4 
al 
fey. 
' 
— 
Méhaute, La Houille Blanche, Subjected to Incident Waves,” by Le 
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the main components are analyzed. Finally, s some on 
‘oo When a wave arrives on a two or three eal discontinuity, as shown © 
in Fig.1, the incident wave energy is partially transmitted, partially reflected — 
. (and sometimes partially destroyed by friction). In the close neighborhood of ‘ 
discontinuities, the waves become deformed, but the deformations disappear 
exponentially with the distance from the obstacle. 
| two-dimensional cases, motion is defined by a potential function that has 3 
the form of a linearized sum of constant coefficients (generalized constants of 
_ Fourier) of which some have an exponential form and tend rapidly to zero. 
os * three- dimensional cases, the motion becomes two-dimensional at a short 
"distance from a discontinuity if the width of the considered domain is smaller 
than a half wave length. If the width of the domain is greater than this value, 
A Jeffreys’ \ waves are theoretically superimposed on the two- ~dimensional mo- 
tion. Only in ‘cases in which the width is large with 1 respect to the wave length, — 
is it necessary to consider their effect. This, then, becomes a pure three- 
On the other hand, if three- dimensional discontinuities are not symmetrical 
about an axis in the direction of the incident waves, resonant two-dimensional 
motion can be excited in a direction perpendicular to the incident waves. __ 


« In this paper, it is assumed that all these three-dimensional effects do not ~ 


“occur, or are negligible, and the deformation of the waves disappears ex- 
-ponentially with the distance from the discontinuity (that is within one or two. 
times the depth). The is consideration of the two- dimensional 
periodical wave motion. _ 


| Consider a discontinuity in a wave flume a as one 1e of those. in Fig. 2) 
* separating the two domains (1) and (2). Incident waves from domain (1) are | 
transmitted to domain (2) and reflected in domain (1). 
in Sufficiently far from the discontinuity for the waves to have the « character- 


7 istics of a two- dimensional periodical wave, it may be assumed that the ratio 


| 


tin. we 


of the absolute value of the wave heights. 
= Similarly, it may be assumed that the ratio of the reflected wave ‘ma a 
he incident wave height has a value 8. These hypotheses follow from the as 
sumption of the linear theory, 
= On the other hand, the transmitted waves are _ of phase with the incident 
waves by an amount @ Similarly, the reflected waves have a certain phase 
} difference >» Afrom the waves that would be reflected by a perfectly vertical wall 
placed in the same position as the discontinuity. — It is assumed that this 3 pre- 


= 


vious result is applicable for a discontinuity that is 3 not strictly located. But, 
e ine this case, it is necessary to introduce an imaginary reflecting ae 
plane related to the discontinuity. It is then possible to designate the phase of 
bi the: waves as being that of their vertical amplitude in front of this reference 
} plane, whatever their direction of propagation. For a theoretical plane bar- 
— it is natural to choose it as the reference plane. When the discontinuity 


ld be ‘chosen as the reference | plane. 
the motion of the free surface at any point can be represented 


a vector in the plane of imaginaries 3), analogous to that that used 
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alternating currents in n electricity. Therefore, we we can | repre- 4 
gent the coefficients of transmission and reflection by vectors, that is by the | 
complex numbers @ and 8, representing these coefficients both by their argu- _ 
Manta: (@and ®) and by their moduli ( @ and 8 ). he be a 
When the phase and amplitude of the incident waves are represented by the 
= complex number: 7 A, the transmitted and reflected waves will be ee 
an the discontinuity is symmetrical about a a plane that is Siilieiailiiiess to the 
wave direction and two identical domains, it is clear that the 


ar 
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Harte 
the coefficients @ a Bi are of the side from wh which the 
However, in the general case, a and B are probably different for the two 
‘directions of incident waves. Therefore, we have to two transmission 
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'~Consequently,a discont ed byeight parameters: Trans- 
mission coefficient in do Reflection coefficient in domain 
; Transmissi mein domain (1)—a, ; Reflection co- 
efficient in domain (2)—, S ait — 
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the continuity and energy conservation 


tween these eight parameters will be established. The reader who wishes to 
use directly the results of this study will refer to Tables 1 through 5 in say 
_ the practical results of this study are given for a number of common cases. 


which i lisa complex of = 1 and argument: =0 
Let ho and hy be the depths of the domains — and nd (1), reapectively, and: p 


7 and 


and, to simplify the ¥ 


2m h 


sinh 2 my, hy 
bet. is pertinent to note from Eq. 4 that A = 1 in the cases in which Ww) a 
two domains, (1) and (2), have the same depth, (2) there exists ‘Shallow water, | 


e conservation of energy principle leads successively 


at + pl = te 


(a, - 8 “ay 
ke: which 0, 1, : 
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a 
to the equalities 
| 
| 
— 
— 2b, equivalent tofour algebraic 
— n of seven independent equations 
ra = = Hence, we are able to express the values of seven parameters as funetions ee 
en @€6=—S—s(«ét:séonly one, a, and the factors Z and A, that characterize the two domains 4 Se. 
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it may be shown that, 
+B 


bi D, so ) that BD = A Z a, and draw the pee of center A sa radius AD, 
re) intersecting CD at the point E. Itis easy = bat that this sketch satisfies | 
all the above conditions. ‘Thus, we have 


et 


— e considered discontinui hape, 4 
quently, it is the others from the followin (8) 
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One ca can also, without further ‘computation, ‘the following noticeable rela-_ 


Ay + By = i | 


(18) ‘a 


riatio ons” of each coefficient are easily s studied by means of 


Fig. 5. "When A is greater than 1, D and E, By and fo, and are 
reversed. (a) If A, Z tend to 0, either bya reduction of domain (2), or by an ‘ 
increase of domain (1), as shown in Fig. 5, then, — to the sketch for 


4 wee 


is total reflection withett change sign 1 (as in a closed acoustic 


pipe); 


ai 


‘There is also total reflection with change of sign, (as open 
_ (b) If AZ tends to 1, that is, if the domains (1) based ) are identical, the re-_ 
Suits published a are found, that i is 
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Even if the has an shape (Fig. 6), the effects 
on waves are > symmetrical. (This is due to the assumption that the discontinu- 
ity is local, and that the linear theory is used. 4 
_ However, this result is also true when domains (a) and (2) are different, 
with A, Z=1, (for example, as in the case shown in Fig. ce But A, Z can be 
- equal to unity | for only one value of the wave period. * vik) 
If a tends to 0 (by obstruction), A, Z being 


ss VALUES OF A AND B AS FUNCTIONS OF THE SHAPE. _- 
OF THE DISCONTINUITIES: EMPIRICAL L FORMULA 
Tt has be aie seen that all the coefficients @ and 8 are known when we are _ 
_ able to express only one (for instance ay) as a function of the shape of the dis- - = 
_ continuity. This value has been obtained with the help of theory, physical con-  __ 
_ siderations, and experimental results, some of which have been previously 
_ published. 7,8 _ (However, the formula obtained can be modified after other — 
Consider the most general possible shape of 
with at the mouth of harbors (Fig. 8), 


The coefficient ay is given by the empirical formula 2 (N) 


b 
) 
in nwhich b is the opening g between domains (1) and (2). However, this empirical 
could be justified by a number of theoretical considera-— 
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e e function on dp is is perfectly some respect to ay. This is pi physically 
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consistent. 
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7 FIG. 7. —A CASE WHERE AZ COULD BE EQUAL TO UNITY DE- all 
‘SP ITE THE CHANGE OF DEPTH AND WIDTH. =, 
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b, and hy # tho obtain 


and the in the case of a ingle change of 


ly 

‘Equating these two values leads to A and Z or =hg when = that is 


In--the case of a discontinuity with deepening and as 

- shown in Fig. 9, when A Z varies from 0 to 1, that is for instance, when the 
dimensions: of domain (2) vary from 0 to those of domain (1), a, varies from - 
1 to 2, By and By from 1 to 0, and Qo from 0 to 1, at the same time. Conse- — 

¢ quently, the wave 7 the domain (2) cannot have a height greater than twice the 


‘ height i in domain (1), and this occurs when the ratio of the dimensions of the 


‘two tends to z zero; that is when — 


It is interesting to note that the height always when the 
depth: decreases ‘Suddenly. 


qj ne Finally, it should be noted that in n shallow water L = VE h, an 
=2mh. Hence, A = A=1 * 
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dle FIG, 9.—A CASE OF MAXIMUM TRANSMISSION. 


10, —VARIATION OF THE DISTANCE BETWEEN OR- 
THOGONALS FOR A WAVE AT AN ANGLE. 
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Some authors, after Lamb, 10 have a ‘that, with the 
| 


the hand, it is difficult to introduce the an obstruction 


Indeed, if for example, one admits, as logical 
y 


"theoretical ‘relationship such as Eq. 8, and giving a. a correct for any 


(32a) 


‘the theoretical relationship of Eq. 8 In practice, 
_ the difference between the two formulas for (ay) (and (ag)) is not very _— 
portant, but it would seem better to keep the proposed formulas that satiety 
both the relationship and the considerations for any par- 
Hence, the formula of Eq.19is valid for any case, that is: 
- tion without change of characteristics of domains (1) and (2) or change of char- 
a _ Introducing the value of ay given by Eq. 19 in the theoretical formulas (Eqs. 
ws through 12) permits us to compute the value of the other coefficients; 4, il -: 


“4 a 9 “Méchanique Céleste,” by Poincaré, Torne, 131, p. 87. 
10 by Lamb, 6th ed 
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Boy ay, in the case e and, ina 1 a number of cases, 


as it is shown in Tables 1 through “re 
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When 1a wave arrives at a sudden : 
depth, a part of its energy is reflected. | The angle of the reflected w wave with 7 eS 
the line of discontinuity is the same as the angle of the incident wave, as is _ 


~ shown by Fig. 10. The distance between orthogonals does not change. == 
via Because of the change of depth, the velocity “ the transmitted wave changes 


ehemela varies. For example, if the distance between orthogonals yy is the 


(1), and is the: shallower domain (2), we — 


in which 6, and are. angles of the waves with the line of 


in the domains (1) and (2), respectively. On the — r hand, be because of the | 


change of velocity of the | 


Mik Ly 
sin 


ey Cantina: these relationships, the value of|— | is obtained as as functions c of -: 


ty) 


and conservation of energy are to 


orthogonals takes place, instead of the ratio of width of 


the domains (1) and (2), Hence, further computation, 
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ofthe 
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this case is similar to the case shown by Fig. 
oid Hence, if the domain tal we assumed to be the shallower domain, the ime 7 
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“Introducing t the ‘values leads | 


ral 


Le is always smaller than 

and a, always” has a ‘definite value. ‘tends _to be infinite if the wave 

arrives at an increase of — at an angle 59 such that (sin 2) tends to 


= 


In the same time, | tends to (Gin 52) is greater than , that is sin 


, there is total reflection at the discontinuity. In practice, such a wave — 


diffracts inthe e deeper domain. anology of this with light arriv-— = 
ing ata medium of a different refraction index is evident.! 
_ This particular solution has been considered to protect the mouth of a har- 
ber from agitation, but it does not offer a very practical s solution because of . 
the cost of maintenance by dredging. 


; _ The most complete verification of this ‘theory has been made by edie 24 
; the agitation in a basin, that is, between two discontinuities. Indeed, atheo- _ 

“Etude de la houle en profondeur finie-Note sur la réflexion totale,” by Carlotti 


Houille Blanche, and Avril, 1947, PP. 112-116, aa 
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corer to) express: completely the value of aj, it will be noticed thata re- = 
duction of depth causes an increase of distance between orthogonals. Hence, 
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‘retical formula for the a mene will we This theo- 
_ retical formula is based on the theory for waves at a discontinuity presented © .. 
= here. If sucha theoretical formula is well verified by experiments, it is con- 
sequently true. that the theory for waves ata discontinuity, the basis for this es 
further computation, is also 


ie results of some tests made at Queen’s University in Kingston, Ontari 


O, 
are presented in Fig. : It is important to recall that the theory is valid only o 


if ‘the quadratic convective inertia is small, that is if (2) (=) 
& all the tests have been carried out with small steepness and large re ative 7 
depth. It must also be remembered that, except in the case of a simple change 
of depth, the width of the channel must be smaller than L/2, = 
Ih addition, it has been assumed that the loss of energy due to friction is 
negligible. T he main way to. achieve this, is to maintain conditions prevailing 
for negligible convective inertia, that is, a small wave steepness and a . 
relative depth. Unfortunately, since the discontinuities used for these tests 
&g _ were made of simple concrete blocks placed in a wave flume, some energy was 
“im dost because of their roughness and because of the leaks oc occurring between the | 
blocks. This loss of energy is evident in the case of a complete closure of a 
flume, because the reflection coefficient is only between 0.8 and 0.9 in- a 
stead of 1. 0. However, in any case the loss of _ computed by the form- 


The theoretical curves proposed by Lamb are with the theoretical 


a values proposed in this paper. ‘Indeed, the fact that a certain amount of energy 
- was lost by friction against the rough concrete blocks and by leaks in between 
= must be considered. This especially effected the reflection coefficient. 
Jee Another occurrence of discrepancies was caused by the fact that the trans- 
~ ies wave was partially reflected on the terminal beach of the wave flume. 
= _ This effect is not negligible if the steepness of the transmitted wave is small, 
as in the most frequent case. Hence, it may be concluded that this theory has ~ 
—" verified within the accuracy | of the tests. However, it would seem rela- 
= simple to obtain some better experimental results. ‘Te 


No have yet to the value the of 


' waves ata discontinuity. The author would be happy to know if there exist any _ a 
13° 6 ver aed results that demonstrate the validity of his theory and further | 
tests” to verify or to improve his empirical formula (Eq. 19) giving ; the value 
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ora change of depth, Eq. 2 gree of ac- 5 
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b seems that many theoretical results are valid, even when the width is 


greater than dimensionless factors such as: could be taken 


may a transverse resonance. 


an int In the case of a loss of energy at the discontinuity, the relationships ob=-. — 
However, 
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‘Physically, this case could be conceived if the is 
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because of the energy loss 
i we In this case, the graphical representation has to be modified as is shown by ee 
12. The tl tend to zero, and 
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_ APPENDIX I. MATHEMATICAL DEVELOPMENT» 

Continuity Principle —We will assume that the discharge of incident 
reflected waves in domain (1) is | equal to the discharge of the transmitted 
wave in domain (2) on the discontinuity at any time. It is natural to = 
this assumption, for if this were not the case, a quantity of water would suc- — 
_ cessively be accumulated near the discontinuity and would issue from it. ‘eee 
_ The discharge is the product of the velocity times the cross sectional peg 
(depth times width). The velocity u is al function of the vertical 1 abcissa > 

Hence, the discharge could be 


Now assuming that the domains (1) and (2) have constant widths 1; and lo, | 
and constant depths hy and ho, respectively, the continuity principle is given 
by equating the the discharge in the at the discontinuity. 


Assuming that wn ‘motions in domains (1) and (2) are —— by the potential 


functions and $0, and the horizontal x-axis origin coin- 
“of 0 or asfollows: 


Hence, the continuity r nship becomes: 


and the potential function of a reflected wave of amplitude 2a x 
“os travelling seawards and having a difference of phase Bye * the ‘discontinuity. 


2 2 


— 
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— ~~ in which 1 is a constant width and h the depth of the con ; ee 
— 
ave 
potential function 18 given by the su 4 
— 


= sin + my a a+ By) . 
‘The potential function $» is that of a wave — amplitude 2a x a4, travelling 7 4 
the same direction as the incident wave a and having a phase difference Ay with Se 
the incident wave at the 


— ht - 


cosh mg (hg + y) 
a my hy ihe ha it cos (nt + + 


= “sin 
these values in Eq. 46 46 since +e 


_coshm(h+y), 
Daye 


Silay an incident wave from the equcetie direction we find 


a Already, one can see the eight will be the 
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Hence, ...60) 4x41 
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tog 
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ia On the other hand, ‘the most important loss of energy at the an is 
due to the convective inertia. Since this force is negligible, one may neglect | 
§ loss of energy and equate the incident wave energy to the sum of the re- 
flected energy and the the writing we intro- 
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"Similarly, fc for a wave even the opposite direction 


directions, with height 2 2 ag and with any 
conservation of energy principle takes the 

+ Ag ly Lo 
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1, + Ay ly Ly + A212 L202, 
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relationship ‘must be verified for any value, - - 
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- aj, successively, we find 
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solution - BQ 
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equations with eight unknowns. it is easy to eves seen 
that Eq. 64b depends ve the so ) Theref re, it is possible to express the 
eight coefficients: ay, By, B a9, Qo, Bo, as a function of only one (for 
instance a), and there exist some other — relations. From the 5 pre- 


(65a) eat 
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ontinuity Principle and Energy Conservation ‘Principle. eo conti- 

relationship we find, ‘seccessively, , for the first incident wave (Fig. (13) 


serting Eq. 69, deduced from the conservation of energy ee ae | 


“and Ba, yields 
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The solution ay a9 is not possible. Indeed, in this ‘case, Eq. 64b would. 
yield By = + 1 (+2nm) and cos py =- cos by. _ Now, from Eqs. 71 a1 and 74, 


_we could establish a relation between ay and Z ue A, and all the parameters 
_ would be known only as functions of Zand A. This is physically impossible. 


"Indeed, these parameters not on characteristics, but also = 
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wave height to incident wave height; 
coefficient of transmission from the domain (1 1) to the domain 
of transmission trom (2) the domain 


oefficient of reflection at the discontinuity: ratio (of the eee 
coefficient of reflection in the domain (1); 
a7 coefficient of reflection in the domain (2); 
phase ott the wave by with the incident 
of relative to the coefficient a1; 


change of phase relative to ) the coefficient ao; : 

phase of the reflected wave with the incident Wi 


relative to the coefficient By; 

£ relative to the coefficient Bos 

complex number of modulus a and argument Q; 

complex | number of modulus B and 


: complex n of i= and argument: 

incident wave-height in domain 


incident wave- height in 1 domain (2)5 


phase o of the incident wave in 


potential function of the ‘motion in the domain (1); 
ences function n of the ‘motion in the domain (2); 


: horizontal a axis; 
3 vertical axis, 
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EDDY FORCES ON OSCILLATING CYLINDERS 


By Alen D. K. Lairdl, and Robert Ww. Walter? 


cause large force fluctuations. 


INTRODUCTION 


— 
ecent 4,5 have that eddies along piles i in 
ocean waves may make important contributions to wave forces on piling. The a 
A present laboratory investigation | of these eddy | forces is part of a continuing 


_ research project being conducted at the University of California in Berkeley. 4 } 


Hie Note.—Discussion open until April 1, 1961. To extend the closing date one month, a iy eg 
ea written request must be filed with the Executive Secretary, ASCE, This paper is part * : 
_ of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
of Civil Engineers, Vol. 86, No. HY 9, November, 1960. 
1 Assoc. Prof., Mech, Engrg., Univ. of California, Berkeley, Calif, 
a Research Asst., Inst, of Engrg. Research, Univ. of California, Berkeley, ‘Calif, a 
3 Research Asst., Inst. of Engrg. Research Univ. of California, Berkeley, Calif. Pre- 
Design Engr., Catepillar Tractor Co. 
4 4 “Ocean Wave Forces on Circular Cylindrical Piles,” by W. L. Wiegel, K, E, me a 
James Moon, Proceedings, ASCE, Vol, 83, 1951, 
“WaterForces on Accelerated 
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“RATIOS OF AVERAGED MAXIMUM DRAG FORCES ON TEST CY LINDER * 
WITH PARALLEL NEIGHBORING CYLINDERS, TO MAXIMUM DRAG FORCE 


CYLINDER WITHOUT NEIGHBORS.® 
(A) Two | neighbors ‘symmetrically spaced to left and right of test ‘cylinder 


Abed To each side, diameter 
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©) ‘Two. neighbors ‘symmetrically ‘spaced to and right of test cylinder path 


1 19° 


. 
D) One ne neighbor | (55 runs used), d 


_ ® Positions of ‘neighbors ar are given relative to the axis of the test cylinder, in cylin- 
der diameters. Ratios were taken at the bottom of the first swing after release from rest. 
Cylinder diameter = 2 in, Amplitude of oscillation = = 18 in, 4 Cylinder = 20 


in, ‘© Amplitude of f oscillation = 6.4 4 in 
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1960 
—_ was the object of this experiment to demonstrate that wave iene ‘ona 
- circular cylinder could be influenced significantly by parallel neighboring cyl- 

: -inders, and shed these neighbors. 
EXPERIMENTAL EQUIPMENT 


sil vertical cylinders were modeled by the snclion of the cylinders p performing — 
B.: simple harmonic motion in otherwise still water. The heavily weighted four- 
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3 FT. “TEST CYLINDER 
2"( DIA. 


‘The test cylinder ‘was mounted on a strain bar 
of the pendulum frame. As listed in Table 1, the neighboring a 
3 _ inders of the same diameter as the test cylinder were spaced, with respect to eh 
& the test cylinder, from touching up to 8 diameters to the side and 18 }diameters 5 
ahead of the test in the direction of the first forward after re- 


For in. the strain bar deflections were by a ca- 


“ferent run by ‘turning the bar through 90°. on the 2- in, 
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FIG, 2,—RESISTANCE AND LIFT FORCES ON 
CYLINDER FOLLOWINGTWO 
NEIGHBORING CYLINDERS FOR COM- 
_ PARISON WITH SINGLE CYLINDER © 
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eee were detected strain gauge bridges and recorded by the oscillo- __ 
a graph. Asimultaneous trace from a 2 gaccelerometer was recorded with each 
_ ‘Tun of the one-directional bars for both the 1/2-in. and 2-in. test cylinders. = P 
About half the runs were taken with a one-directional bar, Introduction ofa | 
two-directional strain bar allowed simultaneous recording of drag and trans- — 
verse forces, but the accelerometer trace was omitted. Calibrations of the 
a strain bars were made before and after each set of runs by applying known La 
forces at the waterline of the test cylinder by means of balanced strings, pul- . 


TABLE 2,—MAXIMUM FORCES | ON 2-IN. CYLINDER, IN GRAMS, REFERRED 


12 175 159 


-125 


-* a The two o parallel 2- -in, neighboring cylinders w were spaced symmetrically to left and he 4 
‘i of the test cylinder path. Amplitude of oscillation, 18 in. Positive and negative — 
lift forces to opposite sides. Forces are maxima with test cylinder following neighbors, ¢- 


but not necessarily on first swing after release. 


To avoid the vibrations caused by suddenly applied accelerations, motion ao 
was started from maximum deflection by relaxing the tension in a hand- held 4 
: cord attached to the center of gravity of the pendulum. 
During many runs , a 16 mmcamera recorded the motion of the eddies shed 
by the cylinders, the traces of the magnetic oscillographs, the pendulum posi- 2 
— tion, , and a timer. Thus, visual impressions of the action could be checked by 
viewing film in slow motion, or frame by frame. 


eae LYSIS OF DATA 


Since the chief aim of this research was to of 
test cylinder, the forces onthe test vis itself were ‘made 
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The drag force histories for runs without neighbors 
were averaged from the oscillograph records andthe extremes were recorded 
Ri. as shown in Fig. 2(a). The curves were not corrected for velocity variations. a 
5 The drag difference curves shown in Fig ” Y _ resulting from 1 subtracting the 
eS average force curve from the force curve for the run under onsideration, a 
emphasized the effects of the neighboring cylinders. Fig. 3 shows typical time | a 
position histories of test cylinders, test cylinder eddies, and eddies from two 
cylinders showing comparison with the corresponding time force 
The neighboring cylinders further the velocity and amplitude of 
- the pendulum > motion over that caused by the test cylinder alone, For each ~ 
7 number of cylinders, the decrement of the amplitude was measured so that the 
- maximum velocity at the bottom of the first swing from release could be com- _ a 
_ puted from the theory of damped simple harmonic motion. A multiplicative — 
correction factor was computed for each ofthe cases by dividing the square of 
_ the computed maximum velocity of the no-neighbor runs, by the square of the 
2 computed maximum velocity of the runs with different numbers of neighbors. 
_ For one neighboring cylinder this factor was 1.03; for two it was 1.08. Each of 
the values in the tables has been multiplied by th the e appropriate ci correction 
The lift forces proved to be non- replicative s0 they» were not averaged. 
Typical curves are included in Figs. 2(b), (c) and @), and maximum values © 
are listed in Table 2. From successive frames of the motion picture films, 
¥ _ measurements of the positions of the test cylinder, and of the eddies shed by ; 
the test cylinder and its neighbors, resulted in curves of cylinder and eddy - a 
from the release ase point of the test cylinder as fu: actions such 


that the neighboring influenced the 


= 


a the test cylinder was in the wake from its neighbors there was a semaedi de- is 
i= crease in the average drag force to minima of 1 of about 2! 20% of single a. 
values at ‘mid- swing. When the test cylinder was not in a wake, there was an 
. increase in mid-swing average drag to maxima of over 150% of single cyl- 
inder values. Superimposed on this general increase or decrease of average 
drag were fluctuating drag and lift forces caused by eddies passing the test 
cylinder. The largest fluctuations of mid-swing lift anddrag forces were over | 
e 90% and 150%, respectively, of the maximum drag force on the test cylinder 
> without neighbors. Since the test cylinder experiences increased average drag 
FS forces when outside the wake of a a neighboring cylinder and nd maximum lift and = 
fluctuations drag forces only when it lies within the wake, there is little 
possibility of complete reinforcement of these effects. Nevertheless, reinforce- 
: eS ment did occur, as shown by the highest and lowest recorded drag forces at = 
nee. mid-swing of 1.75 and minus 0.17 times the ei mid- -swing value for r single 


2(a) and the envelope enclosing all the eneoaanention resistance forces meas- 


A 
ured is shown in Figs. 2(a) and 4, The dra coefficients corrected to infinite a 
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eylinder values 1.3 to 1.4 for vertical cylinders 
~» Gees the water surface. The corresponding lift or side force envelope is 
_ shown in Fig. 5. Data for the 1/2-in. cylinders were the same but for scale. 
- Typical resistance and lift forces on the test cylinder in the company of a 
pair of neighbors placed as in Fig. 6 make up Figs. 2(b), (c) and (@). In these, 
significant departures of the resistance forces from those for the cylinder _ 


- alone are evident. The closer the > spacing, the greater and the earlier a are _ 


6.—CYLINDER AND EDDY MOTION 


ahead of the 2-in can be as from the r ranges of 
Ce pessoa and lift forces found without neighbors. Again, the deviations from a ; 
single cylinder behavior decrease as the distance ce between the cylinders in- 7 
creases, and become unimportant at 4 | diameters. 
a ad Typical differences between the resistance forces on the 1/2-in test cyl- 7 
4 inder without neighbors and with two 1/2-in. neighbors are shown in Fig. 3 
i _ with typical test cylinder and eddy loci corresponding to these force traces. . 
The eddies from the neighbors caused the fluctuations in the force 
‘73 Because the times of shedding of the eddies from the test cylinder. do not cor- Pee 
ie respond to the force trace deviations, these deviations must be attributed to 
the of the eddies from the n the neighboring cylinders, rather than to the 


| 
= 


OSCILLATING CYLINDERS 


_— of eddies from the test cylinder. It is evident that the closest | ap- 


& : proach of the first shed outboard eddy pair to the test atenggere coincides a 


ity of the test cylinder relative to the surrounding water and must therefore be 
the reason for the negative force contribution. Subsequent outboard eddies 
cause similar negative contributions. For the run at 8 diameter spacing, - lle 
— pairs of inboard eddies could be seen onthe film. Their loci precede the — 
_ responding outboard pairs and they must, because of their direction of rotation, — 
cause the increase in drag contributions shown before each decrease. Un- 
E fortunately, on the films the inboard eddies were very difficult to follow 7 : 

. 4 cause of the surface disturbances behind the cylinders. As a result, none were . 
clearly defined for the 4 diameter and 12 diameter ‘spacing runs of Fig. 3. One 
may, however, infer their existence because inall of the runs in which fluctua- 

tions of drag force were noticeable the first deviation was an increase, All 
‘data support the facts that the shedding of eddies from the test cylinders when ) 
alone or with neighbors, caused no significant change in the resistance force, 
and that large changes were caused when the eddies from the neighboring cyl- 
a The maximum recorded lift force of 157 gm is shown in Table 2. This force 
Ls some 80% of the average single | cylinder drag force and about triple the 
average single cylinder lift force. «Lift forces of 50% of the single cylinder _ 
drag force were common. Drag forces were often increased by 50%. Thus, the 
‘lift and the drag force acting simultaneously may augment the single cylinder — e ¥ 
_ force by at least 60% and cause oscillating forces on , Soe Saeed 
me _ Often the motion pictures showed that an eddy was shed almost siaieitinite 
; ously from each side of the test cylinder without causing a liftforce. In many _ 
cases, , the eddy that formed first was not shed first. It sometimes shrank while a 
the opposite eddy grew and was shed. These starts were often followed by ulti- : 


mate shedding as in the single cylinder von Karman street. When neighbors 
_ were close, however, the test cylinder sometimes shed only one eddy, and 
sometimes none, as - shown in Fig. 3. The frequency of shedding | shedding of 
a? eddies from the neighboring cyjinders can be judged, from Fig. 3, to agree 
with that given by Landweber®, ‘The velocities of these shed eddies were 
‘measured in this and other figures for 1/2-in. . and 2- in. cylinders and were 
found to agree with Landweber’s values for eddies shed by cylinders traveling ij 
= a : a The lift force fluctuation periods agreed more closely with the single cyl-— 
inder eddy shedding period given by the von Karman vortex street anal-— 
ysis than with the periods for eddy pair shedding from the neighbors as given 
elsewhere for medium to wide spacing of the cylinders. With these spacings, 
there was little ‘measured difference _ between the lift force records and the q 


on U. 8. The David Ww. Taylor } Model Basin, Le L. ‘Landweber, Report 


| 
3). 
ie 
; of a _ were caused mostly by the shedding of eddies from the test cylinder with bat is 
™ little effect from neighboring cv unless the neighbors were within 
5 


were exerted —— as shown in 


1. The presence of cylinders | may increase the of 
by 75% or decrease ittolessthanminus 15%. 
2. Lift forces of 75% of the single cylinder drag force may result from the 7 


as The average drag force on a cylinder may be increased by as such as = = 


45% if outside, or decreased by 80% if inside, the wake of close neighbors, 
_ Fluctuating drag forces are caused mostly by eddies from neighboring 

One neighbor has approximately the same effect as two. 


6. Neighbors more than a test cylinder | cause only 


‘minor drag and lift disturbances. 
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classifications of flow in | steep, “rough, ‘channels are proposed: 
7 tranquil, tumbling, and rapid. The tranquil and rapid regimes are dominated — 
va by subcritical and oe ag flow, respectively. In the tumbling regime | 


> 


parallel to the bed slope. For such flow the depth, water area, velocity, Ae 
discharge are the same at each section of a channel reach. Whereas these _ 
conditions n may be essentially correct for open channels in which roughness 
a. elements are small, relative to depth, they are hardly correct for channels > 
a _ having large roughness elements such as occur in mountain streams or in al- ye 
luvial channels in which large dunes are formed. Existing formulas relating 
: _ the gross variables of discharge, water area, depth, bed slope, and roughness — 
that are based onthe concept of uniform flow, could hardly be expected to apply 
in the 1 latter case, especially if if regions | of of supercritical flow form over the 


Note.—Discussion open until April 1, 1961. To extend the closing date one month, a 
wetheeh request must be filed with the Executive Secretary, ASCE, This paper is part — 4 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 7 
ciety of Civil Engineers, Vol. 86, No. HY 9, November, 1960, 

_ ® Presented to the Hydraulics Division, ASCE, July 1, 1959, ‘Fort Collins, Colo -_ 

a Prof. of Civ. Enere. and Dean, Utah State Univ., 
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tops. of the this flow in in n nature, the basic 
a can probably no longer be described by simple numbers, but become Statistical. | 
The purpose of this study was to gain some insight into the basic dynamic 4 
i ‘ siaaee of flow in open channels, especially in steep channels, in which 
roughness is an appreciable portion of depth. Since conditions in nature are 
extremely complex, laboratory» flume studies under idealized conditions of 
i geometry w were made as as a first step in trying to understand the problem more 


» _ The design of the experimental program and the apparatus’ was based, among 
other things, on consideration of many informal observations in the field, es- 
pecially of Logan River located in northern Utah. This river is a typical, steep 
(50 to 75 ft per mile slope in its lower reaches), boulder-strewn mountain 
stream with discharge normally ranging annually from a minima of around 
100 cfs to a maxima approaching 2000 cfs. The writers were impressed by the | 
fact that the roughnesses commonly appeared to induce a “tumbling” phenom- a 7 

ss na whereby localized areas of supercritical flow followed by hydraulic jumps — 

 emenih. Fig. 1 shows this condition on the Logan River, that is typical of | 


Tough, steep streams observed at many other locations. 


- It seemed desirable that the flume be capable of rather large slopes, and 

the roughness elements should be such as to induce large variations in 

: depth ; and velocity. Since these latter effects were of primary importance, ll a 
because of a desire to separate them from other considerations and make them 

as simple as possible, two-dimensional roughness elements seemed appropri- — 

a ate. A second step could be to introduce three-dimensional roughnesses, and — 

vs third, to investigate similitude by using similar roughness elements of dif- _ 7 


ferent size. As a fourth step, consideration was given to the introduction of 
roughness» elements having distribution of ‘sizes, however, on 
A number of preliminary theoretical investigations were made utilizing di- 
"nemo and mechanical considerations. These indicated that some form of a 
Froude number, together with slope, some description of roughness height, > 
Be geen and intensity should be the important variables. These analyses could ‘ 
be developed to some degree for specialized ideal ‘cases in the flume, but | 
little information regarding the general case. 
This paper summarizes the principal ‘results and conclusions reached by 
= study. More detailed information, including a complete tabulation of Geta 
- Notation.— The letter symbols adopted for use in this paper are defined | 
oo they first appear, in the illustrations or in the text, and are arranged 
_ alphabetically, for of reference, in | Appendix ‘pork 
; ee Flume. — Opportunities at hand dictated the design of the water supply, 7 
‘ flume, and measuring devices. The installation was made in the out-of- -doors — 
at the site of the Utah Water Research Laboratory (Fig. 2). Water was si- 
phoned from the reservoir above the 40-ft high State Power Dam toa 1-ft- -yy- 


Pee 3 “The Dynamics of Turbulent Flow in Steep, Rough, Open Channels,” by P. K. Mo- 
hanty, thesis presented to Utah State Univ., at Logan, in 1959, in partial fulfilment of — 
the requirements for the degree of Doctor of 
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-ft- -by-7-ft. he: sabox th a 10 in. by a gate: valve at the 
headbox. A maximum discharge of 4cfs was available, 
—— The flume was 2 ft ft wide, (16 in. deep, and 64 ft long. it was constructed of | 
y plywood on longitudinal timber girders supported by ; a hinge at the center and i 
_ by two differentiai chain hoists 9 ft from each end. It could be set on any Slope 
_ between 0 and 3.5%. At the headbox connection, a sliding arrangement consist 
fe ing ofa non-stretchable piece of nylon cloth 16 ft- -by- 14 ft witha hinged and 7 
— sliding door proved leak-proof and otherwise satisfactory. r A sliding tail gate _ 
was provided at the discharge endofthe flume. = © 
‘The discharge from the flume spilled into a stilling pool int the 
earth but lined with concrete. At the end of the pool, a V-notch weir was in-- 
stalled. Depth over the weir was measured at a ‘Stilling well with a micro-. 
j hook-gage. A floating slab of plywood served as a wave Suppressor. As a 
= commercial flow-meter was installed in the Siphon pipe. 
measurements, however, were used in the a 
= - Two- dimensional roughness elements consisted of 0.3 ft- -by-0.3 ft or 0. 0185 
by-0.135 ft ‘square timber bars (or a combination of them). Three-- 


; elements were 0. 3 ft-by-0.3 ft-by-0.3 ft timber cubes. — 


_ were machine planed and fastened to the flume bottom by toenailing. * — 
_ plate was used during installation to assure proper spacing. - Roughness ele- 
ments» having finite width were chosen because of ease of attachment and be- 


cause some closer resemblance to natural , having finite width, 


7 a Measurements. .—The slope of the flume was carefully checked at 2 ft sta- 

tions u using an engineers level. Depth measurements were made by means of 
> a portable micro-point-gage. Measurement of the depth Y, over the leading ; 
a edge of the roughness was repeated at from 5 to 10 roughness stations inor- | 
i, der to check uniformity of the cyclic pattern. The average value was used in ; 
a analysis. Depth was also measured along the longitudinal center line in 1 order ; 
s to obtain the flow profile. A large number of closely spaced piezometer holes b : 
were provided along the longitudinal center line. These were connectedtoa 

_ manometer board and recorded localized piezometric head | at the boundary 7 


that usually differed appreciably from the depth of flow measured at the same — 


3 Discharge was measured using the 90° V-notch weir as described previously, — > 
and the temperature 2of the water was notedforeachrun. 
a Procedure.— After the roughness elements had been installed and the slope 
fixed, the gate valve was opened enough to allow the desired flow. After the "I 


_ flow was established, , discharge > measurements were mi made. The rubber 


meter readings and observed values of ore were recorded after allowing Lm 
sufficient time for attainment of equilibrium. 
Using a reflex camera, at least two photographs of the flow pattern were 
_ taken for each run. A journal describing the physical characteristics of flow 
- noting peculiarities, if any, was kept for each run. - The following shows 


of channel | 
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_ A total of 163 runs were “made 
Dimensional Analysis. -A number | of possibilities for measuring and char- 
_acterizing velocity and depth for the flume studies are apparent. Several of a 
_ these were tried during the investigation. However, if a cyclically uniform — 


For the purpose at hand, the depth and velocity at the front edge of 
element proved most satisfactory (Fig. 3). Letting th the 
ing velocity serve as the dependent variable, 
= K, b, L, t, S, g, 
in which K, b, L and t are, respectively, the height, top width, aay 
_ and transverse spacings center to center of the roughness elements, S is the ' 
_ cel Slope, and g and v are the acceleration of gravity and kinematic vis- = 
cosity of water, ‘respectively. The factor A represents ¢ a roughness intensity — 
factor. Since roughness elements with sharp edges were used, fixed lines of 
"separation were anticipated, and thus the form drag, and consequently, Vy 


Should be nearly independent of Reynolds number. Choosing g and wd 


The roughness factor A may be expressed as 
in which p p is the number of transverse roughness spacings per unit width of | 
the channel, and a is a shape factor. For two-dimensional roughness, p t = 0 


and A= a a K/L. For the two- roughnesses, 


Since the roughness are standard shape, a was not a variable 


of Flow.—Early observations of natural streams focused the 

; x attention of the writers on a condition of flow dominated by scattered regions | 

4 of alternate acceleration and deceleration through critical flow over the large 

—— elements. This condition was termed tumbling flow because the 

water appeared to tumble over the roughnesses. . At lower discharges, in 7 
reaches of lower | ‘slope or smaller roughness, however, , the standing waves 
associated with the a jumps were no longer in evidence, and the flow 
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8 appeared to be in a tranquil state even over the tet of the roughness elements. 


the , and ‘the flow Classified into three 
- major regimes; (1) tranquil flow, which, essentially was maintained through- 
7 ‘ out, (2) tumbling flow—supercritical flow occurred over the roughness ele- 7 
ments with hydraulic jumps and subcritical flow in the regions between, and 
_ (3) rapid flow—supercritical flow occurred generally throughout the streams. i 
We In the flume, the flow could be readily classified into the principal regime 
by visual 1 inspection. ‘The depth Y; also helped distinguish the rapid and tran- 
 quil flow regimes because was smaller than Y¢ in the rapid flow regime 
and appreciably greater than Y, in the tranquil flow regime. No hydraulic _ 
jump between 1 roughness elements formed for e either the the 
% ims In planning the | program, the writers were quite e preoccupied with the tum- i 
bling flow phenomenon and, accordingly, most of the experimental runs were 
consciously chosen to assure that tumbling flow was induced. Nevertheless, 
_ ‘many runs occurred in the tranquil and rapid regimes. e 
sults, the latter are presented first, even though they are not as extensive as 


Tranquil and Rapid Flew Regimes. —For flows” in either of these regimes, 


: obtained were similar for both regimes. Fig. 4 shows a plot of ot K/Y, versus 
A Fi for tranquil flow over two-« -dimensional roughness elements. For this ‘Plot 
an all roughness elements were 0. 3 ft- -by y-0. .3 ft except for those points repre- 
_- by the triangles with the right half blacked, (S = 0.01, L = 15 K) that — 5 
_ were 0.135 ft-by-0.135 ft. The curves in Fig. 4 are typical in as to those a 
_ for various slopes in both the rapid and tranquil flow regimes and for both the 
-_two-and three-dimensional roughness elements. 
Study of Fig. 4 indicates that for intermediate values of roughness Spacing, 
_ the Froude number is smallest, and that for low values of relative roughness | 
the Froude number tends to increase. These observations are consistent with 
_ those made by Morris* in connection with pipe flow. Morris proposed all 


following classifications depending on n roughness spacing: 
1. Isolated roughness flow.—The elements are sufficiently 


so that the wake and vortex generating at element — 


2. Wake —The elements are cl e to 

3. Quasi- -smooth or Skimming Flow.—The roughness e elements | are re suffi 

- ciently close that the flow essentially skims the crests of the elements. _ 

aad The decrease in Froude number for the intermediate values of roughness 


ees could be associated with wake interference, whereas, the increase | 
associated wi ecrease eight and spacing of roughness could indicate 
a d with d d height and f hness could ind 
appears to as is by the for S = 0.01 


— 

7 = swept out by the increased momentum, and the flow appeared to skim over the [i 
© 
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@ L= 5K, S=0.0246 > 
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| 0.30 ftx0. ft bar agen: 
© 0.135 ftx 0.135 ft bar roughness 
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15 in Figs. 4 and 5 similar 


“aan natural eof trangull, quasi are of qt quite high intensity 
their heights are estimated to be of the order of 0.15 Yj. The velocity is 
probably approaching, but is less than i= 1. Fig. 7 shows natural rapid flow, 

a also probably in the skimming region. The flow of Fig. . 7 differs from  . 

Fig. 1 inthat essentially it is completely supercritical; no hydraulic jumps are 

a toform between the roughness elements, and the magnitude and direction 

of velocity ‘remains essentially the same over the roughness elements and 
tween them. 8 shows tranquil and rapid flow in the flume. 
_ Tumbling Flow.—The flow in this regime was characterized by alternate 
eres and deceleration from subcritical to supercritical velocities, and 

vice-versa, in a cyclic order. spacing between two successive e rough- 

Hydraulic jumps formed between successive roughness 

gh shows typical flow patterns in this regime for the. two- and the three- 

dimensional roughness" ‘elements installed in the flume. 1 illustrates 

a Based on considerations from dimensional analysis, a number of correla- 

tions were tried. _ Whereas ‘some appeared to be ‘successful, the simple rela- 


between critical “depth Y= 3 q*/g and ‘was most 


flow characteristics were independent of K/L 

a As in Figs. 1 and 9, it can be seen that each roughness element acts as an 

overfall control. : In the two dimensional studies, the control extends entirely 3 ; 


245 across the channel. For the three-dimensional arrangements, the extent of 
the control depends on the configuration of the roughness pattern and the lat- 
eral spacing of the roughness elements. | In the latter case, the stream t tends a 
_ bo flow at subcritical velocities between the roughness elements and at super a 
; critical velocities over them. The flow, however, is carried by momentum 
¥ both over and around the roughness elements and consequent acceleration 
causes a region of supercritical flow much larger than that actually occupied — ; 
f by the roughness (Fig. 9). A region of supercritical flow induced by a solitary _ 
_ rock in Logan River is shown in Fig.10. The rock is only about one-third the 
— width of the area disturbed. _ One would expect that a relatively low intensity 
a three- -dimensional roughness elements would induce the equivalent of full BS 
ss overfall ‘control in a rough reach of stream under tumbling flow. Thus, the 
y ip discharge under conditions of fully developed tumbling flow would be inde- a 
pendent of K, Land S. For this case, the depth Y, over the tops of the rough- 
ie ness elements would be less than critical depth computed on the basis of aver-_ f 


4 
age that | is, Ye = Va"/e = +c. For the ‘regime, therefore, the 
3 


Be, formula tor the critical flow weir, q = C g Y,° ‘should form the basis for the 


-depth relationship in which and are constants ‘depending on 
_— Fig. 11 is a plot showing Y, = Yq /g versus Yy for the 0.3 ft-by-0.3 ft 
plots for the 0.135 ft-by-0.135 ft bar roughness 
From this figure C = = 83 andc =0, © 
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‘fallweir for a reach of tur“ ling flow. The single overfall for natural vob 
— instead of being located at a particular cross section, is now divided into many 
rs smaller overfalls statistically distributed throughout the reach. Such a reach 
be regarded as a statistical critical depth weir for which the discharge 
fae depend on the relative size, distribution, shape, and intensity of 
the roughness configuration. The depth- discharge relation depends on neither > 
the overfall height for a weir, nor on the — and slope of a es flow 


10, EFFECT OF j 
ISOLATED 


ROUGHNESS 


of velocity and depth, that are now statistical. . This foregoing concept of a : 
virtual weir may furnish a means of approaching the problem. In such a weir, 
_ the crest elevation might be represented bythe proper statistical function de- — 
; - scribing roughness distribution | for a menting: reach of stream corrected for 
separation and roughness effects. 


ata minimum of _ 


in watch Z y Vy are, resp ctively, the height of roughnes s, 


be depth of tad and velocity ata vertical ‘section of the virtual weir, and y i = 
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a y the unit weight. . P is ; the power and h is the elevation of the energy line. For — 
on _ a@minimum, dp/dh = 0. The writers were unsuccessful in using this approach © 
_ since V, is an unknown function of x. . However, it now appears that statisti- — 
cally within the reach q, = Q/w = Vx y dx may be a valid assumption. “Cane 
‘machine computations this would make evaluation of the integral possible. | 
‘The? limited data available need to be supplemented by flume studies of flow 
over re roughnesses having statistical distribution of size and spacing and by ex- 
tensive study of natural streams in order to develop a generalized concept 
a flow and to infer the best means for describing it quantitatively. = 
s Instability of | Flow.— —Roll waves occurred in the flume under certain con- Bi 
ditions; Fig. 14 shows their occurrence in the flume. The leading crest of the Te. 
_ wave is at A. The trough is at B, building up to a crest at C. a 184) ule aact 
ag Ka Roll waves have been reported in natural and artificial channels by a num- 4 
ber of writers, ‘and several criteria have been proposed to ) predict the condi- _ 
_ tions under which they willappear. Depending on whether Chezy’s or Manning’s 
‘formula was used, mast frequently criterion is a lower | limit 


La 


14.—ROLL WAVES IN NFLUME 


_ Froude’s number of 1.5 or 2. .0, respectively, implying that flow above this 

aaa value of Froude’s s number would be inherently unstable, and that roll ‘waves 

ER would develop. Koloseus,° attempting to connect the Darcy-Weisbach resis- 

tance coefficient with free-surface instability, concluded from experimental 
ay studies that for F < 1.6, the flow was stable and the resistance | coefficient in- vag : 
dependent of Froude’s number. F 1.6, the resistance coefficient was 

dependent on Froude’s “number, the flow was unstable, and roll waves would a 

form if the channel were sufficiently long, == 
In contrast with the results reported in the literature, the writers didnot 

a roll waves in rapid flow even though the Froude number frequently - 
exceeded 1.6 and the flume should have been of adequate length. Rather, roll — 


a thesis presented to the State Univ., of Iowa, at Iowa City, in 1958, in partial fuifiiment 
of the the degree of Doctor | of Philosophy. 
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waves | formed in the transition ranges between tranquil and tumbling flow or a 
between tumbling and rapid flow. _ As reported by others, they formed only a. 4 ; 
_ when the depth of flow was relatively shallow. No roll waves were observed | 
using the cube roughness elements . Ordinarily, the large roughness elements 
out instabilities, and a stable condition was quickly established. How- 
_ ever, in the transition region a delicate condition of stability evidently exists. a “3 
= ‘During the present investigation, data on roll waves were collected only > 
incidentally are, thus, too limited for systematic study. ‘Two interesting 
observations were made: (1) For L/K = 5 roll waves occurred during the ~ 
transition from tumbling to rapid flow only when F,/VS = = 6, in which S a’ : 


- the channel slope, (2) For L/K = 10 or 15 roll waves occurred during > 


transition from tranquil to tumbling flow when F, = 0. 9. 
q _ Air Entrainment.—The writers had anticipated that air entrainment might 
create a problem, especially with regard to measurements. Except for the — aa 
ack case of three- dimensional roughness under extreme conditions of slope and , 
4g discharge, however, air entrainment did not affect the measurements signifi- _ 


os cantly. For both grrr tumbling flow, the e approach ivelocity to the — ie 


variation in depth, even though small. Where one draws the line is thus - 


smoot channels at low because any irregularities would cause some 


matter of judgment. For example, the flow shown in Fig. 6 would ordinarily | 
; Z be classed as uniform. However, what about the flow in Figs. 1, 7, and 9? 
; In the flume | Studies under a regular roughness pattern, the water surface _ 
ina regular repetitive cyclic way throughout the length of the flume, 
once equilibrium was established. _ One might, therefore, regard such flow as 
; . a 2 uniform, at least in a macroscopic sense. Over a reach of natural channel one © 
‘ _ would suspect also that there would exist a statistical uniformity of flow, again ey 
in a macroscopic sense. Perhaps a more general concept of uniform flow 
am _ would be more appropriate based on equilibrium between the weight component 
4 _ and the boundary drag (shear plus form drag of roughness elements) over a 
_ - long reach of channel so that statistical description would be con- 
stant. Under this definition, except for the runs where roll waves formed, : 
macroscopically uniform flow was quickly ett f 
Delineation of Flow Regimes .—The relation between 
a “height and spacing of roughness elements determined the flow regime inthe _— 
a . This is, undoubtedly, also” true in ‘nature, although the description of 
"the roughness may be quite complex. Since tumbling flow separates the re- 


which 1 this can occur should lead to criteria delineating the ‘principal regimes. 


‘The distinguishing characteristic of the tumbling regime is the occurrence 


pep of stationary hydraulic jumps. Thus, an understanding of the conditions leading 
we to the formation of hydraulic jumps in sloping gro should be 
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sy + W Sin 6 = 


dy = (11) 


This is plotted in Fig. 


d Gai For hydraulic jumps to be maintained, slope S must be sufficient to accel- a 
erate all or part of nee ory to — velocities, and the height of the 


| 
——+— 
a ‘ 
7 


per 
Poona 


Fic. 15. HYDRAULIC JUMP IN SLOPING CHANNELS 
-roughnesses must be sufficient to induce the required sequent depth. The slope 


on n the ‘channel bottom may be overcome 2 by the component of the weight paral- 

lel to the bed. Form drag depends on the flow distribution around the rough- oy 

ness and, thus, on the flow regime. From rational considerations together » with : 
36 empirical | evidence, tentative criteria delineating the various flow regimes in y 
the flume are proposed (Table 1). Analysis leading to these criteria are in- 


x. “Hydraulic Design of Stilling Basins: Stilling Basin with Sloping Apron (Basin — 
Bradley and A, J. Peterka, Proceedings ASCE, October, 83, No, 10, 
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Hydraulic jumps i Sluping chalneis lave Deen studied Dy a Mumper Of ve 
searchers, and the rest 4 
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rapid 
tranquil or 


=p 


"Skimming 
interference 
or isolated roughness 


cubes 

lz + /( @ 16 j tumbling or tranquil 
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wt will 


paratively low relative ‘roughness, use in the past has been extended | to 


4 z pass through the critical value regardless of S, L, and K. Thus, in Eq. 12 the 
s - Values of C', m, and in | would tend to become 1, 1/2, and 0, respectively. Ih 


a Within the transition zones, ‘the values o of ,m, , and n apparently will vary 
pending on of the flow balance to ‘tumbling flow. 


pressed in the general — 
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FLUME STI STUDIES 


of and L/Y1, discussed previously, 
suggests that the important parameter representing roughness geometry 


- In this form, A is the ratio of the projected effective area of the sidealen a 
«perpendicular to the direction of flow ; to the total | area . of tne channel bed. a: 
It is believed that both and n are of the relative roughness 


evaluation of Eq. 16 and 17 for even the simplified conditions of the flume 


: he Roughness in Hydraulic Design.. —The writers believe that use of roughness a 
elements have not been fully exploited in hydraulic design. With better under- — 7 
standing, their rational use could very possibly result in improved perform- 

_ ance and economy of certain structures. Some e specific possibilities might ae 7 


Energy Di Dissipator. r.—A chute with a single large stilling basin is fre- 
quently very expensive. _ By proper use of roughness elements it might be 
practical to establish tumbling flow throughout the entire length of the chute. ra 
The energy would then be dissipated by a number of overfall sinks throughout 

_ its length rather than by a Single large ge one. Velocities would never ‘greatly 
Use of Steeper Slopes in Channel Design. —In earthen steep 

slopes resulting in scouring velocities are usually undesirable. In lined chan-— 

nels of steep slope there is the danger of instability, impact waves, , and forma- 4 
ton of roll waves. Use of roughness elements to control the flow at a suitable 
p> velocity might permit the use of steeper slopes where this is desirable. 
Increased Sediment Transport. —By using to induce greater 


However, lack of adequate data in all sub-regimes of flow prevents quantitative, =} ‘e 


flume to jo develop 8 sufficient turbulence to ee the bed load ofa stream, thus &§ 


“Total Sediment Load Measured in Turbulence Flume,” by P. C. Benedict, 
and D. Q. _ Matejka, Transactions ASCE, ' Vol. 120, 1955, PP. epibess 


— _ Due to the extre as compared to its sides, the [i 
_ channel should behave essentially as avery wide Channel, and m should be very “4 
— pe 
to 1/2. This strongly suggests the possibility that a general open channel 
— 
— 
1) 


making p possible the measurement of total transport. This same idea could 7 ” 


3 utilized to insure transport of sediment loads through channels. ARIS Fie 
; 4. Stream Measurement. — —Through a San understanding of the role le of 


Rayan Also, because | uniform flow is rapidly induced by large roughness 
elements, a ‘roughness flume might be developed for water measurement. This 

19 _ would seem to be a promising device for intermittent streams transporting e 
large quantities of sediment. A properly designed tumbling - -flow flume might 


‘permit use of a flume having ‘approximately the same natural slope and rough- 
ness as the stream, thus gay the problems inherent in providing a — 


open channels comparable to that of Nikaradse for is desirable, but 
motbeenfullydeveloped. 
6. Flow Analysis. Understanding of the ‘around 
roughnesses in open channels under various flow conditions dy- 


} 


ms. 


equation ¢ of the type = C'S", in ‘which and n depend on the roughness 
In the tumbling regime, the jumps form over the. 
_ elements to the degree that a reach of stream acts as a control. Such a reach 
j. may be regarded as analogous to a critical flow weir. 
s = 4. For shallow depths under certain conditions, roll waves | may form a’ at 
_ the transition zone between tranquil and tumbling flow and rapid flow. ie 
= a * . Based on better understanding of their action, the use of large rough- 


, ness elements could result in improved hydraulic design for energy dissipators, 


as part of 2 a continuing project on flow in steep, channels. 
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the data, this was true within 2%. Considering 1 ft of 


applying the principle of “ie 
in n which. FF, is the pressure ‘on the downstream face of the ceeds, ‘ 


For _supereritical values: of , if the right side 21 is greater than ‘the: 


/ +23 2 ¥) 


‘ 
= 0. 3 Cp = 1.16, approximately, 
(% + / (0.35 +2. ¥2) 
sie, 
(F, exceeds ‘of 1/2 or of (yy + Y2) / (0. 35 + 
rapid flow would be be expected; not, thar or flow v would» 


the sequent depth ‘exceeds K + Yor rapid | flow pre- 
| Fy dy (S + 0.153) /0.183 > 


Spacing.— In Fig. 16(d), TT? distance to the. at which the 
centerline of the jet intersects the line joining the tops of the roughness ele- 
= it L hes 1, the flow will skim. Let x and Zz be the horizontal and vertical | 
— —Keutner8 found that for H 0. 85, Fig. 1646), the nappe 
a/ 
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hren aus der Geschw — 
Berchnungsverfahren fur by C. Keutner, Die Bautechnik, 


FLUME STUDIES 


"Length of Jump.—If d is the uniform flow depth established en channel 
excluding the roughness elements, and da is by Ye + K, 100 


2/2+pavy + W Sin = 
me: 28 may be used to compute W. Assuming an average a ani ofjumptobe §f 
(dy - + te + K)/ 2 (or any other appropriate average depth), compute the length P 
of jump 1 j necessary to provide the weight. If this length exceeds the ‘rough- 
Residual Slope. —Part of the channel slope may be thought of as providing 
gravity force to overcome the drag offered by the “roughness element. - 
si this be designated Sp. If S - Sp < S,,; in which S,, denotes critical slope, the 
velocity cannot reach critical, and tranquil flow will 


Roughness Height. —Let the velocity of approach, Fig. 16(c) be V, : 


and rearranged toform 


= 
‘Values of K1 less than that given by Eq. 29 will be too ‘small to induce mnie: ; 


‘The following letter symbols, aiiaianies use ir in the paper and for the guid-— = 


of discussers, conform with American Standard Letter Sym- 


/ 
upstream head over the weir; 


= downstream head over ‘the weir; 


aaa + = transverse width of roughness per unit width of channel; 
= height of roughness; 
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= 1960 : 
= number of ‘spacings unit wi 
channel; 


= slope. used to gravity overcome form 


drag; 


width ‘of flume or channel; 


vertical distance; 


shear; 
= angie of 
= kinematic v 


7 
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— = discharge per unit width of c : wae 
ity at upstream roughness crest; — 
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"TESTS ON CONCRETE EMBEDDED CYLINDER PIPE 


2g Hugh F. Kennison,! M. ASCE 


=a of stress cracks in the concrete and on strain recordings from 


te INTRODUCTION 


introduction in 1942. In this s original p pipe, the sueaknennall wires ance st wrapped 
- directly on the steel cylinder encasing the concrete core (Fig. 1). “The pipe is 
"manufactured in diameters from 16 in. to 48 in. In 1952 the first prestressed J 


i “concrete embedded cylinder pipe was commercially produced. In this pipe, as 


distinguished from | prestressed concrete cylinder pipe, the prestressed wires _ 
are wrapped about the concrete core in which the steel cylinder is pene be 
- (Fig. 2). This pipe has been manufactured in diameters from 24 in. to 120 in. 
and is particularly applicable for large diameter ; and high pressure pipelines. 


} a Both types of pipe are described in the American Water Works Association , 
Prestressed concrete embedded cylinder pipe consists of a steel cylinder 

P ed with bell and spigot joint rings | welded to its ends. This steel membrane is Ml 


a... hie! Note.—Discussion open until April 1, 1961. To extend the closing date one month, a a 
7 . a written request must be filed with the Executive Secretary, ASCE. This paper is part 
= 4 of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
2 ciety of Civil Engineers, Vol, 86, No. HY 9, November, 1960, a 


1 Vice Pres., Engrg. and Research, Co., Bast mt Orange, N. 
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MORTAR , PLACED AFTER INSTALLATION PIPE DIAMETER~ 


catia FIG. _PRESTRESSED CONCRETE CYLINDER PIPE 


STEEL CYLINDER STEEL BELL RING 


NOMINAL PIPE 
DIAMETER— 


LINDER PIPE 
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_ After the core has been properly cured, it is helically wrapped with ieailhohaien ia E 


_ wire under a high uniform tension. The wire is then protected with a concrete = q 
To obtain realistic information onthe physical behavior of prestressed con- — 
crete embedded cylinder pipe, the American Concrete Pressure Pipe . Associa- _ 
tion (ACPPA) sponsored a series of tests ata laboratory in Wharton, N. J. 
This three- -day test program was witnessed by consulting engineers and major 
Consumers from New York, N. Y., Philadelphia, Pa. -» Baltimore, Kansas 
_ City, Mo., Dallas, Tex., Fort Worth, Tex., Ames, Ia., Denver, Col., Los Angeles, — 
Calif., and San Prancisco, Calif. “These tests were designed to termine the 
a structural behavior of 60-in. pipe when subjected to combinations of external 
‘o - loading and internal pressure, since under normal use these simultaneous con- 


pare PREVIOUS 1 TESTS CONCRETE CYLINDER PIPE” 


ao pipe have indicated a cubic parabola relationship between three- — “2 
bearing load and internal hydrostatic pressure, as shown in Fig.3. 
The basic curves resulting from these earlier investigations are the “design” 
- curve, and the “ultimate” curve. The “design” curve has been the principal : 

control for a combinedanalysis of prestressed concrete (lined) cylinder pipe 4S : 


4 The combination of normal operating pressure and dead load has been Sian 


ie 


within a curve which is 90% of the 
a a, 2 The “ultimate” curve represents the combinations of pressure and load which ; a 


‘consequently seldom used. Under the three- -edge type of loading this curve docs 
a not give an accurate representation of the pipe’s performance under field con- be 
he eer . As the pipe nears its ultimate capacity it behaves as flexible pipe,and 
under field conditions this capacity | would be far greater than indicated, due to - 
the extra side support given by passive earth pressures. In actuality, a much —- 
higher bedding factor should be with the ultimate curve, to obtain real- 


ae The terminal ordinate and abscissa values | of these basic peri can be cal- 
culated and/or verified by relatively simple tests. These terminal values are: a 
_ 001 in. = three-edge bearing load producing the first 0.001-in. crack(in- _ 
 cipient cracking) ; = =  three- edge bearing load producing ultimate failure 
of the pipe; Po = pressure causing zero concrete ices in the core; and Pp 
met DESIGN FOR 60-IN. LOAD TESTS 
“3 for this | test since it falls approximately in the middle of 
the range of sizes and pressures for this type of construction. Seven identical 
test specimens were constructed and tested at various combinations of inter- 
nal pressure and external loads. These test pipe were 4 4 ft long, 60- in. inter- _ 
nal diameter designed for a combination a pressure of 160 


___hydrostatically tested for water tig 
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FIG, 3. —~COMBINED LOAD CURVES FOR.  PRESTRESSED CONCRETE 


Ultimate. 


ae 
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BEARING 


FIG. 4. COMBINED LOAD DIAGRAM FOR TEST PIPE <2y 
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_ PRESTRESSED PIPE 


and a three- -edge load of 6, 175 lb per ‘lin ft. common design was 
- chosen for all specimens to minimize the many variables which occur in such 


‘concrete coating in.; cylinder = #16 U.S.S. hot Tolled 
“¥ (ASTM A245); prestress wire = #6, MB untempered (ASTM A227); steel cyl- ‘ 

- inder area = 0.718 sq in. per lin ft; steel cylinder diameter = 63 in. (outside 

diameter); prestress wire spacing = 0.392 in., c. toc.; prestress wire area 

_ 0.883 sq in. per lin ft; Py = normal operating pressure = 160 psi; Po 1‘ zero aa 

_ compression pressure = 216 psi; P) = bursting pressure - = 535 psi; Wy = al- 

_ lowable three-edge bearing load at operating pressure = 6,175 lb per Yin ft; 
Walt = load to cause three- -edge | bearing failure at zero pressure = 27, 200 lb p : 
ss lin 73 ; gross wire w wrapping stress = = 140,000 psi; resultant wire 

= 105,000 psi; and wire stress at P,, = 111, 
it 4s shows the design curve, the 0.9 design curve, the ultimate curve the q 
ormal operating pressure, and the normal three- edge bearing load for these Pa g 


po specimens. These curves will be reproduced on all ees aan of 


te 
the test results for reference purposes. «la ae 


For these tests the ASTM standard three- edge bearing load was , chosen as 
- the method of applying the external load because it was felt that this a 


TABLE | 1.— ~THREE- -EDGE BEARING LOADS TO 0 EQUIVALENT 


DEPTHS OF COVER fi 


Bearing Load, in 

Foot = 100 y= 100 pot é w= 120 pet 


with Bedding, i in Fet 


This table of equivalent depths of co cover ris based on on the studies of Marston, ¢ et al, 
: assumed “ordinary” bedding (a bedding factor of 1.5 when compared with three-edge — 


soil constants of KH and of 0,130, the ordinary maximum for 
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__tensive research work — 
oreviouslv performed at Iowa e College (Ame 
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interpretation mn of relatively simple and ‘routine crushing tests. nile 
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5.-SYSTEM OF 6.—60-IN, PRESTRESSED pet 


7. SCHEMATIC VIEW OF PIPE SPECIMEN ASSEMBLED 

_ been prepared for converting | recorded three- edge bearing te test 

= to the more realistic trench conditions, in feet of cover over the pipe; 


in Table 1 ‘Bg is the trench width, Do Canctes the outside diameter of the Pipe, 


4 Fring the three-edge bea additional 
informatio ions. The use of this excellent material elimi 
iti 
— 
— 
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the ‘combined sof ‘three- edge bearing external load and internal 
hydrostatic pressure under which these tests were carried out is illustrated 

os in Fig. 5 . The 60-in. test pipes, , with 18-in. long thin steel cylinder extensions, 
were placed over a tubular internal bulkhead. The steel cylinder extensions 
° eo reduced the supporting effect of the bulkhead on the three- -edge bearing load. J 
The test specimen completely assembled in the combined i load testing facility 
— is shown in Fig.6. The arrangement of the test specimen mounted in the tubu- “a 4 
ane lar bulkhead is illustrated in Fig. 7. The support which was provided by the — 

hydrostatic bulkhead was carefully det determined all three- -edge b bearing loads 
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TERIOR VIEW OF TUBU-_ 
LAR BULKHEAD 
THICK GLASS WINDOW, 

SUBMARINE LIGHT HOUS- 
AND LONG FOCAL 
LENGTH MICROSCOPE 


this “The tests and results for evaluating this 
mae! _ Throughout the tests it was necessary to observe the inside concrete lining as 
of the pipe to detect and measure the development of lining cracks. To ac- 
-complish this the tubular bulkhead was constructed with a thick glass window, 
through which the inside lining of the pipe could be surveyed with a 40X micro- hae: 
was placed opposite the bottom of the pipe, which would 
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marine lighting device was constructed within the bulkhead to illuminate the 
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of the glass window, the microscope, and the submarine light is 
PROCEDURE USED TO OBSERVE TEST RESULTS 


One of the major indications of adequacy of design in prestressed concrete, 
; al in reinforced concrete, is the degree to which stress cracking of the con- ‘a 
crete is controlled under load. Reinforced concrete pipe must necessarily 
a = tensile stresses and resulting — of the concrete whenever the ten- — 


b) 
margin exists. between the load or pressure 
¥ a causing first visible stress crack , and that causing ultimate failure. a 


10 (a), 


= 
th concrete pipe, the stress crack occurs in the exterior 
Roe coating and is conservatively defined as having a a sustained maximum width in 
‘ ‘pix These tests report on the cracking of the lining and eonting caused by the 
"combined load application. The cracking behavior of the pipe was obtained - 
a careful examination of the concrete lining and coating with a 40X microscope *_ ar 
- equipped with a calibrated grid to accurately | measure crack widths. The con- 
of cracking or first tensile failure can recorded by trained 
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visible to the eye a and is defined as ac crack, 0. 001 width 

ek: and 12 in. long. The first visible crack is defined as acrack 0.002 in. wide and sd] 

al 12 in. long and can usually be observed with normal eyesight if a meticulous - 

a ae of the surface is made. The first visible crack supplies a consis- _ , 
+ ~y tent criterion of the stress cracking of the concrete and is so used in the ngs 

- cording and interpretation of the results of these tests. Fig. 9 shows the sur- 

vey of the pipe coating by microscope to detect and measure stress } exacts 4 
occurring in the coating while the pipe is under a combined load. 

Ls further study the structural behavior of the pipe ‘under combined loads, 7 

strain recordings were obtained through the use of SR-4 gages. These gages 

were installed in three test specimens at various locations as shown in Fig. - 

10(a). ‘During the tests the pipe were placed with the strain gages at specified 

- locations (Fig. 10(b)) relative to the points of application of the external load. Pe 


Strain readings were then taken during load applications. 


a.) seven test specimens were numbered 1834 through 1840 for purpose of 
record. Table 2 gives the history of these seven test 


‘AND CONCRETE 


Tested 4 Conerete Strengths, in 


Core 


ane #1836 we was used in the raeti conducted to »evaluate the effect of the ses 


2500 lb per lin ft; 
were measured to the nearest 0.001 ‘in. at each of load. This 
_ was repeated to record any effect of concrete tensile failure. The lower curve = 
_ shown in Fig. 11 resulted from these two tests. The 18-in. long steelextension — 4 
: cylinders: were then welded to the ends of pipe #1836 and the unit assembled in 
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Date | Date | Date 
Cc] Prre- | Coaat-| bination | dro-| Test 
Cast stress | ed | Loading |staticl) — 
— issale/s | 7/9 | 1/29) 9/28 | 9/30) - 4,200}5,700| 4,100] 5,000 
6/14 | 7/23 | - | 11/7 | 4,150]5,400] 4,000| 5,180 
6/29} 7/9 | 8/3] - 9/9 & 9/13| 4,200] 5,100 | 4,400 
887 | 6/23 8/1 | 4,650) 5,625 | 4,000 5350 
7/9 | 8/2 | —-9/30 | 9/30 | 3,775|5,150 | 4,250] 5,050 
7/9 | 7/26 19/29 | 4,400| 5,750 | 3,800 4,800 
1840|6/16 | 7/9 | 8/5|  - {9/28 | 4.27515,150 4.050) 4.900, 
— 
aa 
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PRESTRESSED PIPE 
curve of Fig. 11. . The difference in three-edge bearing load between the two 


curves, for a given deflection, is then takenas the supporting effect of the hydro- 
static bulkhead at that load. Fig. 12 shows the resulting correction factor which 
used to obtain the true three-edge bearing load. 
TEST PROCEDURE AND RESULTS BASED ON CRACKING BEHAVIOR 
‘The following discussion contains a separate section on the | tests and re- 
sulting data, based on cracking behavior, for each of the individual test speci- 
_ mens. Results obtained through the use of the SR-4 strain gages will be pre- 


sented separately following this section. A variety of test procedures was 
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(1000 Ibs. /lin. ft) 
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—Two Additional “incipient” Cracks 


"First "Visible Crack 


wo Cracks Increased to.002" 

Crack-.0 


PRESSURE ( (Ibs. / 


4 


ag 


é 


incipient” cracking (0.001-in. x 12-in. crack); (3) subjecting the specimen 
at to combined loadings sufficient to produce the first visible crack (0.002-in. “te 
: 12-in. crack); (4) subjecting the specimento combined loads in which both the aa 
‘a pressure and three-edge load were increased by a predetermined 

(5) subjecting the specimen to an increasing three-edge bearing load 


= the a variety of combined loadings si sufficient to produce 


(no internal pressure) resulting in ultimate failure of the pipe; and (6) sub- 
a - me the specimen to an increasing internal pressure (no three- — bearing 


sulting on a combined load curve this design such as was shown in 
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Movember, 1900 


i: cause of ‘the ‘very adequate margin between the theoretical design, or incip- 
ient -erack, curve for the pipe tested and the actual pressure - load combina- 
= @€©—s«éPipe: #1834.—Pipe #1834 was first subjected to a series of combined loads — 
that traced the 0.9 design curve. These combined loads are | plotted on Fig. 13. 
It was next subjected toa series of combined loads that traced the design curve. 
_ No microscopic cracking of the lining or coating was detected during these two pt 
- series of combined loads. Following this the internal pressure was held at ‘es 
160 psi (normal operating pressure) and the three- -edge bearing load was in- 
et until the first “incipient” crack (0.001-in. x 12-in.) formed. This first be 
“incipient” crack appeared at the top quarter point of the spigot end. The last il 
- test consisted of recording the cracking behavior of the pipe as the ‘internal 4 
_ hydrostatic pressure (no three-edge bearing load) was increased to the ‘burst 
ing pressure. All loads applied to pipe #1834 are plotted in Fig. 13. he 
‘a fi #1835.— This test specimen was subjected to an increasing three- edge» 
ne load (no internal pressure) until it reached ultimate failure. The a 
+ ~ tical and horizontal diameter changes were recorded along with the behavior 
a the crack formation. The relationship between three-edge bearing load and 
the diameter changes are shown in Fig. 17. The various cracking loads and 
- information is plotted on the ordinate of the combined load curve in 7 
«Fig 14. The results of the hydrostatic tests on test specimen #1839 and #1840 
4 are also plotted in Fig. 14 as will be discussed subsequently. ait, 
baile Pipe #1836.—This test specimen was used to develop the correction factor — 
Fs to be applied to the three-edge bearing machine | load caused | by the supporting © : 
— effect of the hydrostatic bulkheads. This test was discussed and the results 
presented in Figs. 11 and 12. 
#1837.—This test specimen was first subjected to various combinations 
ig : of internal pressure and external three- -edge bearing loads which produced = 
} = first visible crack (0.002-in. x 12-in.) condition, _ The test procedure was to 


hold pre-established three-edge bearing loads and increase the internal pres- 
sure until the first visible crack occurred or reoccurred under subsequent loads. 
\ The combined loads resulting in the first visible crack are plotted on the com- 
bined load curve in Fig. 15. In this series the first visible crack occurred in ; 
a the pipe coating under all load combinations. The controlling crack under in- 
» _ ternal pressure and no external load occurred at the top of the pipe while the © 
controlling crack under external load and no internal pressure was located at © 
the springline (90° from the topofthe pipe), 
Next this pipe was subjected to three combinations of internal pressure and 
three-edge bearing loads in which these | loads were increased simultaneously 
to 2.0, 2.5, and 2.75 times the normal operating loads of 160 psi internal | pres- 
_ sure and 6175 Ib per lin ft three-edge bearing loads. The maximum crack 
with at each of these three combination loads was recorded. inf informa-— 


Pipe #1838.—As in pipe #1837 this specimen was ‘first to 

“combinations of internal pressure and three-edge bearing loads which produced % 
the first visible crack (0. 002- in, x 12-in, ) condition. The procedure on this” 
en 


__ inerease the three-edge otnrtin load until the first visible crack formed. ‘The 

: ~ combined loads resulting in the first visible crack are plotted onthe combined _ 
he load curve in Fig. 16. The controlling crack, in the case of pipe #1838, occur- 
in coating and started at the upper quarter from the of 
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Next, as in pipe #1837, this specimen was subjected to three combina- 
ions of internal pressure and three-edge bearing loads in which these loads 
were increased simultaneously to 2.0, 2.5, and 2.75 times the normal operating 
loads of 160 psi internal pressure and 6175 Ib per linft ft three- -edge bearing loa 
_ The maximum crack width at each of these loads was recorded. This inform a 
tion is also recorded on the combined load curves of Fig. 16. ‘This specimen ai 
Ne was then subjected to a series of increasing internal hydrostatic pressures — 
to 300 psi for the purpose of strain recordings, the results of which will be agi 2 


19. FAILURE OF PIPE NO, 1840 UNDER HYDROSTATIC 


TEST PRESSURE OF 605 PSI 


final test on pipe #1838 consisted of increasing the three- -edge bea ing 
oad (no internal pressure) to ultimate failure and recording the cracking be- _ ; 
havior throughout. This information is likewise recorded in Fig. 16. Previous | 

to this last test to ultimate failure this | pipe had been subjected to numerous ees g 

= high combination loads including one above the theoretical ultimate i 

of the pipe. This resultedin a substantial increase in the flexibility of the pipe Bs 
tie as is shown by the perecene diameter deflection in Fig. 18. In spite of this the — # 


Pipes #1839 and #1840.—These twotest specimens were subjected to an in- -— a 
creasing internal hydrostatic pressure (no three-edge bearing load) upto burst-— 
pressure of 620 and 605 psi respectively. The cracking behavior of the 
Reeossei was observed and is shown on the combined load curves 3 of Fig. 14 _ 
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Bee 
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along with the results of the three- -edge bearing load —_ on pipe #1835, Fig. 


7 a Pipe 2 #1834. _—Pipe #1834 was constructed with SR-4 strain gages distributed 
‘ oe radially through the walls of the pipe at three different circumferential sections. 
The radial distribution of the SR- 4 gages is shown in Fig. 10(a) and the 
of the circumferential sections is shown in Fig. 10(b). The circumferential 
ha sections were spaced at 45° so that the strains recorded during the combined 


_ loads were those occurring at; (1) the top of the pipe directly under the three a 


edge bearing loads; (2) the quarter point, 45° from the top of the pipe; and vedi 
the springline or 90° fromthetopofthe pipe. 
at ‘The strains occurring at these locations in pipe #1834 were recorded during 
the combined load applications which traced the 0.9 design curve and the de- _ 
sign curve as shown previously in Fig. 13. Since the strain recordings in the 
) re two series of combined loads were of essentially the same pattern but slightly — 
d 


sign ¢ curve will be presented. . Some: of the gages v were damaged during con-— 
"struction of the test specimen or for some other reason failed to record properly — 
during the test. For this reason there is not a record for eachof the SR-4 gages 
shown in ‘Fig. 10(a). 20 shows a plot of the strains recorded during the 1 
3 application of the combined loads on pipe #1834 which traced the design curve. 
i: The letters on the abscissa of the strain curves denote the particular load | ap- 
plication, as shown on the combined load curve. 
The solid lines in Fig. 20 represent the strains as. recorded by the SR-4_ 
gages. The dotted lines are computed values of the strain obtained froma theo- 
retical analysis of that particular section under a combined internal pressure 
external three- -edge bearing load. ‘For the computed curves an elastic’ 
modulus of 28.5 x 106 psi was used for the steel and 3.7 x 106 psi for concrete. 
_ As seen on the curves of Fig. 20 the computed and recorded strains in general = 
follow the same pattern, The displacement of the curves representing gage #7 
at the top of the pipe is difficult to explain unless it is a zero shift of the re- 
Corder that occurs only when the pipe is loaded. The larger variation between rs 


a _ Pipe #1837.—Pipe #1837 had SR-4 strain gages mounted on the prestress wire 
and the concrete coating at the springline on both sides of the pipe. There were 


_ two gages on the prestress wire and two on the concrete coating at each location. 
_ Type A-7 SR-4 gages were used for mounting on the coating and Type A- 9 for 


2 produce the first visible crack as was previously ‘shown in. Fig. 15. The : strains — < 
- resulting from this series of combined loads is shown in Fig. 21. The strains _ 
shown in the upper diagram of Fig. 21 are plotted directly above the combined _ 
_ loads, plotted in the lower diagram, which produced those strains. Since ~ 
; "aye was reasonable agreement between the four SR-4 gages on the prestress wire. 
the average strain recorded by these gages has been plotted. The average of 
coating on each side of the ave been plotted separately. 
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a ae shift in the zero readings of the coating gages is shown by the © 
= measured strain following removal of the last load (load K). 
_ This strain data was recorded when the first visible (0. 002-in, in, x 12-in n.) 
i 2 ‘crack formed, which, in all cases occurred in the coating. Since any Beet. 
a — concrete coating will necessarily relieve some of the tensile stress in 
‘4 a. the coating, the strains as recorded by the coating gages would be affected ted by 
location and general characteristics of the coating cracks. 
#1 #1838.- - —This test specimen was instrumented with strain gage 
a section. Again, as with pipe #1837, strains were recorded during the ap- 
plication of the combined loads to produce the first visible (0. 002- “in, x 12-1 in, Le) 


4 during this series of tests are shown in Fig. 22, “Again the recorded = 


meme —HYDROSTATIC LOAD STRAIN DATA FOR PIPE NO. 1838 : 
a Pion e the four SR- -4 gages on the prestress wire have been averaged and 


is nearly the same as ‘that exhibited by pipe #1837, However, the indicated 
strains of the concrete coating on pipe #1838 are considerably greater thanon 
a pipe #1837. As discussed in the previous section the strains measured by the 

o coating gages could be affected by the cracking characteristics of the coating. © 
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November, 1960 
both cases #1837 and #1838) the strain in the remains sub- 
stantially constant while the prestress wire strain decreases with a — 
Following the increasing combination loads on pipe #1838 (shown in Fig. 16) > 
strain recordings were taken onthis specimenas the internal hydrostatic pres- 4 
- gure was increased from 0 to 300 psi in 30 psi increments with no three- alee. 7 ia 
: bearing load. The resulting strains are plotted in Fig. 23. Since there was 7 | 
agreement between th the various gages, the average, of the four prestress. 
_ wire gages and the average of the four coating ggaes are plotted. Shown in Fig. a 
23 also are; (1) the theoretical strain - pressure relationship for the prestress iy 
_ wire with no coating effect with and without all the inelastic losses occurring 7 i 
> after the core is prestressed; : (2) the theoretical strain pressure relationship | i> 
oF for steel only. The displacement of the coating strain curve to the left of the , 
prestress wire strain curve can possibly be explained by the stress relief in = 
. the concrete coating caused by the cracking of the coating . It must | be remem- 
_ bered that previous to ) this test pipe #1838 had been — to numerous load — 


The prestressed concrete embedded pipe in this” program 
were designed for r operation at 1 160 psi internal pressure and 6175 Ib per lin ft 
-three- -edge bearing load. This r represents an equivalent trench condition of 
ft of cover with ordinary bedding in a trench whose width is 3 ft greater than — 
a the outside diameter of the pipe. These normal operating conditions fall on the = 
0. 9 design curve which is 10% below the design curve, i 
_ The test results illustrate the conservative design of - prestressed concrete — 
embedded cylinder pipe, and clearly demonstrate the adequate stress crack _ 
- control provided by the present design. No stress cracking was observed on 
the lining or coating for any combined load falling inside the design curve. Ino 
fact as pointed out under the heading “Test Procedure and Results Based on 
_ Cracking Behavior” there is a considerable margin of safety in the oe 
- erack condition. A combined load of 400 psi internal pressure and 15,400 lb 
B per lin ft three-edge bearing load was required to produce the first critical © .¥ 
= (0, 005-in. ) in the concrete. This is a combined load equal to 2.5 times 
the design operating conditions for both internal pressure and external ca 
Upon removal of the excessive loads, these cracks closed to the point where 
only visible the 40X n microscope. This illustrates that minor 
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temporarily exceed normal design conditions is not serious, since the q 4 

@ stress crack will close when the excess load is removed. | 
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“Ina an attempt to bring to hydraulic engineers an up-to- listing trans - 


pre- 


Pit = suitable items. catewebtel readers are urged to ‘submit discussions to (a) add 
. ss this list, (b) offer suggestions for ‘improvement, and (c) otherwise assist the 
at task force in fulfiling its aims. The present list is the third issued since the a 
_ publication, in 1957, of ASCE Manual 35, entitled “A List of Translations of 
_ Foreign Literature on Hydraulics.” The first list appeared! in the August 1959 — 
_ Journal of the Hydraulics Division and the second in the June 1960 journal. - 
This list is to be consideredas anaddendum to Manual 35 and will be included 
_ in a revision of that Manual when the time is deemed appropriate. lie ae ae TN 
fallsverk pp. 75- 76, 1951. (On file at Army Map | Service. ) Ref. ‘Ro. ‘TC 89, a 


Experiments in Hydraulic Engineering.” Statens Vatten- 
- fallsverk, pp. 73- 74, 1952. (On file at Army Map Service.) Ref. No. TC 89, a : 


ANONYMOUS. Model in Engineering.” Statens Vatten- 
= fallsverk, pp. 80-81, 1953. (On file at Army Map Service.) Ref. No. >. TC 89 _ 
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“ the "individual papers in this symposium. To extend the closing date one month, a 
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— pp. 46- 48. Tr. by Technical Information and Library Services, Minstry 


ofS Supply, London, England. (ASTIA Order N No . AD 215819). ‘Saar 


CHECHEL, «The Effect of the of Blades on the Power 
- Output of High Speed Variable Pitch Blade Turbines.” Transactions of the 
Lenin rad Polytechnical Institute of M. I. Kalinin Ministr of Higher Educa- 
tion of the USSR, No. 193, pp. 38-50 » oe (On file at DSIR, LLU.) Ref. 
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[RESISTANCE EXPERIMENTS IN A TRIANGULAR CHANNEL! 


RALPH w. POWELL, 21 F. ASCEand CHESLEY J. POSEY,22 F. ASCE.—The 
authors wish to express their appreciation for the contributions of the discus- 
sers. They also wish totake this opportunity to report on additional tests which 
have been made since the original paper was prepared. Table 3 gives the re- 
‘sults of 19 more runs on the smooth channel and 27 runs on the channel with 
the walls coated with paraffin. This was a mineral wax with small — 


content, suchas is used for waterproofing paper milk cartons, which was heated 
and rolled onto the inside surface of the flume with painting ‘rollers. As the 


contact with the e relatively cool surface of the ~ 


= to 1 measure this roughness ‘and analyz ze the “meas- 


“urements, was built, but so far it has not been perfected. 
5 _ The numeration of runs in Table 3 is not continuous with that of Table 1. : 


pezoidal "The re: results for the trapezoidal shape will be reported else- 
one of the 19 runs on the smooth triangular hada Froude 
ae ber: over 0.85 (No. 202 with 0. 88), and it showed no greater resistance than the aa | 
others, so that these runs are all classed as tranquil flow. The average n by = ; 
eid Manning’s formula was 0.0092, the same as for the 23 runs already reported. _ 
The average discrepancy v was s slightly more, however, , bringing the oe a: Ly 
entire 42 runs with tranquil flow up to 3.8%. 


& 


Mr. Ackers ‘Says that for the channela better fit is obtained the 


is changed to 2.2. (He He Lon uses 14.8 instead of 14.83 and 2.51 instead of 2.52. ee 
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DATE TIME R 
92-58 11:00-11:23 PARAFFIN V 0490 12.0 -0153 0.334 ia er 
92:58 10:31-10:46 0-686 11 8 0163 0-439 0%9 1.89 0241-0106 
9-258  9:56-10:14 V 0.889 11-4 0153 0-560 1.852 1-97 1075 
9258 9:19-9:35 1-084 10.6 0153 0.694 3.04 1028 0243 0114 
9-2-8 8:46-9:06 V 1-236 10.0 0153 0.762 422 207 tae 0106 
92-8 14.5 0152 077% 4.55 214 63 


9-258 4:23-4:37 040 147 0152 1.006 6.18 102-8 0213-014 
92-58 357-413 


0.600 149 0152 1.165 12-27 2. 09 1207 1144 0196 
9-258 390-34 


(0-740 15.0 0152 1.261689 220 497 1205 0177 


184 I2B 2:59-3:16 14-9 0152 1.442272 230 1812 (25-9 
a 

V 0-349 14.6 006% 0-313 0.1818 0-83 59-6 
V 0-500 14-4 008% 0.405 0.42 1-06 1233 763 -OM2 


98-58 1102-11-28 123 0.683 2.32 1-50 1085 0108 


9-6-58 10:08-10:30 V 1-275 11-3 0-891 4.56 -1093 0215 
95-58 2:49- 313 40-499 13-7 0067 1.068 7-65 155 0208 
P55 1:55-232 0-623. 13-7 008% 1-33513-00 1-58 10% 113-5 -0200 
9-558 1051-11:35 4 0-836 13.0 0087 | 036 | 1194 0180 
73059 11:20-11:26 SMOOTH V 0-352 135 004 0.343 018% 067 59.8 119-4 0180 
7202) 10-41- 1:00 12-8 0516 0-556 073 1170 “0155-0087 
73059 9:39-10:05 00144 0.649 |. 049 0-77 137-0 -0137 
72959 2:52- 1-228 13-2 -00144 1-106 4.15 0-81 0126-0089 a 
747-59 0:20-3:28 0-349 12.1 0047 1-7 541 0.78 
7209 3:04-4:24 084 122 2-01 1955 0-85 1056 «148-5 “0117-0095 
202 748-39 0677 11-9 2-1 22-69 0881139 0093 


9-09-94 0-266 11-0 00069 0-330 01131 045 35-4 1166 0189 
OMT 12-7 00068 0-509 O45 1169-0188 -0096° 


213 821-59 PARAFFIN V 0.590 36 0-699 0-665 100-6 
10:40-11:30 V 12-5 -000% 0:90 1-426 0-46 1696 113-4 -0200 - 0108 

216 91-59 4:59- 0313 132 1-396 0169 0107 ofl 
453 268 142 000% 1-421 450 047 355 115-8 0192-014 

218 4:03- 13-7 “00077 72 845 (052 534 162 0109 
219 4:35- 40-664 13-9 -O0077 2-15 053 743) (128-1 on 


13 199 26 43 
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_Ackers’ formula, and have ,Fecomputed much of our this formula. It 


The 19 new in the smooth channel were computed by this 
3 € = 0.00030 ft, which was the value Mr. Ackers found for the 23 runs we had 
~ already reported. _ The average discrepancy of the computed f from the ob- 
7 served f was 5.4%, which corresponds to an error of 2.7% inn or V. Combined 
with the 23 runs with an average error of 3.5% as found by Mr. Ackers, this” 
j makes a weighted average discrepancy of 3.1%. This is slightly less than the - 
» _ average discrepancy |! found by Manning’ s formula or the formula § suggested by 
the authors under the heading “ Discussion of Results,” but the ‘difference is 
- For the tranquil flow in the triangular channel with the rolled-on paraffin 
_ surface, the average value of n was 0.0111 withan average discrepancy of 0. 0003 
- or 2.7%. Using Ackers’ formula an average value of € of 0.00269 ft is found. 
~The f’s computed with this value agree with the observed f’s with an average 
_ discrepancy of 5.5% which is equivalent to a 2.7% discrepancy in n or V. This" 
is just the same as found for Manning’s formula, so it has no advantage. These a fi 
data are for the 8 runs (Nos. 213-220) for which the Froude number did not o> 
exceed 0.53. No runs were made for this roughness with the Froude number > 
between 0.53 and 0.83, but from the results with rectangular roughness it is 
"assumed that tranquil flow would exist up to a Froude number of at least 0.70. 
nee When fitted to Nikuradse’s formula the ¢€ for these 8 runs is 0.0014 ft and 
: the average discrepancy of the computed from the observed value of f is 6.3% 
= corresponds | to a 3.1% error in V. But it is not surprising that this for-_ 
Sai mula does not give good results here, because the flow is in the transition 
rather than the channel range. The expression <~——~ had no value 
“higher than 15.1 instead of 100 as it needs to have for “fully developed turbu-- 
ad The 16 ‘runs (Nos. 175-184 and 188-193) for which the Froude number was 
- more than 1.27 gave an average n of 0.0111, just the same as for the » tranquil | 
- runs. The average discrepancy was 0.0004 or 3.7%. Ackers’ formula gave an 
_ average € of 0.00341 ft and the f’s computed from this value agree withthe ob- 
a ‘served f’s with an average discrepancy of 6. 9% which is equivalent to a 3. 4% 
- discrepancy in 1 inn or V, not significantly better than Manning’s. The ratio of the @ 
_ observed f to the f as computed by Ackers’ formula with ¢€ = 0.00269 ft was also | 
: investigated. The average of these ratios was 1.049 with an average discrep- 
a) ancy of 7.3%which is equivalent to 3.6% in n or V,not significantly better than 
a Manning’s formula. It is surprising however that the f’s averaged 4.9% larger 
oS for ‘rapid than for wh flow, while Manning’s n remained the same. _ a q 


from 0.00685 to 0.0332 ft, and the ratio of f to the f computed ps Ackers’ for- 
mula with e€=0. 00269 ft varied from 1.32 to 2.61. If we | can| believe t the results: 


23 “Hydraulic Practice,” Engineering No. 0.18, U.S, Bur. of Re- 


lam., D » Colo., 1953, 
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of these three runs, the ‘resistance to critical flow. is much greater sie for ‘ 
tranquil or rapid flow. Runs 21 and 22 also fall in this group. Considering the a ~ 7 
waves which form, this high resistance is perhaps not surprising. It was sur- © 
prising however that the run at a Froude number of 0.83 gave more resistance _ 

than the one with a Froude number of 1. 06. (These are the Froude — 


mii 


ors 


17, ASA. AS A FUNCTION OF OF FROUDE NUMBER 


aie However, Mr. Ackers ’ Fig. 1 indicates that even dividing the runs into three nd 
catagories, tranquil, critical, and rapid, is not sufficient. Fig. isa 


— 
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* While | there is a great deal of scatter and the lines drawn are only tentative, 


it is evident that n n is low ai at low Froude numbers, rises to a maximum as crit- ; 
ical flow is approached, and then falls off gradually rather | than suddenly as dl 
Froude number becomes more than one. Similar curves would result from — 
€ against F F as in Mr. Ackers’ Fig. thank him for calling this 
— regard to his last question, the authors planto other of rough- 
ness, but sut probably will make no further tests on on rectangular batten rouganegs. a 


ogi | TABLE 4.—VALUES FOR a, b, AND 


Zegida 
‘Trapezoidal 
Powell 


"Waterways | Ex- 
= periment Sta- 


ton 


Powell & Posey 


= -alog RO andwhen R Racomen ves very | 

when e= 0, — = a log c. Int the literature we the 

_ It seems, therefore, that while the value of a 

= depend upon the shape an the cross section, it depends more on the judg-— 
ment of the experimenter in fitting values | of a, b, and c to his data, and, of 

Sots also on the accuracy of his data. Also, it must be remembered that 
Colebrook’s formula is at least partially empirical and that it still remains an 


‘open ie as to whether it is the best formula to represent the facts. todas + 


‘toon fits our date: as wellas some of the others. This is shown by the ‘ie 
of the in the computed of V 


— 
— | 
2.03 | 2.210 | 0,993 
— 
— 


ne Nov mber 1960 


from the « observed tra quil flow given in Table 5. 


‘ning’s. ‘But forthe 12 in., 6 in.. , and 3 in. spacing of battens it gives the poorest ae 
= 4 3 We have not tested the fit for rapid flow except for the paraffin lining. y 


tr Here the average percentage discrepancy inthe velocity by the three formulas — a 
ig was 3.7 by Manning, 3.6 by Nikuradse, and 3.4 by Ackers. tees agua a 
We will turn now to Mr. van Malde’s discussion. He is very kind to refer 
. to our work as accurate. We of course tried to be as accurate as we could, 
= tt the sources of error which we have already described must have had a siz- 
able effect, as the scatter of points in Fig. 17 indicates. ens age ial 
: eis ‘Varying the slope w while keeping the depth constant is ‘a good w way of varying 


; vs the Froude number. Since we have shown above that the resistance to flow de- P 


wide variation in our slopes for some of the roughnesses. We thank Mr. van 


Malde for calling attention to Kirschmer’s papers. We had seen only the. earlier 

"one in its French form,” and that not until we had torn out our 3- -in. batten = 
i ~ spacing. Since none of Kirschmer’s papers are widely available in the U.S. we 

are reproducing Fig. 14 from his French paper as our Fig. 18. The observa- . 


by Varwick, a student of Kirschmer’ s. ordinate ‘is our f/4an 


— — 


For mula Smooth 


the ahecisen is our R /4. With our apparatus it would have been rom eae 
4 to observe a large number of slopes with the same depth, but we are gladthat 
Varwick was ableto do it. Similar results had already been in Russia 
5 We agree with ‘Mr. van Malde th that these tests give no certain answer as | to 
what formula for channel resistance is “best.” About as far as we can go is. 
to say that for the range covered by these experiments, Manning’s formula — 
_ with n constant for any given absolute roughness and any given Foude number, a 
a discrepancies no larger than - might be expected due to experimental er- — 
rors. Other formulas investigated do not fit the data well enough to warrant © 
- their use. This is especially striking since the authors had begun their study 
F with the belief that a formula definitely better than Manning’s could be found. 7 \ 4 
wide variation of n for. rough channels with changes in the 4 
is perhaps the most important finding of this study. “ee 


ale de I’ Hydraulique, Paris, France, Vol, 15, No, 51 (May-June, 1949), pp. 115-138, 4 a 
24 *Teoriia Podobiia i metodika rascheta gidroteknicheskikh modelei,” (Theory of | 
and Methods of Designof Models for ) by A. P Zegjda, 


For the emooth 
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pends on the Froude number, we agree with we did not make o 
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DISCUSSION 
J We believe one 1 tte van Malde’s remarks in the first paragraph of his es 
discussion might be misinterpreted. _ For tranquil flow in the smooth channel 
4 we found Y = 0.93 which means that f averaged 7% less than by Prandtl’s for- 
mula. This means that our observed V’s averaged 3. greater than Prandtl’s 
¥ formula would give, o or that if our data are correct, , the average error in V by | zi " 
-Prandtl’s formula would be -3.5%. But for the rough channel, using the € which 
gave the best fit, the average discrepancy in V was 2.18% or say 2.2%, with no 
minus sign, as this is the average of the of the discrepancies 
a private communication to the authors, Koloseus, M. ASCE, ha has 
_ made a very thoughtful study of why the two equations given under the heading — 
7 Results from Rough Channel, x ” seem to be about equally satisfactory while the © 
‘first makes f depend only on the relative roughness ¢ /R (but with € a function | 
of the batten spacing A ), and the second makes it depend on A andfp which is 
af function of the Reynolds number. 


ameter which he calls by | the the square root of 


twice the expression (3,2 2). used by Bretting. 25 Koloseus shows that 
f= =a that f could be expressed in terms of R and S, instead tead of — 
4R /e . Then if ina given set of data S« is approximately constant,f seems to 
be a function of R alone, or of € /R alone. _ He points out that with the 12 in. 
batten spacing there was very little variation in S,. The same was true with 
the 24 in. spacing. _ However with the 3 in. and 6 in. spacing there was quite a 
- little variation in Ss. . and he submits Fig. 19 which shows an interesting rela- 
pepo _ He also points out that Nikuradse’ s formula for rough Pipe with € = 3 


= reported for his modified Colebrook formula. It 1 may be added that i 
Koloseus26 has made a very careful study of resistance to flow at high Froude 
_ numbers in a rectangular channel roughened by closely spaced cubes. _ oak’ 
a: The authors wondered why Buckley called his formula Manning’s formula. 

- This matter has now been cleared up by one of the authors. epi 
a. authors are sorry that there was no discussion of the ‘capillary rise” i} 
4 _mentioned under the heading ‘ Results from Rough Channel. ” It seems that this 

er known phenomenon may ‘account in whole or in . part for the apparent in- 

crease in roughness in critical and rapid flow. 
__ During the summer of 1959 a very large number of pitot ‘measurements of 
the velocity were made for v uniform flow at a relatively large depth. In order 

a to do this the flow had to be held steady for many hours at a time. _ ‘This took | 
less man-power than establishing new steady flow rates at frequent intervals. 

Advantage was taken of the long duration of the uniform flow condition and of | 
_ the availability of ‘personnel to determine water-surface and bottom profiles — 


more thoroughly than before. —* were read every 2 ft along the length of © 


— “A Set of Practical Formulas Based on n Recent Experimental Research,” by yA. < 
_ Bretting, Report, 2nd Meeting of the Internatl. Assn. for Hydr. Research, 1948, pp, 399- 


86 Dissertation, by H. J. Koloseus, Presented to the State University of Iowa, nt 
City, Iowa, in partial fulfilment of the requirements for the degree of Doctor of — 
-ophy, August, 1958, Presented at the October, 1958, ASCE Convention, = 
27 of s Formula,” ” by Ralph ‘Powell, Journal of Geophysical R 
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rom the observed value for flow given Table 5. the 


_ Here the a average percentage discrepancy inthe velocity by the three a 
3.7 by Manning, 3.6 by Nikuradse, and3.4by Ackers. 
Wes will turn now to Mr. van Malde’s discussion. He is very kind to refer is 
to our work as. accurate. We of ‘course tried to be as accurate as we could, 

_ ‘but the sources of error which we have already described must have had a siz- 
able effect, as the scatter of points in Fig. 17 indicates. = = © Bare 
ee the slope while keeping the depth constant is a good way of hee 7 
a ‘the Froude number. Since we have shown above that the resistance to flow de- 

pends on the Froude number, we agree with Mr. van Malde that is influenced 
it by S. Unfortunately, we had not foreseen this effect so that we did not make a 
_ wide variation in our slopes for some of the e roughnesses. We thank M1 Mr. van 
Malde for calling attention to Kirschmer’s papers. | We had seen only the earlier 
i one in its French form, 9 and that not until we had torn out our 3- in. batten — 
_ spacing. Since none of Kirschmer’s papers are widely available in the U.S. we 
are reproducing Fig. 14 from his French paper as our Fig. 18. The observa- ot, 

tions are by Varwick, a student of 8. ordinate lis our and 


a - fit. * We have 1 not tested ‘the fit for ‘rapid flow “except for the paraffin a. 


in 


the abscissa is our R /4, ‘With our apparatus it would been difficult 


_ to observe a large number of slopes with the same depth, but we are glad that 
8 Varwick was able to ¢ do it. re ‘had already been obtained in Russia 
agree with Mr. van Malde that these tests give no certain as 
"what formula for channel resistance is “best.” About as far as we can go is 
to Say that for the range covered by these > experiments, Manning’s: formula 
With n constant fora any given absolute roughness and any given Foude “number, as 
P- - gives discrepancies no larger than might be expected due to experimental er- ei 
‘rors. Other formulas investigated do not fit the data well enough to warrant | 
thelr use. This is especially striking since the authors had begun their ae 
_ with the belief that a formula definitely better than Manning’s could be found. 
wide variation of n for rough channels with the number 
“Pertes de dans les conduits forcees et canaux decouverts,” Revue Gen- 


a de l’Hydraulique, Paris, France, Vol, 15, No. 51 (May-June, 1949), pp. 115-138. Phd 
24 *Teoriia Podobiia i metodika rascheta gidroteknicheskikh modelei,” (Theory of 
Similarity and Methods of Designof Models for Hydraulic A. P. 
Union of Soviet Socialist Republics, 1938, 
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= 0.93 which means that averaged 7% less than by Prandtl’s for-_ 
4 mula. - This means that our observed V’s averaged 3.5% greater than Prandtl’s a] 


-Prandtl’s formula would be -3, 5% ‘But for the rough channel, using the which 
gave the bast fit, the average discrepancy in V was 2.18% or say 2.2%, with no 
- minus sign, as this is the average of the absolute values of the discrepancies ; 
_ In a private communication to the authors, H. J. oe ASCE, has 
made a very thoughtful study of why the two equations given under the heading ~ 
“Results from Rough Channel,” seem to be about equally satisfactory while the © 
7, first makesf depend only on on the relative roughness €/R (but with ea function 
4 of the batten spacing ), and the second makes it depend on A and fp which is > 


a function « of the number. introduces another | dimensionless par- 


ameter which he calls = Se (which by the w: is the of 


twice the expression X = gs s used by Bretting. 25 k shows that 
so that f could be cada in terms of R and S, instead of 
Then if ina given set of data S, is constant,f seems to 
‘be a function of R alone, or of €/Ralone. He points out that with the 12 in. | 
batten spacing there was very little variation in S,. The s same was true with | 
the 24 in. spacing. _Bowever with the 3 in. and 6 in. spacing there was quite a 
little variation in Ss. and he submits Fig. 19 which shows an interesting rela- 
; tionship. He also points out that Nikuradse’s formula for rough pipe 1 withe = 
0.00040 ft fits the smooth channel tranquil data with an average discrepancy in 
2 f of 6.6% which is somewhat better than the 3.5% for Vv (7.1% for f) which Mr. _ 
* Ackers reported for his modified Colebrook formula. It may be added that: Mr. 


Koloseus26 has made a very careful study of resistance to flow at high ian ry a i 
numbers in a rectangular channel roughened by closely spaced cubes. a 
"The authors wondered why Buckley called his formula Manning’ > sentry 
This matter has now been cleared up by one of the authors. 27 
-_ The authors are sorry that there was no discussion of the “capillary a - 
mentioned under the heading * Results from Rough Channel.” It seems that this a 
+ little known phenomenon may account in whole or in part for the ‘apparent in- 
in roughness in critical and rapid flow. _ 
- During the summer of 1959 a very large number of -pitot n measurements of | 
, the velocity were made for uniform flow at a relatively large depth. - In order © 
to do this the flow had to be held steady for many hours at a time. . This took x 
less ma 2n-power than establishing new steady flow rates at ioeanbae intervals. 
Advantage was taken of the long duration of the uniform flow condition and of 
the availability of personnel to determine water-surface and bottom profiles 4 
more thoroughly than before. Values were read every 2 ft along t the length = 


25 “A Set of Practical Formulas Based on Recent Experimental Research,” by A. 


‘Bretting, Report, 2nd Meeting of the Internatl. Assn, for r Hydr. Research, 1948, pp. 399- 


26 Dissertation, by H. J. Koloseus, Presented to the State University of Iowa, Iowa 


City, Iowa, in partial fulfilment of the requirements for the degree of Doctor of Philos- 
ony August, 1958, Presented at the October, 1958, ASCE Convention, _ 
of s Formula,” Ww. Powell, Journal of Geophysical Re- 
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DISCUSSION 
the flume with the exception of only a few points that sceiahindead not be seen from 
— results are shown in Fig. 20 compared with the values. ee: the mean A 
water-| surface and bottom profiles as ordinarily determined, that is, by the 4 
— method with bottom elevations determined every 8 ft and water- -surface | 
E ¥ elevations every 24 ft. It can be seen that for this flow condition the ordinary 
_ method determined the slopes fairly well, but underestimated the depths by © 4 
_ perhaps 0.06 ft or about 2. 85%. The Froude number for this run was 0. 88 which 
meant that surface waves hada relatively large. effect. ‘Together with the large 
- deflection of the bottom due tothe heavy water load, the underestimation of the - 


depth in this test nearly an an extreme forthe pres- 
‘Sa _ Attempts were also made to get transverse water-surface profiles, but set- 
- ting gages to the fluctuating water-surface proved too difficult a task. With _ 
: the slight irregularities of the flume obviously c causing much larger water- er 
face waves, they would be of doubtful value. 
__‘The fact that the channel showed less resistance to flow than the e correspond~ 
: 8 ing smooth pipe, has troubled both the authors and the discussers. However, 
; A ay and Robertson28 reported that tests made at the University of Illinois, 


= = W. M. Owen, showed f for a smooth triangular flume less + for the cor- ; . 
responding 


c el channel by Varwick and on a rectangular one by Raju. ‘These indicate 

‘uted a, about 1.04 for the trapezoidal channel and 1.12 for the rectangular one. Kir- ; 
he schmer feels that this is to be expected because the effect of secondary flow 
would be greater in the rectangular section than in the trapezoidal. It would 

. Pe seem that a a triangular section should give y at least as great as a trapezoidal Fa 

our channel bottom ‘may have averaged as mech as s 0. 01 ft 


below the values measured at the 8-ft panel points, so that the tabulated values 


The elevation of 
am”, % Of yn may have been as much as 0.01 ft too small. The 42 runs for the smooth 


flume with tranquil flow were 1 recomputed with each Yn increased 0.010 ft — “8 
the value given in Tables 1 and 3. The resulting y was 1.007 instead of 0.936. 
Bhi It was also found that the average discrepancy from the Prandtl- Nikuradse — 4 
formula for smooth pipes _ given by these assumed values was only 5.4% for . ¥ q 
which is equivalent to 2.7% for V or n, compared with 3.3% using bottom ele- 
i vations as recorded. But the close agreement with Owen’s s res results en — 


as to whether this this should be 


28 Discussion of Flow at Small Numbers,” by Wallace 


M. Robertson, Transactions, ASCE, Vol, 123, 1958, Pp. 707, 
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‘Therefore Owen’s would be 0.295/0.316 = 0.934 which is almost exactly the 
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THE SETTLING PROPERTIES OF aeaenar 


by N. Claes H. Fischerstroem 
CLAES H. FISCHERSTROEM, TROEM,! M. . M. ASCE.—Mr. McLaughlin’ S paper isa 
- valuable contribution to our knowledge of the sedimentation process, ,especially — 
; 4 since it includes results from determinations of the settling of suspensions, 3 a 
well as practical hints regarding the execution of the tests. ts” 
= _ Although, according to the introduction, the paper is limited to the settling 2 
‘ properties of suspensions and tothe we of the results in computing removals, 
the author, in section 4, also discusses ‘ practical engineering problems, ” He 
__- makes some statements regarding the influence of ‘other factors” onthe design | 
4 a settling tank. The writer has some questions to ask regarding the author's 


opinions, and presents some further conclusions from the results. 


amount of suspended solids which passes a certain plane, calculated by = 
mean concentration and mean velocity of the solids under non-disturbing flow 
- conditions. The equation arrived at covers the settling in both the free and the 
. 4 hindered zone. The main value of the equation, according to the author, is to 
= _ show what research and development are necessary for complete calculations. 
But are there really any new items of research mentioned in the paper \ which 
are not already known However, the analysis isa good basis for +!:- oxten- - 
a: The flow of water ordinarily takes place mainly inthe free settling zone and a 
“the equation is not very useful bo the purpose of drawing conclusions | regarding © i 
this part of the settling basin. Such an equation is still lacking. == — 
the end of section 2 the author discusses scale models.  Sertling studies 
va * models cannot ignore the flow conditions. The writer believes that a 
actual conditions ‘it is more important to study the flow conditions in the 
scale model and the settling in the tubes. If the flow is satisfactory, and only 
une then (not necessarily ideal) will the sedimentation in the prototype approxi- ¥ 
Section 3 gives valuable information regarding results and discussion of 
- tests by the pipette and the multiple- -depth analysis. The latter will be =e 
In Fig. 5 the settling in a tube is illustrated by a concentration dagrem. rt 
the” dashed lines represent concentration at constant z/t, shouldn’t co orl 
start at @ (z,t) = 100%? Figs. 7 and 9 represents the settling in a tube of ‘ee 
suspensions. Shouldn’t the full lines (time curves) start at ¢ (z,t) = 0 


the dashed lines (zA curves) start at (z, t) = = = 100%? As is seen from 


4 


“limit values are taken into consideration. 5 This first part of the curves can x 
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never be e vertical but will always be or more or less 
which is important when drawing conclusions 
pe: From ‘Figs. 7 and 9 we can see that the settling “properties of the ai “a 
: sions are bad. If there has been no flocculation period there will be a chemical 
reaction and we are not studying a settling process alone. It is known that 
- Fitch’s sedimentation test is carried out after dosing of the (whiting and) oo os 
-risul by a rapid mix, but without flocculation « of any kind. Is the _ author’s test = 
* with betonite clay and alum also carried out without preceding flocculation ? 
Or if it is flocculated, what are the details (stirring, time etc) ? Is the ma- _ 
:- terial tested really a suspension ? Isn’ t it a suspension under formation by a ‘al ay 
a), process which is not only time-consuming but also influenced by turbulence, 4 
concentration, pH a.s. -o. The writer does not mean to say that such suspensions a 
would not occur, but nevertheless these tests might be too to allow 
section 4, part b, the author states ‘This approximation willoften bebet- __ 
a ter than making an elaborate study of the hydraulics of the tank while ignoring 
= the properties of suspension. ” Why ignore any of the important factors? The | 
"underestimation of the hydraulics seems tobe not uncommon for the laboratory 
- fellow. In contrast to the field engineer, he gets no experience from the distur- Bi 
‘bances from density currents, side currents, wind, turbulence etc. The settling 
properties are easily determined, but the hydraulic properties are not 
a perhaps why so many existing basins, and also so many basins under con- 
struction, are unsatisfactory in hydraulic respects. 
‘In the same section, part d, the author states, “It will now be demonstrated 
whether ‘overflow rate or detention time should be used depends upon ‘the: 
nature of the suspension. Indeed, for a single tank the removal ratio may wa 
closely related to overflow rate for one suspension, to detention time for a 
second, and to neither fora third.” This sounds simple, but Mr. Fitch was 
a < perhaps nearer the truth when, in a discussion “regarding influence | on settling 
9 of Class II suspensions by the overflow rate, namely, the detention time, he 
a wrote, both are working.” Will it ever be possible that a sedimentation oy 
depends upon time ? what does the author mean by “ neither for a third 
.. illustrate the preceding, t the author starts with a theoretical discussion - 
Pa of the concentration diagram. Overflow rate, he says, is determining if z/t _ 
¥ curves are straight and sential tothe z (depth) axis and detention time is de- 
q ons termining if the concentration (time) | curves are parallel to the same axis. If 
we deal with suspensions which settle, the latter will never be possible because = 
_ all time curves but the one at t = 0 must start at the concentration ¢ (z,t) = 0. 
Inthe thin surface layer, sedimentation results immediately in a complete 
} ‘moval. In the surface layer there is no flocculation. Also in the next, shallow 
_ layers, clarification is more rapid than in the deeper layers. Take, for in- Be 
a stance, the author’s curve for T = 700 sec (Fig. 7). Ata depth of 10cm we i 
have a removal of nearly 50%, at 20 cm of about 40%, and at a depth of 100 cm 
23% (at the same time or volume load). At the ‘same 
Ee will always give a higher removal than a deeper basin. It i 
_ consequently, even with a suspension with poor settling properties, more ef- 


a ficacious to u use, ,for example, a 0.2 in. basin with 5 times as large surface than 4 
. use a 1.0 m basin. The shallow basin furthermore will have lower turbulence = 
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o illustrate the relation iaiiaiia ait rate and time for certain kindof 

ne we can study the writer’s Fig. 1, based upon the results published 

_ by Fitch. This diagram, which is of the same type as the author’s though re- 

arranged for | technical purposes, shows the residual after sedimentation: 


ta a. at constant settling time (volume — — variable depth (overflow rate, 


| 


ESIDUA 


PERCENT R 


DEPTH OF SETTLING CHANNEL, ft 


aca The time curves show us that at constant time, 1.5 hr, for example, the re- 
: sidual increases with the basin depth, or with a deeper basin with a smaller a 
H area. From these curves we may draw the conclusion that the best way - ~- -- from a 4 


a sedimentation point of view - use a volume is to give the basin’ 
large an area as 


— 
— 
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The rate curves indicate that the residual decreases with increas- 
_ ing depth, which ‘means with increasing time or a decreased volume load. The a 
_ increased efficiency must be bought bya greater basin volume. The effect o 
. served, which depends upon an accelerated settling because of flocculation, 
f shows us, however, that there is a certain relation between the overflow rate 
and settling time for every required settling efficiency, = 
_ We can, for example, take a recovery process with 90% settling efficiency bys 
residual) and a purification process with 2% residual. _Corre- 
sponding overflow rates and detention times will 


Residual < rate per detention time 


1/2.5 hr 4/2 hr 0. 0.2/1, Shr 
From the we can confirm the often expressed relation: “the lower 
_ the overflow rate the shorter the time,” in order | to reach the same removal. _ 7 
_ ‘The only possibility to increase the volume load and to reduce the basin — 
- is in fact by decreasing the overflow rate, environmental factors being un- 
_ changed. Very low residuals, near 0%, which we now require in certain pro-— 


cesses, seem only to be attained at ; very low overflow rates and wer small 


= per basin unit volume and other factors, such as available space, stability 


We can also consider tha overflow rate and time talleetinde in this way. a 


_ free particles of a suspension at the time t = 0, for example, before any 
flocculation can take place, have certain n settling velocities, which we represent 
in a settling diagram, we can always determine the overflow rate for a required 
removal at depth = 0 or time = 0 for this suspension. _ Because we have a cer- 
tain depth or or time, , the actual removal will be increased, and to keep the re- - 
‘quired removal we can increase the volume load. Because the incluence of 2 
flocculation will increase the settling properties only up to a limited depth aa 
Es above which the particle speed will disperse the floc again - this improvement — 
-~ increased time or depth will take place only to a certain limit. Increasing , 
}. the time above this point will not increase the removal ratio, and any attempt S 
; a to increase the overflow rate load will be a failure. From the diagram we can | 
Bs see that the constant time curves will always lead to zero at decreased over- 
_ flow rate; curves for constant | overflow rate would finally be parallel tothe = 
ye To summarize, it will always be possible to get a required removal, ,even 
a by reducing the time ~0, by lowering the overflow rate, but this will not be — 
a possible by too high an overflow rate, | even if the time is “increased to — 0. 


; Although both overflow rate and time can be said to be limiting factors, it fol- _ 

™ lows that the overflow rate might be sa said tobe of primary importance when an 

.. seems most probable that the very shallow basin will have a future, pro-— 


| 
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& 
le story basins or multiple story basin 
be achieved by shallow sing 
t 
— 
vided it 18 designed on technica ry. Todetermine the settling of a 
ime settling is unsatisfactory. iod will require new 


DISCUSSION 
methods. 


lence sais enough in the “tube test.” The writer believes that we perhaps — 
must change to a continous flow method. 
Fhe investigations made by the author, Fitch, and others during recent years, 
“n. given us an idea of the relation between the overflow rate and time for __ 
- some special types of suspensions. Further investigations regarding orl 
te types are desirable. What seems to be especially valuable is what has pre-— =. 
F - viously only been assumed or felt intuitively or possibly doubted or denied, that i 
for flocculated suspensions of Class II there is no direct proportionality 


ik drawn from the overflow rate theory, are  oeekan ti best basin for all ae 
suspensions is oul the one with depth and time 0 and surface o. The in- 


- that the only constructional way | to increase the. ‘volume load is to decrease the 
7 _ depth, not to increase it. The conclusion that a settling process can be depend- 
_ ent only upon the time, seems to be drawn from concentration curves lacking 
the important part just below the water surface. | 
The prior discussion is not meant as a criticism, a, only to give a more » com-— 
plete view on the conclusions drawn from the author’s settling theories and 
“S - experimental results. We must be grateful to the author for his excellent work — 
the writer agrees with his" conclusions, with the exception possibly for 
es _ The writer believes, namely, that when a choice must be made, it may . 
nek more er. to study the flow in va" basin than the properties of the sus- 


rate, the speed depending upon the type of suspension. The investigations also 
q 
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IVAN J. “MOSKVITINOFF, ASCE.— This pa paper is a valuable contribution 
to the history of engineering science in the United States in providing the col- 

_ lection and systematization of the information scattered in many publications | 7 
and magazines. | The author paid a tribute to the American pioneers in hydrom-_ 
Be and to the engineers and professors who built the scientific and technical | 
base of modern hydrometry. The paper contains important information on a 
parallel development of of the e hydrometric c methods ar and devices in some soning oe 

— -‘Deseribing the investigation on velocity pulsation at Burlington, the author | 
concludes: * The influence of pulsation cannot be eliminated by short observa- 
tions at a great number of points as it is generally believed.’ * The resu'ts of 
writer’s measurements of the velocity distribution in many irrigation 
nels of different sizes in the Syr- Daria in Turkeston, Russia, in 


a short ‘period of “observation, were sufficiently “accurate in comparison with 
the rates of flow obtained by two- point method with a long. period “ cnet 
tag The writer vigorously supports the author’s proposal to use a more alent. e 
_ The writer hopes also that the author will be able to fulfil his intention to — 
publish the history of the modern hydrometry after the First World War, for 
_ this would be a valuable contribution to engineering 


wt 
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by F. P. Linaweaver, 3 John Geyer, and Jerome B 


A. M. ASCE, JOHN C. GEYER,““F. ASCE, and ; 
"JEROME WOLFF? M. "ASCE. —There are many problems that continually _ 
. arise in the operation and planning of a water distribution system that require 
a. a hydraulic | analysis. If the basic analyses of average daily, maximum n daily, 
ie 
peak hourly conditions are available for each pressure zone of a water dis- 
tribution system, some of these problems can be answered by the useof propor- a 


2s loads as presented by Mr. ‘McPherson. For systems without storage or og 


with storage at or near t the pumping station, proportional load hydraulics ac- 
curately produces the desired answers, _ The author has used test data to show ; 
that fora single pumping station, and a reservoir elsewhere on the system, the “ag 
: head loss between the pumps and any given point in the system is afunctionof _ 
the ratio of pumping rate to demand rate. The discovery that such 
s exist should prove very useful, particularly when nanalyses are made by long } —_ 4 
methods or by use of digital computer, 
The method appears to have usefulness for solving problems other than those 4 
"involving pumping versus storage relationships. One recurring problem : is to. 
_ determine the effect of proposed housing developments onthe distribution sys- 
rc tem. This would normally be a concentration of loading at a point and, there- 
ay fore, walle be a non-proportional load similar to a fire flow. The flows dl 
a or out of equalizing storage facilities are in effect non- -proportional loads, and, 
_ therefore, a procedure might be determined whereby Eq. 2 can be used for - 
e - In making analyses involving a storage reservoir, the fact of discontinuity 
in the situation when pumps are changed or when the reservoir either fills or 
-_ pet vote must be dealt with. This requires, in particular, that the time of fill- 
_ ing andof emptying the storage reservoir must be determined, because at these 
<4 times there occur drastic changes in the hydraulic situation. . Following of the fs 
\j cyclical variations in storage levels would appear to require the reverse ap- ¥ 


In the caseof a more as Fig. 7, the determination 

-_ the appropriate Qp/Qa ratio to account for various pump curves of a par- 

ticular pumping station, and the changes inlevels of tanks andstandpipes would 

me - appear to be very difficult. The author states that “The computations for the 
Recessive hourly balancing of pump-network-storage functions can now be 

much more simply using Eq. 2,. ” In Eq. 2, Qq, n and would be 


in order to determine Qp must be determined. If there is more 
4 


January, 1960, by M.B. McPherson. “it 
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Prof, and Chmn, of Sanitary Engrg., Johns Hopkins Univ., Baltimore, 
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i of the quantities of water being supplied by each source would not be known. . 


would appear that in order to determine the ratios of Qp/Qa as shown in col- es 


—umns (3) and (10) of Table 4, an analysis of successive balancing would be re- - 
but accounting for 


still be required. It would be interesting to know ‘the author’ S ‘procedure _ 
_ determining columns (3) and (10) of Table 4 especially in the Number 4 series — 


In analyzing complex systems on the Baltimore Analyzers, a procedure for 


balancing is used, whereby the operating levels of the storage tanks or reser- x 


voirs are determined on the initial test, and then the system is operated through — 
r the day just as th the pump operator or an | automatic station would operate. ie 
stead of using an hour-to-hour balancing, the period between tests is ion 
needed to obtain results within the desired accuracy. 
:~ a The operating level is defined as the high water level haben the early 2 a.m. — 


reading the difference in head between the sources and the storage facilities. 


_-- The question, “Is it not feasible to consider writing a design program for a 
digital computer which would permit direct determination of the best combina- 4 
_— of pipes to satisfy prescribed conditions defined by means of Eq. 2?, is a 
very challenging one. It could be expanded to adetermination of the best com- 
i bination of pumps, pipes, and storage facilities by modifying Mr. ‘McPherson’s Ss 


a 


__5 “Distribution Design Considerations,” by M.B. McPherson and J, V. Radziul, Jour- 

AWWA, Vol. 51, April, 1959, pp. 489-502, 
“Peak Demand George G. ‘Schmid, 
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SCOUR AT BRIDGE 


cussion by D. Joglekar 
. Rylands Thomas Mushtaq Ahmad, : and Pier Luigi Romita_ 3 
is agreed that each river and each reach must be 
Studied to understand it: its individual, , almost personalized, characteristics and 
that scour at bridge piers is closely related to (a) the river conditions 


stream and downstream (such as indicated by flood hydrographs from year to 
a year);  (b) the meandering tendency of the river that, in turn, depends on the 


“ing in its alluvial plain or has its sides and bed resistant to scour, 
= river is flashy or has sustained floods, and (e) the constriction of the riv- > 
er section caused by the bridge. The reach of the river also has an important a : 
: ‘a on scour at bridge piers, for example, bridges on hill torrents, bridges — 
on alluvial plains with mild slopes, and bridges on the tidal reaches of the ae 
river. In a braided river near a hilly region the river changes are so rapid : i 
that it is very difficult to estimate the waterway required for bridges across 
various streams of the braided river. An illuminating case of violent move- 
- ment of Manas River!7 in Assam can be cited. bie 1909, three bridges were — 
provided on the three tributaries of this river: 
4 on the Manas River which had 10 spans of - 100 Ree: 
4 2 on the Bholookadoba Branch - 2 spans of ' 75 ft, 


on the Beki Branch - 4 spans of 12 ft. 


Due to changes in river courses, the b bridge had to be rebuilt several times 0 on 
the Bholookadoba Branch between the e period 1909 to 1945: 1945: 


of 250 ft and 2 of 150 ft, 


4xt 12. ft, 


OER 20 ft, 


river is is still uncontrollable ‘in spite of various ‘training m measures. 


February, 1960, by E,M. Laursen, 

a 16 Adviser, Central Board of Irrigation and Power, Poona, India, st f a 
“River Railway Board Technical Paper No, 334, 
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a. Station, Poona, India, in 1938 and 1939 for finding scour round a single | pier 
placed in the center of a parallel sided flume, embedded in sand of mean 
eter 0.29 mm. The following relation was worked out: 
0.78 


pier below water level. As it was difficult to correlate this with the depth of 
at prototype piers, due to q (intensity of discharge per foot) depending 
upon the curvature of the river upstream - #3 which varies from river to river, | 
it was considered desirable to study cases of actual scour in prototypes and : 
work out a general, empirical relationship. Besides, it has to be remembered > 
that the angle of repose of bed material in the model and the prototype is the 
neces hence, the extent of scour in plan in the vertically distorted model is 
feed always relatively greater than in the prototype. This in effect reduces 
fe discharge intensity at the pier due to greater dispersion of flow and hence — 
the depths of Scour obtained the model would be relatively less. Datal9 


were, therefore, collected for scour round piers of various bridges 
structed in India and relationship! was worked out as follows: 


in which Qis flood cu sec; dg is the maximum 
of scour below highest flood level; f is = 1.76 Vm, and m is the mean grade of 
P bed materialin millimeters. In Eq. 15, a representative f value has to be used. 7 
A From bore data, values of f for each strata is to be worked out to ascertain 


that the anticipated co is not based on the f value which is higher than that 


- depth of scour to provide adequate grip length. Where this cannot be done, one, high 
level stone protection, though costly in the long run, has to be employed. An- _ 
other advantage of deep foundation is that, because of increased side friction 
onthe pier sides embedded in sand, the load bearing capacity of the pier in-— 

a creases considerably as compared to that of a pier with shallow foundation. 
Generally this additional load bearing capacity is not taken into yen sr all 


in the design but is kept as a margin of safety. _ The previous railway practice — 
was to work out the depth of bridge foundations according to the = 


Spring? 22 and Gales.** In the case of shallow 


Pe i9 “Behaviour and control of rivers and canals with the aid of models - part II, » Re- 

gearch Publication No, 13 » Chapter VIII, Central Water and Power Research Sta., Poona, 
20 Technical Annual Report, Central Water and Power Research Station, Poona, — 
: “Stable Channels in Alluvium,” by G, Lacey, Journal of the Institution of Civil a 

“River Training and Control of the Gute J. Spring, ‘Rail- 

way Board’s Technical Publication No, 153. 


a 

Technical Annual Report, Central Water and Power Research Station, Poona, 


stitution of Civil of Civil Engineers, December, 1938, Paper No. 5167, 
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DISCUSSION 
has to be pr provided — stone. “pitching and nd if this is at high level, a lot | of the 
- stone (stone used generally weighs 80 to 120 lb and is called one man stone) ~ 
is washed away downstream due to turbulence and has to be replenished, even 
meen to ensure the safety of the bridge. _ Due to this turbulence very 
a: _ deep scour occurs downstream of the piers as in the case of Hardinge Bridge24 a 
the river Ganges, depth of scour being of the order of 200 ft. 
oe _ The current railway practice is to provide a grip length for the pier equal — 
to half the depth of scoured bed below H.H.F.L. so that the total length of a 
seal Experiments25, 26 were carried out at the Central Water and Power Research 
Station, ‘Poona, India, for | testing the design of training works, waterway, and 
>= of piers of arailway bridge at Mokameh on the river Ganges near Patna 
(Bihar State). The pier foundation of this bridge has been taken to a depth of 
about 200 ft below H.H.F.L. which is equal to 3D(Lacey). As the 
are deep enough, protection by way of stone pitching round piers is not pro- 


4 Various cases of bridges, for which rivers had to be trained to ensure the 
safety of piers, have been experimented on. These experiments are des described 
in the Technical Annual Reports of the Central Water and Power Research Sta- j 

_ tion, Poona, India, for the year 1937-38, 1938-39, 1939- 40, 1940- 41, 1944 to 41, ; 


Inthe case of flash flood type rivers, fixing the waterway is snais more dif- 
ficult. In such cases the flood rises and falls so rapidly that the river has no 
_ time to scour its bed with the result that the afflux (difference in water sur- 
“In face upstream and downstream of bridge) increases enormously; the bridge is 
é likely to fail by outflanking. A railway bridge on Luni River in Rajasthan State2? 


a. In the case of inerodible bed material, it is difficult to estimate the maxi- a 
mum depth of scour. Hydraulic model experiments are unable to reproduce ~ 
_ this scour due to obvious limitations. Recourse has, therefore, to be taken to 
field data. In some cases, the maximum flood level is underestimated and the 
- waterway provided is insufficient. If the bed is inerodible afflux increases be-_ 
yond the safe limit, with the result that standing wave conditions prevail down- ¥ 
} stream of the bridge, which lead to undermining of bridge piers. The railway 
; ae bridge on Yeshwantpur River28 in Andhra Pradesh failed and collapsed for 


= The writer agrees with the asia that afflux caused by a bridge depends so _— 


wwe 


much on the erodibility of the bed material. Thus in the case of the railway — 
bridge on the River Ganges25 at Mokameh, the afflux caused by the bridge was 
iv only a couple of inches for a flood of 18,000,000 cu sec. This was because the _ 
flood + was sustained and the river scoured its bed as the discharge increased. es P 
_ 24 Investigations carried out by means of models at the Khadakwasla Hydrodynamics — 
- Research Station, near Poona in connection with the protection of the Hardinge Bridge 
which spans the river Ganges near Paksey, by C.C. Inglis and D. V. Joglekar, East Ben- 
gal Railway, Public Works Department, Bombay, India, 1936, Technical } Paper No, 55, 
em _ 25 Technical Annual Report, Central Water and Power Research Sta., Poona, India, ie 


2 6 “Manual on River behaviour, control, and training,” Ch. VI, Pub. No, 60, Central 
Technical Annual Report, Water and Power Research Station, ia, 
Technical Annual Report, Central Water and Power Research Matton, Poona, India, 
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In the « case of rivers with flash floods or with inerodible bed, ‘full afflux has to 
be allowed for in the design according to standard formulas. 
As emphasized by the authcr, more field experiments are necessary to im- i 
prove our method of estimating scour at bridge piers. 
OW. J. BAUER,‘ 29 M. ASCE.—The writer was inattendance at the University 
bs of Iowa, when the research described by Mr. Laursen was getting under way, 
- and is, therefore, aware of its pioneering nature. The purpose of this discus- a r 
‘sion is to consider the application of the results presented by the author to 7 : 


= 


= The fundamental purpose of a bridge over a stream is to separate the 
hu uman traffic and the flowing water. It is commonly accepted that provision © 
wai for the passage of the maximum conceivable discharge through the waterway 
beneath the is economically undesirable, It is good practice to pro- 
_ vide for the safe, although infrequent, overtopping of the roadway by a rare 
flood. If the design is adequate, after the passage of such a flood, the high- © 
-* is immediately ready for service with only minor repairs being required 
during subsequent routine maintenance. 
7 In order to achieve this objective, the backwater produced by the struc- 
ture must be small at the stage of incipient overtopping. The writer has used 
values of 0.5 ft to 1.0 ft as being reasonably small for the backwater at this : 
stage. ~The design problem then becomes how to provide for the passage of _ 
| corresponding to the stage of incipient overtopping without exc exceed- : 
Some of the flow area required will exist beneath the elevation of the hatte 
7 flow bed of the stream in accordance with the reasoning set forth by Mr. 
4 Laursen under the heading “ Local Scour at Piers” and Abutments: — 
x at Bridge Crossings. ” The extent to which such area is ; available to the flow 
during flood may be controlled by suitable scour protection, 
_ At bridge sites in West Virginia, at which scour protection of broken rock 
readily | available, the writer has used a design velocity beneath the bridge 
of between two and three times the - typical velocity in the stream during 
_ flood. This was accomplished without ‘exceeding the allowable backwater at b. . 
the stage of incipient overtopping. Such a relationship, between the velocity _ 
¥ in the contracted section and that in the approach flow, gives sufficient as- JA 
= jj surance that the material deposited over the broken rock fill during low flow _ 
&§ ci will be scoured out during flood, provided the depth of scour required is not : 
A! excessive. Fig. 6 indicates the relationship between depth of scour at an cor 
croaching | abutment and a parameter dependent on the degree of contraction 


a of the natural waterway by the structure. Fig. 6 indicates the attaining of _ 
‘ _ depths of scour equal to the depth of flow at relatively minor amounts of = ; 
_ contraction. It is, therefore, not difficult to imagine circumstances in which | 
a it would be possible to provide | as much waterway area beneath the low- flow 
In the instances that the writer has in mind, the waterway area provided 
Ni beneath the low-flow bed of the stream was about 20% « of the total waterway _ 
area : area at the stage of incipient overtopping. The depth of the waterway — 
¥ beneath the low-flow bed was only — 40% of the nominal depth of flow at | 
the stage of incipient overtopping. Nevertheless, significant cost savings 
were achieved compared to a longer bridge ' with equal waterway area entirely i 
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which the length of the bridge might be out in half, the AE, 

maining the same. This has obvious economic significance. 

oe <a The natural filling of the scour hole during the recession of the. flood ie 
to the. depth of the pool water that left under 


Sediment, 
Flow stream 

_ bed 


= 


rodible sediment, 
flow bec of 


(6) Section : 


Fic. —PROPOSED 


flow. As the beds of are commonly 
up of alternating pools and bars | anyway, the addition of another small pool | 
beneath the is of little consequence. 12 type of 
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J. TISON. 


of scour at piers nin eaten for the case in which sediment is Supplied i. 
scour hole, inextremely interesting, 


The author presents an analysis in which he compares the flow without 


ontraction and the flow due toa long contraction. He then uses, for both 
“flows, an approxima e form of the total sediment load relationship recently _ 
a4 s proposed by himself and consideres a bridge crossing as a long a by 
4" ———— in such an extreme that it has only a beginning and an end. & 
_ The writer has found that a scour at piers and abutments also took place | be 


when the contraction of the flow was without significance and introduced the 

_ idea that the scour at piers may be affected by the curvature of the stream-_ i 

a pier with an arbitrary shape is placed in a it produces 
(Fig. 13) a first curvature with a center 0; and a radius pj,followed bya 

) curvature 02 Po. In the neighborhood of the bank of the river, the pier “igs 7 

ae oe no action on the ‘direction of the streamlines, cd. But, considering th the it 
region of the first curvature in the neighborhood of the surface, a line such > 


, _AB, tangent in each point to the principal normals of the streamlines, the pak 


B' are on a same line perpendicular on the bottom, 


The variation of the B follows the hydrosatic 
that: 


vis the velocity in the neighborhood of the bottom). 


~ 
— 


A  eeneeanence of. Eq. 19 is that the motion | cannot be parallel tot 
‘This diving motion will have ‘another consequence : a local attack of the 
bottom under the influence of the first curvature 01 1. Evidently, the im- 
. — of scour will depend on the value of the vertical component of of the b 
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ag 


| 
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Gi oe, and this vertical component will depend on the value of =e 
‘example, small values’ of p will the e depth of scour. 

ments with different models of piers, with the same length and breadth 


with» 
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(reduced values of P) produced a scour of 113 mm with 


a ‘. tie of 301 per sec in a flume with a breadth of 70 cm and a height  - 
10.5 cm, The length of the pier was 24 cm and its breadth was 6 cm. cc 
the simple rounding of the edges of the piers reduced the scour to 91 mm —. 
yhereas a triangular shape reducedit70mm. 
ct An aerodynamic ‘shape gave no further reduction, but the shape of a lens 
produced a reduced scour of 54 mm. (large values of p) 
The protection realized by a single pile before the lens-shaped pier re- 
duced the scour to 38 mm. This pile — — values of p and ee 
a reduction of the width-length ratio. 
A change in the ‘repartition of the velocity will have an ialladai on the 
_ gecondmember of Eq. 19 and,therefore, on the value of the scour. A higher 
aa roughness" of the bottom on a distance before the pier (with pebbles on the - 
bottom) is realized, In Eq. 19, was reduced whereas v was increased 
with the conclusion that the erosion had to be increased. That is what was os 


s @€3—-- found with an erosion of vale mm with the lens- shaped pier (54 mm with the 


The consideration of Eq. 19 shows that it would be possible to ‘consider- 


ably reduce the erosion by adapting the ee of the radius p with that e 


corresponding with the values of ‘The result was a non- 
¥ pier with a larger base. The erosion oe this new model was was sis gear il 
It is easy to see that the second curvature 09 P92 will | give a 
Ea. 19) with a negative second member. A rising motion will, therefore, _ 
; — the first diving motion, This was observed on all the models. For a 
rectangular pier, the small radii of curvature around the upstream edges 
ies. Behind the pier, with the lens-shaped section, the curvatures of the type 
2 was followed by another curvature of type i, and the formation of a sec- : a 
smaller scour was to be observed just behind the pier. 


Many y other results” could be deduced from three first ent 


s also from the consideration of the spiral motion induced by the diving motion = 
‘tine w by the reaction of the bottom. — results can be found in some — 
ti 


eau, 
1988, 
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| 
— © MOUON OF SeaMMents Gerivations, 
sffoui de piies de ponts en riviere,” Acad 
4 ison, L. J., “Affouillement autour dep de ponts 5° Serie XXII]. 1937, 
jue, Bulletin de la Classe des Sciences, Ser 
ire sur l'Erosion dan les cours” 
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letin Société belge — — 
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bras naturels ou artificiels,” Association internationale des 
res de Navigation, Congres de Rome, 1953. 
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CHITALE. 


of model and pro- 
3 totype data indicated that the depth of scour ‘could be treated as simply another _ 


_ length and that equilibrium depth of scour obtained in the field. This a 7 


of equilibrium depth is due to the author | and more clearly enunciated by the 


2 
ee least a as a first approximation the equilibrium scour depth, with 
= qualifications as to the flow conditions, appears to bea function — 
ie of ‘geometry, i.e. the relative | depth of flow, the shape of the pier and e 
the angle of attack. . . velocity of flow and sediment size (and, therefore, — 
rate of transport and intensity « of boundary ry shear) d do not influence the = 


” 


Some investigations in models have been made on the subject of scour ont : 
at Pporrsss piers inIndiaat the Central Water and Power Research Station, Poona 
_ Which sachet thought would be of interest in context with the author’s find- 
me The first series of experiments were conducted i in connection with the Hard- 


_inge Bridge over Ganga River and were reported in the Annual Reports of the _ 


Station onfor the years 1938 to 1942. A geometrically similar replica of Hardinge 
Bridge pier was embedded in Ganga sand of mean diameter of 0.29 mm in a 
parallel sided channel. The results of these experiments gave the _— 


=" in Seog = maximum depth of scour below H.F.L.; d = depth of flow in the © 

a ™ flume; b = widthof pier; and Ic = discharge per foot run upstream of the piers. 

7 _ This relation is not dimensionally correct and cannot therefore be — for 

Basic experiments were subsequently conducted in 1941 with the object of 
testing the influence of upstream depth and sand diameter on scour round piers. 

The pier tested in these experiments was also 1/65 scale model of that of the 


Hardinge Bridge. It was es in section, of length 1 ft, width 0. 6 ft and a 


31 Research of fficer, ‘Flood Control Central Water and P 
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2 “Scour Around Bridge Piers 4, May, 1956, 
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_EXPERIMENTS. 


V= MEAN VELOCITY U/S OF 
pens és = MEAN DEPTH OF FLOW IN THE CHANNEL. UIS OF PIER] 


DEPTH 1 OF ScOUR at | OF 
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was laid with sand of 0.32 the following used just 
around the pier in successio 

e of san ean size, inmm 
4 screened dam sand 


= 1, 
‘ : 
‘The the piers was laid with | upstream bed level. 


at cfs per ft was run and water level was adjusted to a 


diameter 
sand 
around 


69 
80 
87 
01 
18 
.69 
.80 
.95 
18 
.54 
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0.3012 
0.4762 


0.1261 
0.1631 


0.2254 
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mm 


0,103 | 
0,304 
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In the 8 ft channel with sand of m= 0. 


prs experiments: 
“Seale of pier: 1/65 b= width of pier = 0.57 ft = 37 ft 
1 = length of pier = = 63 ft 


constant q=1 cfs per ft 


ds = maximum depth | of scour =: Se the pier 


a particular depth, the depths varying from 0.5 to1.45 ft. Each « 
continued until final maximum scour was obtained round the pier. 
oo In afew tests in which the upstream depth was lessthan stable depth the up- By 
‘stream bed scoured and blanketed the pier. In such cases, 


— 
— 
#8 
— 0563 | v 
| 0.647 | sand 4 
— | 0.9700 
0.1005 0 | 0.069 | 
0.1631 O 0.428) nila 
| 0.52 | 0.417 | 0.16 ¢ 
3 | 0.65 | 0.68 | 1, 
— — 


the maximum depth of scour rat the : nose was just before | depositjon 
inthe scour hole of sandfrom upstream occurred. 
: Te In experiments when sand around pier was coarser than the bed material ee 


| me the bed around the pier was laid higher than upstream level, to get p 
scour round the pier for the upstream depth laid. 
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REFERENCE = STABLE STABLE DEPTH 


Cdia=o.16mm) 


o.- os 
Cdia=0-6am.m.) 
Simm.) 0-77 
Dern OF SCOUR AT 


WITH NALA SAND.... 


SCAL SCALE PIER ds: FOR q=1 


able 2 shows the results obtained as 


nl With axial flow, maximum depth of scour | was always at the nose of the 
_— 


r at sides being dese} by dd to 15%. 
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2. ratio of scour the nose and depth of flow in the 
simple relation with the approach velocity in the channel (see Fig. 14). 
4 3. The depth of flow on the onan) has aes an influence on scour at the — 


MA 


REFERENCE, 


OF SAND AROUND pr 


SCOUR BELOW NOR 
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DEPTH OF 


MODEL PIER. 


better criterion. n. Fig. 16 shows plot of depth scour the Froud 


| 
— 
= | + — 
| 
— 
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15) 4 t annear to he satisfactorv. The writer therefore further _ 


ber. scatter of is in this plot, th the 


is remarkable, statistical equation obtained over the range of experimental Ae 


ao is thus seen that experiments at the Research Station do not lend support 
% to the author’s equilibrium scour theory ‘y, and it appears that further rinvestiga- 
_ tions both in the field and in the laboratories are desirable before it can be 
Grateful acknowledgment is made of the kind permission given by M.G. | 
‘Hiranandani, Director, Central Water and Power Research Station to refer to. 
q experiments previously conducted at the at the § Station and also for useful suggestions 3 
offered by him, 
A. RYLANDS THOMAS, 32 F. ASCE, — Fig. 17, the author’ s design curve 
ali piers (Fig. 9), witha factor of 0.9 for application to piers with semicir- 
- cular nose form, is compared with results of experiments with models of piers 
carried out by the writer in association with Sir Claude Inglis. a0,08 These ; 
piers were models of the piers of the Hardinge Bridge over the Ganges River — 
— (now in East Pakistan) which were 37 ft wide and 63 ft long, including the semi- 
~ circular nose and tail. The model piers were fixed ina parallel- sided channel . 
a a bed of incoherent Ganges sand about 0.3 mm mean diameter and the pro- 
- cedure was to runa constant discharge without sediment load until scour had 
i. It will be seen trom | Fig. iM that the author’s design curve, and still more 
his curve from the Iowa data, Fig. lie well above the Poona results. 
med appears to be too great to be explained by sediment load, 2 tages 
_ this cannot be ruled out as the upstream depth may | be reduced more by sedi- 
- ment load than is the depth of scour at the pier. It would be of value if 
_ author presented the data on which he based his design curve, showing partic- 
--It is also most desirable to compare small-scale results with observations 
made under full-scale conditions. Such data are difficult to obtain because of é 
* the need to take observations during high floods which are very often of short 
duration. It is easier to measure depth of scour at the pier than in the chan- 
nel upstream, where a number of observations must be taken to average out 
the effect of bed waves. Even when this is done it is not certain that thebed — 


% thata reductionin channel depth, due for laecaie to resistence to scour, would — 
increase the depth of scour at the nose of the pier. 
_ The relationship between depth of scour at the nose of the | pier and the depth . . 
of channel upstream is, therefore, perhaps not the most suitable one for com-_ 
_ parison of full-scale data. . Nor is it generally the most suitable for practical _ 
f applicationof adesign formula, because the upstream depth during a maximum _ 
flood is not often known beforehandand would have to be calculated. An — 
_ in assessing this depth would lead toa corresponding a and Gromer error inesti- 


mating the level to which scour is liable to occur. 


33 Annual Reporte Central ety and Research Station, 


‘and nd Navigation Research Station, “Poona, India, 1949, 
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oe more basic relationship is that between ais of scour below water level, 
7 Dy: width of pier, b, and upstream discharge per unit width, q. The  relation-— = 
ft-sec units, scour sand within the limits of q q2 
PS at least from 1 to 6 (see Fig. (17) and not greatly affected by the grade of . 
sand, It is: necessary to ‘take into account the effect of obliquity 
flow, which may be done by using for b the projected width of pier ona plane > Cr 
: normal tot the di mid iy of of approach. A _A factor of safety should also be used in | 


Fig. 3 _ Design curve, “Bulletin 4 


4 (Semi circular cutwater, axial aan 
Has Ix 
Semi- -circular cutwaters pret 


Value of 


- 
= FIG, 17,—SCOUR AT BRIDGE PIERS 
mt ‘The ailehaiiai per unit width to bet used for design m must, in most cnses, ™ 
= than the mean Value because inrivers the 4 channel meanders from 


in cu pec, and f is the Lacey silt 1.8 the where 


m = mean diameter of bed sand in mm). Eq. 23, derived from data of several 
bridges in India and Pakistan, takes into account obliquity of approach and con-_ 
“a centration of discharge but not the width of pier, which is clearly an important =) 
factor. _ The data related to normal widths of pier and todepths of scour ranging — 
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MUSHTAQ AHMAD. problem of correct estimation of depth, 


_ and extent of localized scour is important from the point of view of establishing * ft 
sound design practice a and safety of hydraulic structures such as bridge piers, < 
7 _ abutments, spur dikes, groins, or pitched islands. This problem is much more ci 9 
important for hydraulic structures in alluvial rivers of West Pakistan 
Lj fine sand is found for hundreds of feet in depth, and scour depths of 40 ft to 80 “a 
r ft below water level are common. The author has, therefore, dealt with a sub- _— 
of special importance and great utility to this region, 
™ §_—iIt is proposed to discuss the author’s approach to the problem with special — 
in reference to field and laboratory experience in West Pakistan. The author me 
assumed the depth of scour as another length and has related it to the normal 
: depth or the width of a pier, and maintains that the depth of scour does not e 
pend on the degree of concentration until scour holes around neighboring ob- 
_ structions s overlap. He maintains that there is equilibrium or limiting depth — 
and believes that for a given mode of sediment movement, the depth of — e 
depends only on the geometry of contraction and the approach depth. He holds af] 
that the effect of velocity and sediment size can be neglected as that of secondary i 
order. _ The writer agrees that there is an equilibrium and limiting depth of | 
: ae and that the effect of sediment size can be neglected. Tests made on 7 
depth of localized scour at spur dikes36 also showed that localized scour depth 
does not vary with grain size in the range (0.1 m to 0.7 m) usually met with 
alluvial plains of West Pakistan. However, ‘concept may not be valid 


for the entire ‘range of bed material sizes ranging from fine sand to gravel an and 


_ The author’s view that scour depth is another length connected with normal = 
of pier that it varies with the or more cor- 


_ where the ooaag varies between 0.64 to 0.69, or on the average the relation a 


2 


Dir., Research Inst.,West Pakistan, 
fing 36 «Experiments on Design and Behaviour of Spur Dykes,” by Mushtaq Ahm 


‘ceedings, Minnesota Internatl.  Hydr. 1-4, 1953. 
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DISCUSSION 


in which K may ‘bea a of geometry the ‘contraction at the 
bridge, abutment shape and thickness and shape of the nose of the piers, ,and 
a so forth. Because the depth is a dependent variable on discharge intensity, <a 
td latter may be used in preference to ‘depth. Here only lies the difference in 
7 proach between the author and the writer in studying this problem. As | 
the two are not very different. The writer prefers to 
= K as a function of f boundary geometry, shape 
a pier 1 nose, and abutment, characteristics of bed ‘material, and distribution of a 
velocity in the cross section at the piers representing the concentration of — 4 


flow. The functional relation for for study may be of the type: 


~~ width upstream of the bridge; B; is distance between the abutments; n 
the number of piers,; andt is thickness of each piers. 
The writer does not agree that the depthof scour does not dependon the de- 
gree of concentration. It has been found that by changing the concentration or 
distribution at bridge site bytraining works upstream, the scour depth 
(can be considerably reduced. The variationof scour depthas a result of change 
ms in velocity and flow distribution due to the training works is explained in Figs. 


Ss 19, and 20. , Figs. 19 and 20 depict a a | model of the Ling Stream showing = &g | 


. at the bridge me to —- more uniform velocity and lesser scour as shown 
: In the method of plot adopted by the author in which scour depth in relation a a 
to U/S depth has been studied, the non-dimensional scour perimeter may not a 
_ have significant relation with velocity, for, in open channels as the velocity or Y 
discharge increases 80 does: the depth, although not necessarily by erosion of 
- the bed. For a uniform approach velocity, a mean depth can be assumed for _ 
use in the author’s relations. But, generally, abnormal scour depths at piers 
a or abutments are associated with concentrationof flow resulting from the var-_ 
-* iation of velocity and depth in the approach section and the selection of repre- — 
sentative approach depth for use in the author’ 's relation and the estimation of 
Z ’ discharge distribution in the channel and on the berm becomes very difficult. ee 
a In fact, the author’ 8 term (Qc + Qo)/Qe in Eq. 9 is, also, a measure of con- — 


centration of flow. However, the estimation of Re | or Qw difficult for 
practical use in the computation of maximum scour depth in case of a curved 
‘approach. The writer has shown?” that change in flow concentration above a 
= dike, can be depicted in terms of | velocity distribution and common types 

are depicted in Fig. 21 by I, Iland II. The type III velocity distribution gives 
‘the maximum concentration onthe outside of abend as dueto negative velocity 
onthe inside, a part of the water way is blocked by reverse flow. For the three 7 a 

—— of velocity distribution depending on different types of approaches, the “ 


values of V and K as determined from s on spur dikes 
may be applicable to abutments), are given ‘in ‘Table 3. Since maximum scour 4 
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DISCUSSION 


by the river curvatures likely to occur, keeping in view the restraints aieed 
on river by the upstream training works. The method commonly used in West 


Scour below a severe bend on the ou concave side ac- 
-companied by a swirl on the convex bend. 


r Straight ot obstruction placed at any angle of 30° to 90° | 


Straight at an angle of 90° - 150° 


TABLE 4, SCOUR DEPTHS OBSERVED AT BRIDGE PIER AND 
ABUTMENTS ON DIFFERENT MODEL STUDIES AT IRRIGATION RESEARC i: 


Scour depth 


= 


“Maagla 
New 

S- 
Kasur 
Dhok P Pathan 


Indus: | Thatha 
Sujjawal 


proposed 
bridge, 


Computation from Fig. 9, Bulletin No.4. 
b Bed material - five sand. 


- ture likely to occur. - This will give scour depth below water level, , and the ‘depths 


a below bed level can then be worked out. Maximum depth of scour from model 
studies at bridge piers and abutments of rivers and 


— 
We... 
| 
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are summarized in Table 4, In this table, ‘scour depth has been n computed from 
the author’s design curves of Bulletin No. 4, in Fig. 9. ‘The scour depths 
tained after correcting for the model scale distortion are generally higher a 

those of the author. In computing Ds value by his method, the actual mean © 
value of .. ,has been taken from the cross section above the bridge. The value 
of K does ‘not exceed 1.7 because severe bends are not in 
to the presence of guide banks, = 


= 


42.9} 651.2. 


12.8 
5. 


oo 


bi, 
12 x 90 
= 1080 
15 x 90 
= 
1,52000] 691. 9} 14 x 90 | 
@ Part of the water way is marked by a semi erodable island at the bridge and hence 


e di fference in the water way for different years. * No training works a river bed u/s. a 
A spur constructed u/s to current approach, 


Fig. 22 shows abnormal scour ata guide b _— head. _ A heavy embayment is 

i noticeable on the right and construction due to silting caused by roller on the 

left. K value of 2.3 can be obtained in such a case, 
B a _ The heavy embayment resulting in abnormal scours illustrated in Fig. 22 


4 with K value greater than 2 are possible only on spur dikes or guide bank heads Be 
, 37 “Effect of Scale Distortion, Size of Model Bed Material and Time Scale on Geo- 


‘metrical Similarity of Localized Scour, ” by Mushtaq Ahmad, Proceedings, Internatl, 
Assoc, of Research, The he Hague, , Netherlands, 1955, 
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DISCUSSION volt ice 
ag and are ‘not to be allowed on bridge piers. or abutments. In fact, for alluvial 
. ‘viveta or streams of West Pakistan, guide banks at least equal to the length 
of the bridge, with curved | heads and expanding water wayon the upstream side, — 
are provided to shift the maximum ‘embayment and abnormal localized scour 
tothe guide bank heads instead of allowing it to occur near the bridge with 
short abutments. Under these conditions, the values of K for estimating the © 
maximum probable depth of scour can safely be taken be between 1, W! and 2, 0. In 
Sng = of ashort abutment length before the conditions for abnormal scour r where i 
- K greater than 2 are obtained, the approachroads or railway will be threatened 
sg cut by the embayment formed by an alluvial meandering river. It is there- 
for necessary to design pier and abutment depths for scour calculated from 
= 1.7 to 2.0 and ; provide proper - guide bar bank to keep the road or railway ap- 
proaches safe from embayment and to keep abnormal scour away from the 
main bridge crossing. Actual scour depths observed on the railway bridge on 
OR Ravi near Lahore and the computed values of K as recorded in eel 5 = 
supported the recorded values of K. iy 
_ PIER LUIGI ROMITA.,® 8_The author’ 8 long lasting efforts to cast light upon 
(the pleated of scour around bridge piers and abutments are to be — 
commended because of the great practical importance of the problem. ‘™ 
clear and condensed presentation in this paper of the results of these ctforts ; 
is a substantial step towards a generalized solution of the problem, and — : : 
of great use to designers and tothe administrations responsible for the con- — ; 
- ‘ struction and maintenance of roads and railways. To this enda coordinated — 
7 effort should be made, in order to aaeued the necessary verification on proto- a 
_ There is no doubt that one of the most important factors of scour, when — a 
piers have a sufficiently high length per width ratio, is the angle of attack be- 
: tween the pier and the flow. Even small deviations of this angle from the zero 
value are responsible for fast increases of the scour; ‘this usually represents 
a much greater danger to the stability of the pier than any underevaluation of — 
the maximum possible flood. The collapse of many piers during floods is, in 4 


cross- currents and deformations of the river bed, than to the unexpected 
In view of all this, the writer carried out, some years ago, a systematic — | 
del investigation o of the influence of the angle of attack upon the depth of 
scour, The experiments were carried out in a glass-flume of the Hydraulic 
_ Laboratory of the Polytechnic Institute of Milan, supplied with clear water, 


_ there was no general b bed movement but only localized scour around the pier, 
‘The shape of the pier was not particularly studied, and it reproduced a rather 
widely used type its length per width ratio was about 5. Two series of 


‘The data obtained at a zero angle o: of attack are in rather good agreement _ 
with the curve representing Eq. 11 Fig. 9, but for the fact that they consistently _ 
lie a little below this 3 curve. . An increase by 60% in the discharge s (which cor- 
responded toan increase by 20% in the average velocity of flow) brought around a 
only a 20% increase in the maximum —_ of peur, at the same angle of aha 


| 
| 
— 
hat 
__uging uniform sand of 1 mm diam as bed material, and in such conditi = 
4 
4 


The tests with varying angle « of attack have shown the basic importance of 
R. this factor. For an angle of is the maximum depth of scour was 1.9 times 
x that corresponding to a zero angle; for an angle | of 30° the increase of scour 
_ depth was in the ratio 2.6:1 in respect to the zero angle, for an angle of ll 
£ the ratio 3.0:1,and for an angle of 60° in the ratio 3.3:1. Increases of the angle 
of attack beyond 60° and up to 90° did not bring around any further er appreciable 


iG increase in the scour depth. These values are considerably - in excess of those - 
_ Another interesting ‘observation was that the point of maximum scour depth — 


iy aways: occurs 8 very | close tc to the pier wall, so that, af the depth o! of the foundat n n 


long stretch upstream of every bridge crossing. 


the possibility that an angle of attack occurs should always be taken into ac- _ 
count, and the pier foundations accordingly. 


the river banks should be stabilized with adequate measures fora easy 


In braided rivers, however, 
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As a conclusion, the fundamental importance of avoiding any angle of attack 
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‘CONSERVANCY DISTRICTS AS FLOOD CONTROL ORGANIZATIONS? 


Discussion by Samuel A. Greeley 


5 

“SAMUEL A. GREELEY, F. ASCE. —This is an and timely paper. 

In my opinion, it has important bearing not only on 1 conservancy districts, but 


"portance ‘of organizations of this kind, which may be generally referred to as 
Authorities, was emphasized in fine fashion by Mr. Robert Moses i in his address 


: before the American Public Works Congress in New York, on August 17, 1960. 


It is pertinent, therefore, to consider Mr. Chambers’ — as applying to more 


kinds of public works activities than flood control. 


It is of interest to refer to the Progress Report dated January, 1945 of a 
Committee of the Society on Organization, Financing, and Administration of 
- Sanitary Districts, of which the writer was Chairman. Relative to the organ- 


ization and creation of sanitary ¢ districts (or authorities), the following is ab- — 


er In general, the Committee feels that the steps required ieee form, 7 
and organize a sanitary district should be taken wi within the locality which 


et - The Committee favors the origination and creationof districts by local a 


petition followed by ageneral vote by qualified persons withinthe desig- — 


‘The Committee favors the determination of the extent of the district 
by: a “competent engineer or board of engineers after an n adequate — 
and investigation of the locality, 
oe e Committee favors, in general, the election loc of t e memt ers 
- Thec ittee f ing al, the election locally of th ieee 
of the administrative board. The Committee recognizes, however, that 
many appointed boards have been eminently successful and suggeststhat 

ll appointment of board members rs by some locally elected official may ae 

In Item 3 of the ‘Summary, Mr. Chambers’ outen to the needs « of “ “such de- 


y © need as compared to other interests should be included in the original studies. 


-.. is further suggested that the paper does not give sufficient weight to fi- 7 


nancing ‘procedures. In two | elections under the Illinois Conservation District 


April, 1960, by Cloyde C,Chambers, 
Partner, Greeley and Hansen, Chicago, Tl. ‘ts 


Act, organizations of the districts have failed. The writer thinks one reason 


_ was that the enabling act limits all financing to taxes derived from assessed — 


valuations. It seems likely that properties in different parts of the district | 
will have different views relative to the need and validity of the proposed dis- a 
trict or authority to be organized and financed. he ottaining of funds the 
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i = charges or assessments by zones carefully related to the resulting benefits an 
or uses of the flood control works. This is a good example of an approach to : 
j fair financing. - Having in mind the conception that the paper as a whole is use- > ; ; 
= ful in the consideration of other kinds of municipal organizations such as au- 
_thorities and districts, reference should be raade to the so-called two-part 
po rate as developed by a Joint Committee of the Society and the American a 
- Association with six other National organizations. The report of this group» 
_ stated the following fundamental principle relative to afair procedure in finan-_ 


- 


7 “The needed total annual revenue of a water or sewage ‘aaa shall be i . : 
eontributed by users and non-users (or by users and properties) for 

_ whose use, need and benefit the facilities of the works are provided ap- = : 
ae proximately i in proportion to. the cost of providing the sue and the b ~scal , 


4 _ This statement means that the needed total annual revenue (operating “s 


debt service costs) shall be derived from the users- and in 
tion to the amount which to be spent. 
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COMPARISON OF STREAM VELOCITY /METERS* 
Discussion by Harold G. Golden and I. Trotter, 
HAROLD G. GOLDEN,!° A. M. ASCE and I. L. ‘TROTTER, this 
yy interesting paper, the authors concluded that the screw type andcup type > 
= meters gave identical results in U.S, Lake Survey flow aie 


‘PABLE 4,—COMPARISON OF DISCHARGE MEASUREMENTS 
AINED BY PRICE AND OTT METERS 


varia 
tion of discharge ie 

obtained by Ott 


Ott meter 


discharge obtain- 


9-5 
6-30-58 


7-21-58 

4-58 7 2 
8-25-58 624,000 "624,000 

9-26-58 399,000 | 397,000 

11-10-58 | 206,000 | 205 | 


percentage variation of 1958 measurements — -0.69 


4-18-60 1,047 000 1,059 ,000 +1, 4 
600 546,000 547 


PR 


‘Average percentage variation of all measurements 


April, 1960, by F . Wayne Townsend and F, A. Blust, 


10 Hydr, Engr., U. 8. Geol, Survey, Jackson, Miss. 
11 Engr. Tech,, Hydr., U. S. Geol, Survey, Jackson, Miss, 
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Ott current meters at the stream- -gaging station « on the Mississippi River 1 near 

= Miss. The results of these measurements are shown in Table 4. - 
. All measurements were made from the downstream side of the highway bridge _ 7 


near Vicksburg. ‘In all cases, , the Ott meter discharge measurements were — 


made after the completion of the Price ‘meter discharge measurements, with 
7 an average time lag of 


4+ hr. The cross sectionalareas used forthe Ott meter > 
- measurements were based on the soundings obtained during the Price meter 

measurements. For comparative purposes, The Ott meter discharge meas- 
urements were adjusted when a significant change in river stage be- 


tween the two measurements. 
An equipment change was made between the measurements in 1958 
and in 1960. The measurements in 1958 were obtained using a 300 lb Vicks-— 
Baty -type lead sounding weight and a single pin attaching the Ott meter to the — 
- hanging-bar. The measurements in 1960 were obtained ; a 150 lb Columbus- se 
type lead _ Sounding weight and two pins attaching | the Ott meter to the a 


bar -. The ¢ effect of this — change is not evaluated but is notedas a mat- _ 


— November, 1960 — 
Of particular interest eter, 
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teat LOSSES IN LINES WITH SERVICE CONNECTIONS 


Discussion by M. Diskin, Annabel L. Tong, M. ‘McPherson 


and J. ve and K. M. Yao 


MH 1. DISKIN. author has presented a solution, based the Hazen- 
— formula, for the problem of head losses ina pipe with a non- -constant 

7 discharge along its length. The solution is based on a number of assumptions 
_ that should be critically examined before the solution can be adopted for gen- 
eral main assumptions on n which the analysis is based are: 


an The head losses per unit of pipe, at any the 
‘same whether it is a constant i tal or whether it is variable along the 
2 The local head» ees in the main stream due to the side outlets are 7 
¢ _ 3. There is no redistribution of energy content, per unit weight of fluid, 
between the main flow and the side flow. 
‘ e first two ‘assumptions are _ mentioned or implied in the orginal paper 
and require no further explanations. The third assumption refers to the pos- ‘ 
sibility that water arriving at a side outlet (Fig. 8) with a total head (energy > 
per unit weight) Hy and discharge Q, will be « | divided into two streams Qj and 
Q3 in such a way that one of them will have a higher head (He > Hy) and the 
other a lower head (Hg < Hy) than the original head ( ( Hy). » The total energy | 
‘halance will, of still be in with the he principle of conservation 


in which y is the unit weight of the liquid. 
: That sucha redistribution is possible can be demonstrated witha a simplified — 
BF eet of the conditions at a side outlet using the momentum equation. It is 
_ assumed for simplicity that the side outlet is in a horizontal plane. The free — 
body to which the momentum equation will be applied is defined by sections | 
bs 1-1, 2-2, and 3-3 in Fig. 8. Section 1-1 is taken just before the side outlet, i 
and the other two sections are taken sufficiently downstream for the velocities | 
to be uniform. The momentum equations for the free body 


a —e Head of Hydr, Lab., Technion, Israel Inst. of Tech., Haifa, Israel. 


alk 


— 
2 

| 
and based on the usual assumptions for the momentum equation, A 


-@ 


& 


¥! Dividing Eq. 22 b by Y Ay, it it is is transformed aa 
which may also be written as: 
: Since the ie at section 2-2 is smaller thanthat in section 1-1 (V9 < Mi) 
- the total head at section 2- -2 is higher than the total head at section 1- 1: : 4 


— 
bie 
sin which Ay and Ag are cross sectional areas of main pipe and side outlet : 
‘spectively, pj, pg, pg are pressures; V1, V2, V3 are velocities at the corre- 
sponding sections and pq is the pressure at the main pipe walloppositethe side 
— 
— 
a 


Dividing the se second mo. momentum Eq Eq. 23 by y Y and assuming that the pres- 
sure at the pipe wall opposite side outlet (P4 is equal 


A 


from which see is seen 


at is aie claimed that the preceding analysis is exact; S 
rather to show that a redistribution of energy at a side outlet is possible and 
that it could influence the results. Evidence that such a distribution takes place 
a can be seen in the experimental results reported by Jaeger® for head loss co- 
- efficients at side outlets that, at cortain combinations of of flows, beeaeesal negative — 
al Consideration of the possibility of the energy redistribution, together with | 
possible errors introduced by the other assumptions given previously, leads to | 
a * conclusion that more experimental work is desired on the subject of head = 
losses in pipes with varying discharge, before the results presented in the pa- 
can be accepted as the actual conditions in such pipes. 


ANNABEL L. TONG,’ M. ASCE. - —The author’s methods in determining the 
friction losses in pipelines with service connections will give exact solutions 
for single pipelines of uniform size and constant coefficient of friction in the 
_ pipe. In Tables 2 and 3, equivalent discharges 28 Qe are given for ideal cases of 
inifinite and finite numbers of connections versus (4/2) ‘andalso versus num- 
ber of connections “n” for the case of finite number of connections. But, in 
actual practice, water » network systems in service are mostly constructed of 
: enclosed loops and sometimes combined with a number of dead end lines. It ’ 
would be interesting to know to what extent the author’s methods are applicable 
, In order to make use of his tables, one may ‘convert a single loop system 
. ates of one source of Supply into a ‘Single line by taking an equivalent — 
= Le to represent the complete system asa “unit. When the pipelines in 
one system were built at the same period in an area of similar environment, 
their physical condition would be assumed alike and a constant c value would 


pp 480, Mechanics by C. Blackie & Son Ltd., London, England, 
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be used for good accuracy. With these se sim-_ 
_ Plifications, the author’s tables can be readily applied to a simple loop system | 
An example is shown in Fig. 9 to illustrate the degree of accuracy of his — 
‘eee applied to the loop systems. A water distribution network of 12 pipe-— 
lines, forming four enclosed loops is shown and the size and length of each line 
are marked in parenthesis adjacent to the line. ie ‘frictional coefficient C of © 
- 100 is assumed for the entire system. Seven service connections and the cor- 
_ _ responding discharge rates are indicated in junction points B, C, D, E, G, H, 
and 


I. if a | hydrant test flow of.300 gpm (21) was make at junction si and the _ 


| 
head loss (Hy) through the system during this flow is measured to be 4 ft, -_ 
what will be the corresponding head loss (H2) when the hydrant flow increases aa 
_ With the information obtained from the hydrant test the head loss loss” Hg) may 


a . Method 1.—The system is ; assumed to be a ‘single line of infinite number of — 
connections where q = 700 gpm and = 


(a) gpm. From Table 2, with 

q/Q; = 0.7 and Qe /Q = 0. -676, one may find: (Q,) = 676 gpm. To determine — , 

gy of athe system, the Hazen-Williams formula in a simplified form, may be 
. 


“Nomograph For Equivalent Lengths of Water Pipelines, by . L. Tong, Wa- 
ter and Sewage Works 1960, 
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where Lg is the length (in of the pipe for the 
unit for Qe is gpm and the unit for H feet. 


When Q; changes to 1700 gpm and from Table 2, (Qe) is com- 
Method 2.—If seven connections are assumed inthe equivalent | pipe, Table 3 
can be used to compute Qe. Using n= 7, q/Q = 0.7 then (Qe)1 = +134 x 1000 


gpm. The length is is computed uted aS foes sient 


Similiarly Q = 843 x x 1700 gpm = 1430 gpm, and 


‘Head loss Hg is then computed by the Hardy | Cross Method, , the 
answér is found to be 14.0 ft. Using this value as a basis for comparison, ‘the 
_ percentages of error from the preceding methods are approximately equal to — 
 +6% and -2%. If the number of service connections reduces to 3 by discharging © = 
200 gpm at junctions H and D and 300 gpm at junction G only, the result | by 
Method 1 is unchanged and by Method 2, Hg equals to 13.0 ft. Again using the 
' Hardy Cross Method, the head loss Hg corresponding to 3 service connections 
is found to be 15 ft. The degrees of accuracy by Methods 1 and 2 are -2%and | 
--15%. A change of location and rate of service discharges shall effect the value 
S of the equivalent length for the system slightly, and subsequently shall effect — if 
_ the head loss under the same fire flow also. But both the methods are satis 
For a complex. system consisting of pumping g facilities and two or more 
«sources of water supply, the assumption of using an equivalent, length to rep- 
resent the system will not be ge the above should be 


B. McPHERSON, 9 M. ASCE AND J. V. RADZIUL°—Before c commenc-_ 


overall design factors is necessary to place the limited scope of the subject — 

United States metropolitan municipal water works, itis common practice 

to supply me Getrintion system, or separated service districts, via an arte- - 


prof, of Hydr, Engrg., Civ. Engrg. Dept., Univ, of [linois, Urbana, Ml. 
10 Water Dept., Philadelphia, Pa, 
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November, 1960, 


rial feeder) network, The smaller-pipe | network is connected to the 
arterial system at relatively few points (about 1000-1500 ft apart). In a strong © 
=— a greater number of grid connections is not nocnneery while a mini- 


ar 
connections are normally made into mains. In systems 
. oe the arterial system and grid system a are often inseparable or indis 
tinguishable. Analysis of large systems is normally concentrated on the arte- 
rial network. Some of the major design considerations for arterial network a 
Se been presented previously. 11,12 The author’s paper is con- = ’ “a 
cerned principally with grid “water mains.” In closing, the author states that __ 
_ his procedures are “particularly suitable for use in network analysis... _ a 
eg some of his procedures might be adapted to the skeletonizing of a grid 
network, they have no place in arterial network analyses. in, ee 


bes Only the more progressive water works maintain a program of field meas- 


‘urements of pipe friction loss coefficients. Usually time and fund limitations © 
restrict measurement programs to the arterial network; only scattered sam- ; 
pling of grid network mains can be made. In t the city of Philadelphia, for ex- 
- ample, there are over 2,600 miles of mains equal to or smaller than 12-in 7 
_ dia, but only about 360 miles of mains larger than 12-in. With few exceptions, © 
the arterial mains are larger than 12-in. The auihor refers to the head loss | 
“which actually exists” as a point of reference in defining a “range of error, ” 
 . it is seldom that a head loss coefficient is known to any resonable degree of . 
exactitude for the grid mains with which the paper deals. In fact, one might — 
ask how a flow test can be — properly in the presence of active ser- 
Let m represent the p power to the flow is raised to compute head 
loss. The author has assumed m = 1.85 throughout his paper. Arguments ad- is 
vanced i for the use of m = 2.00 for “old” pipes in arterial networks are not a 


12 with the little inform ation ‘normally available on grid 
> “mains it would seem quixotic to argue whether m = 1.85 is more suitable than 
2.00 for grid mains, and this discussion includes both powers in relations to- = 
? be presented. For the low Reynolds n numbers characterizing many grid mains, 
use of either 1.85 or 2.00 and assuming constancy of friction coefficient, re- 
_ gardless of Reynolds number magnitude and variation, , are both unrealistic. — 
Further, there are local losses fittings, valves, and specials) in grid 
mains than in arterial. Field flow coefficient values for arterial mains include _ 
the influence of all local losses over the reach of measurement . Although local 
losses in grid mains may be neglected, asain the paper, the consequent exacti- 
tude of the results is even further ‘dirainisned, 
While the paper implies that grid aesign is based on domestic in 
the United States, fire flow requirements are often more severe. It is oo . 
~ implied that the local domestic demands, q, are known as a precise fraction | 
_ of the imput to the local main, Q;, which is not likely. as Se, ae 

_ Retiring from the preceding reservations, the remainder of this discussion F 
will be directed towards the development of the equations in the paper. As in : 
_ the paper, a constant pipe coefficient (embodying friction coefficient and dia~ is 
meter) and equally spaced, equal-valued drawoffs wiil be assumed. 
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— of of Eq. 9, has as appeared in /perhaps a 
- dozen American hydraulics and fluid mechanics texts, but not in such general 
_ terms. One of the earl’ =r presentations was by Dougherty. 13 Fair and Geyer!4 © 
have used the same principle in describing the hydraulic losses in filter under- 
@rain manifolds. These two references use m = 2,00. The author refers to the 
yar of Eq. 9 as being for an “infinite number of connectiuns” which is 7 a 
ae? with the idealized uniformly decreasing flow (or velocity) assum-_ a 
* along with constant diameter and friction coefficient. When the flow = 
ne service connections is large and the number of services are few (or 4 
example, n less than 10, q/Qt more than one-half), performance is more com= 
plex. It t cannot be described properly by means of the simple pipe friction re- 
7 ~ lations | to which the paper is restricted because the hydraulic behavior of : 
‘manifold system15 is being approached. It is not clear why the series com- — 
prising Eqs. 4 and 5 were introduced because Eqs. 1 and 6 are readily inte- ‘ 
as shown in Eqs. 7 and 8. % 
In the section on a “finite number of connections” (n), there are typograph- 
vs ical omissions and an error in Eq, 12. Eq. 12 has been corrected and expressed 


in t K' i 


idealized case the number of pipe lengths between cross-mains must be one y 
plus the number of service locations (two services ‘served ; at each location; | 4 
there are thus six. service locations in Fig. 10(a)). In Fig. 3 the author has 
either one too many connections or one too few pipe lengths. Eqs 10 through ee 
& 14 and much of the remainder of the paper ane based upon this error. Letting © . 

_ the number of service locations be n' (that is a differentiation from a of the 

p paper), ‘the | corrected, general expression for K is given by Eq. 35. ‘, 


= 


‘Table 9 (a) lists’ of x’ for 35 n' and m = 85. Ea. 9 was, 
therefore, used for an infinite n ", Curves which should replace those in Fig.6 __ 
are given in Fig. 10(b). It should be noted that for a realistic number of ct 
_ vice locations, K' is practically constant, as opposed to the large variation ia 
indicated in Fig. 6 or Table 4. The most divergent case is for a q/Q of one Bs 

which would be for adead-end main or a main in which flow was supplied fom 


pot h ends . Eq. 35 is not proper for this case because one too many lengths 
u is the instance in which 3 and Eq. 34, would 


4 13 «Hydraulics,” by R. L, Dougherty, McGraw-H 


15 «Engineering Hydraulics,” » ‘Hunter 
1950, p. 436.0 
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Values from Eq. 35 for 1 m = =2. 002 appear in Table > 9(b) for comparison, Note 


= The most partinent section of the paper deals with “a “ratio removed at the — a 
beginning of | of the line,” r. ” r. Eq. 19is not affected by | the error in the development — oa 


i Eq. 34, and has been ex; expressed in terms of K' in Eq. 37. besten hm = 


38 is for case m = 2.00. 


ar in Table ees and from 
While — leads to a dependence of r on n, 
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Fig. 7 and 8 ae 10(a) an and | 10(b) indicate an 
_ proximate value of about 45% for most of the realistic oombiantions. The 
author says that “European practice t has been to assume.. 45%... , which ae 
pears reasonable. On the few occasions in which this type of has 
arisen, the writers have used a 50-50 division, , partly to allow for fiow from | a 
either direction but mostly because further refinement could not be justified. 
‘This practice, to the best knowledge of the writers , is common in water and 


as works practice; the > writers are not aware of any utility where “a single + 
juantity at the end of the line” is. used, 
a The assumptions” leading to Eq. 35 are little better than those leading to 0 
Fig. 10@), it may be seen that the assumption of equally ‘Spaced 
loads | involves a series of premises: building lots must be equal in frontage, -_ 
_ services must be located in the center of each lot, distance from center “ ~~ : 


“frontage. In the opinion of the writers ‘these are negligible 
4 pared with the assumption that the distributed among 


9.—VARIATION OF kK’ K' FOR EQ. 35 VERSUS NUMBER OF LOCATIONS n', Sat 


= ,No. of 


0.442 
| 


0,820 0,580 0.520 “500 


0.816 | 0,662 | 540 
| (0,656 
0, 813 0.653 


0. 00 
0.527 | 0.425 | 0352 
520 4 0. 0.520 | 0.414 333 


in the bargain, is known to an exactness consistent tent witht the numer- : 
ical values which issue fromitsuse. 
elegance is desired, a paper by Brown and Markell6 be consulted, 
They have mathematically treated the case of varying diameter, varying flow oF 
coefficient, unequal service flows, , unequal lengths between services, illustrated 4 
for an m of both 1.85 and 2.00. Their development leads to a two-load equiva- # 
lent. The only assumption for an m of 1.85 is “that the sum of a MacLaurin’ a 
series (of a specified form) canbe nearly approximated by the sum of the first 
three” terms of the series.” No| assumptions are involved for an m of 2. 00. a 
Be their development leads to apair of equivalent loads ; which give iden. 
tical losses regardless of direction of flow from one end or the other. Of - 


Cunmummmainan? by Wilmer T. Brown and Robert A, Markel, AIEE Miscellaneous Pa- 
per 52-334, September, 1952 (Presented at Toledo Meeting, October, 1952), ee ; 
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Mcllroy Network Analyzer, is germane to this discussion. Because his dis-_ 


" sented in the ‘subject paper and for conceiving that it might exist. The 
method should often permit studies of Pipeline n networks to be 
out. accurately on. linear or non-linear electric calculating boards that 
_ would otherwise have inadequate capacity, and it should be generally 
helpful in studies conducted on paper. The opportunity of 
all the loads and pipeline sections between two adjacent network k junc- 
tions by means of two pipeline elements, one intermediate load, and one oo —— 
terminal load is a very significant contribution to the art of pipeline- ae 
“The suitability of the method for flow in either direction, subject to 
Bo: iy: stated limitations, and for exponents of any value encountered in 7 
fluid- flow problems, gives it wide usefulness. 
4 = “In some studies of pipeline networks, exact values cannot be as- | 
a5 signed to loads at definite locations along the pipelines. In these in- ae 
tances, the assumption | is often made that the total load delivered to = 


i 


0 485 
0,490, ag 0.478 0.462 0.445 
0,490 482 470 0,456 


485 468 440 | 0,412 


F users along the pipeline is uniformly « distributed between network junc- 
tions and may be represented by two half loads at the ends of the meee 
. tion. This ) rough equivalence is false, asthe authors point out, but under 
_ many common operating conditions causes an error no greater than that — : 
_ inherent in the nature of known data. The error is not significant when- im * 
“s ever the f flow delivered to ot to other network elements at the downstream tial 
> terminal of the pipeline equals or exceeds the total load distributed ica - By. 
along it. Wherever reasonable doubt exists regarding the ratio of dis- _ 
_ tributed load to terminal load, or wherever definite values are known for i 
: toads” at specific locations along a pipeline, the authors’ method is — 
powerful and accurate tool for reducing the number of quantities to be | 


considered in an analysis. The procedure is a technical contribution of at 


_ _ While some of MclIlroy’s statements would appear to refute several 2 argu- 


the it should be noted that the 


— 
ks 
(a) Via Eq, 37 using K' from Table 9 (a), m=1.85 
i | osc 
0.293 
0.376 
— | 
— 
A 
— — 


rown and Markel was phrased in terms of gas network analysis. err ie 
tudes and variation in local gas loads are normally known to a much closer 4 
3 y tolerance than with water services (namely, the major share of the design load 
as aa often comprises the winter heating load which fluctuates rather predictably a 
with temperature). Mcllroy’ s discussion was, therefore, directed principally 

- In summary, the common practice ee splitting the service connection loads m 

: “ane between the two ends of a grid main is entirely consistent with the de- — “ 
aaa gree of accuracy to which other hydraulic factors can be evaluated. If grid we 
main connections to an arterial network cannot be consolidated in an aes 
and equitable manner, the procedure developed by Brown and Markel is ispre- 


ferred, but can be applied only between arterial main junctions. 
_ Some of the major design considerations for grid mains and networks in- 
modern residential developments have been reviewed by the writers. 

3 Fig. 2, the term Qe is incorrect; perhaps Q, was intended ee 

M. YAO, 18_the arrangement of service shown in Fig. 3. 

Pie ‘is a bit different from the actual pictures. Usually, pipelines in a pipenet er 

- from one street intersection to another and it is unlikely to have service _ "i 

nection at both ends with possible exception of dead end lines, generally not 
- considered in pipenet analysis. Therefore, one service connection on the pipe- — - 
line would mean a service connection at the midpoint of the line instead of at the +o 


end of the line as assumed by the author. If this is the case, Eq. 10 will be: 


j : q 
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Design Considerations,” by M. B. McPherson and J, V. Radaiul, 
AWWA, Vol. St, April, 1960, pp, 400-602, 
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‘SEDIMENT PRC PROBLEMS oF F THE LOWER COLORADO R ‘RIVER? 


ASCE, —It is suspected that the authors, like the writer, 


og acquired most of their knowledge of river behavior without the aid of formal — 
aad 
instruction, for the simple reason that, | even now, practically no colleges pro- 
a ay vide any for engineers. Accordingly, it is hoped that they will not misunder- — 
7 _ Stand when the writer expresses the opinion that the knowledge required to 
_ forsee and arrange for the general trendof the events they have described ex- 

Rf hay isted before Hoover and Parker Dams were built, but was not available to civil 
engineers generally; and he hopes they will agree with him that the absence 
North American college curricula of some form of instruction on regime 

behavior of rivers and canals is as deplorable as would be the absence of soil © 
mechanics. Actually, there is a little college activity in the matter already. 
A The University of Alberta has given river andcanal instruction at posuraniae 
level since 1949, and is now (1960) introducing it in the Hydraulics, Soils, and 
Municipal Engineering options at undergraduate level. The University of Sas- 
_katchewan has introduced a postgraduate course. The writer is aware of a few 


_ the July issue of the Engineering Institute of nae Journal” at the 1959 Con- 
_ gress of the international Association for Hydraulic Researchh = 
pe That article Sis relevant tothe present paper bi because it (1) gives six hypo- 
7 a ‘thetical typical situations in which civil engineers ‘could make major e! errors (in- 
eluding lack of foresight), because they lack the knowledge of specialized col- 
leagues and geologists, (2) draws attention to official figures of the cost of | 
_ sedimentationtroubles in the United States, (3) | discusses the magnitude of Dal 
rorsthat canbe made with the aid of modern construction equipment, (4) gives 
y = opinion that “ignorance of river behavior is not universal; the trouble > 
i - the it is the prerogative of the civil engineers who build and plan water pro- 
_ jects,” (5) argues that the only effective method of | removing this ignorance is <7 
to “offer civil engineering students the opportunity to study appropriately 
- signed courses in river engineering, just as they now study soil mechanics 
(6) discusses the need for suitable textbooks and refers to the writer’s own, ie 
@) emphasizes that the tind nature of river phenomena prevents any single 


April, 1960, by Whitney M. Borland and Carl R. Miller, 
prof, Civ. Engrg., Univ, of Alberta, and Pres., T. Blench and Assoc, Ltd., Cons, 
¢ 15 -o “River Engineering as a College Course for Civil Engineers,” by T. Blench, En- 


gineering Institute of Canada Journal, July, 1959 Bal 
ecleen Behavior of Canals and Rivers,” 


nada, Blench, Butterworth’s Scientific 
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time. A Canadian of Indian extraction “drew attention to the value of teaching 
Fare 


Gar 


person from expertness in terms of nothing but the experiences ot 


formation without which, of course, even the most cannot be 

d 4 i _ applied satisfactorily. The writer, therefore, would appreciate fu further eluci- 

@ 1. In view of the facts that (a) Fig. 5 shows the bed material sissies: size 
% - about 0.3 mm) to follow closely the logarithmic Gaussian distribution of — 

cle size that sedimentary petrologists seem to have found general for sandy 
sediments, 17 and that the writer himself has found in river bed sands; ws and 
(b) the bed ‘material is following this distribution even after several years of 
a -_retrogression has occurred, what is the factual or theoretical justification for =e 
ps believing differentiation will now occur as implied i in the statement “It t was de- 
_ termined that bed stability would occur when a6-in, armor layer of 4 to10mm 
- material was obtained.” Also, assuming that this differentiation of coarse ma- <3 
terial is possible, in general, cana computation and description of transport : 
tox mechanism be produced to show how 6 in. of it could differentiate out of the 
particular | sand in the depth estimated as yet to ‘be degraded? Does the All- 

American canal show such differentiation, and is its behavior relevant? 

Ras _ How were channel widths found in Fig. 7 from Leopold and Maddock’s 
Ss when they did not propose any design formulas there, but instead re- 5 
az corded relationships with very different constants found for very y different rivers 7 
% whose bed-material sizes were generally not observed? How were widthsfound _ 
oe by any of the methods when none of them concerned erodibility of banks? —- é 
oa is the width for Q = 15,000 cfs and 0.3 mm material about 2,500 ft when the 
All-American Canal and afigure shown elsewhere!9 for 20, 000 cfs in the 
« Yuma show that nature is satisfied with about 450 ft? 

Why was canalization adopted as ag solution to the ne difficulties in | degrading 

a reaches when the effect of straightening : a ‘meandering ‘channel is to reduce its 
regime slope radically and shorten its | 


length, thereby causing it to degrade 


oe 4, Is the canalized channel expected to ‘meander? If so, what amount of 
stone protection is foreseen for the sides? ive 
Is the proposed section in Fig. 9 related to the channel dimensions 
= oe served by Leopold and Wolman!? to have developed at Yuma since about 1940? 
ies _ Their work seems to show that the channel is incised at discharges as _ 6) 
1s 000 cfs and had acquired, for that discharge, a breadth of about 450 ft from 
1940 to 1950, and a somewhat fluctuating depth that changed from 11 ft in 1940, — 
to 13.5 ft in 1943, and back to about 10.5ftin 1949.00 
«6 . Considering that the section just mentioned is for an appreciable bed- 
' toad, obtained by dumping the total sediment of the All-American Canal into 
- river (whose discharge is probably less than that of the canal), would the 
sectionne needed for vanishing bed-load not be considerably narrower and 
what computations was “dominant discharge” decided? 


x 17 “Size Frequency Distribution of senate and the Normal phi Curve,” by W. .c 
Journal of Sedimentary Petrolo 1938, pp. 84-90. 
18 “Normal Size Distribution Found in i of River-bed Sand,” * by T. Blench, 
Civil En eering, February, 1952, (Answers to discussion, February, 1953), e< = 
: “The Hydraulic Geometry of Stream Channels and Some Physiographic Implica-_ 


‘tions, by Luna B, Leopold and Survey, Professional 
252, 1953, Fig. 28. 
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s 8. Dothe authors attribute eany significance: to the following res results from re- vi 
Be gime theory!6 taking a bed factor 1.05 for 0.3 mm sand with — — * 
charge, and a side- -factor of 0. 20 pebeeceondn cohesive ‘sides: 


(a) canal of 15,000 cts = 1/11,00 000. 
"slope of meandering canal approximately = 1/5, ) 
Regime breadth for 20,000 cfs = 310 
Regime depth for 20,000 cfs = 16 ft 

Note that slope is very insensitive to discharge : and side factor pro- 
portional to bed-factor. Breadth varies as the square root of the bed-factor a 
divided by the side-factor, and depth as the cube root of side-factor and in- > 
versely as the 2/3 power of bed-factor, so b and c are not intended to be more 

9. Why was the problem r not not solved by model originally? | The reasons fa- 

voring a model include (a) there would be no bed-load injection, (b) field cata 
: ss every kind are available, or easy to collect, (c) models using natural bed- a 


aerated at the ‘U.S. Waterways Experiments Station that the meandering of riv- 

ers, including degradation and aggradation, could be reproduced 

Bese! even artificial sediment. (f) Leopold has reproduced the details of braid-— 

ing by using a model with suitably chosen sand.2! The conditions seem remark- 
_ ably favorable for a quick accurate model solution, provided regime principles | 
were used to eliminate unnecessary trial-and- error with scales. he 

_ In the present state of river engineering, the writer feels that the oo 


+ of. quantitative field information in accessible and discussable form is of high % 
value, and he congratulates the authors on presenting so much on 


Laboratory Study of the Meanderin, Alluvial Rivers,” by J, F. Friedkin, U. 
8 Waterways Experiment Station, Vicksburg, Miss., 
*F, _ 21 “River Channel Patterns; Braided, Meandering and Straight,” by Luna B, Leopold 
; andl M. Gordon Wolman, U.S.G.S, Professional Paper 282-B, 1957. 
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S BACKWATER AND DROPDOWN CURVE TABLES* pe 


oie Discussion by Praxitelis A. 


pon 
PRAXITELIS 21 F, ASCE.- _The writer will 


“4 on the formulas and tables by G. Tolkmitt for the computation in the rere 


practice of backwater and dropdown curves and to give practical and theore- x 
tical supplementary information on this method. 
3 _ The Tolkmitt’s method is well known in many countries of Europe as i : 
as in Greece. On the other hand, the Tolmitt’s formulas and tables for the ui- ff 
form prismatic channel of broad parabolic cross section form are not 
found in n the | German texts. 22 
: “ae very usual and important in Hydraulic Engineering to o determine ow 
fc a at which, in a given artificial or natural channel, the flow takes place. 
On the other hand, it is well known that when a certain point of the channel is — 
-_ obstructed by a dam, sluice gates, weir, or - other construction the raising of 
_ the latter will cause the water to set back up the. stream. In this case, the 
water depth is increasing, anda newsurfaceisformed. = 
= me Usually, in the case of a backwater c curve, two questions arise: how much = 
the water is raising at a given distance upstream from the point of obstruction _ 
and what is the distance for which weg influence of the curve is felt. The ques- | . 
“ss in the case of dropdown curve, are analogous. = | aa 
- Formulas, based on the assumption of channels possessing regularity in ~ 
"shape of cross section and having a uniform and constant bed slope, have been — 
4 = Also, in these formulas the roughness of the lining is assumed 4 
to remain constant. the other hand, to the computations, 
= tables were prepared. 
Generally, in the case of a prismatic. channel, m of computing 
backwater and dropdown curves can be d divided = 
or semigraphical methods; and 


‘The methods of Riihmann, Bresse Schaffernak, ‘Tolkmitt, Bakmeteff, 
puit, Chow, Von n Seggern,23_ Mononobe, 24 Ming Lee, 29 and Francesco 


pr, of Civ. Engrg., Special Scientific Researcher, 
_ 22 “Hydraulics,” by P. A, Argyropoulos, Ed, of the Technical Chamber of Greece, a 
Pe 23 “Integrating the Equation of Non-Uniform Flow,” by M, E. Von Seggern, Proceed- _ 
s, ASCE, Vol, 75, No. IR 1, January, 1949. 
= 4 “Backwater and Dropdown Curves for Uniform Channels,” by Nagaho Mononobe, — 
"25 *Steady Gradually 1 Varied Flow in Uniform Channels on Mild ild Slopes,” Univ, of Il- 
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- Ramponi,26 are noted. The method of Flavio »Scano2? and Escoffier-Raytchine- | 
Chatelain28- are graphical methods, while the Silber’s method29 isa method 
usinganomographh 
laa On the other hand, in the case of an irregular natural stream with varying - 
slopes of bed and different cross sections along its length, the solution of a 
é _- of backwater and dropdown curves is more difficult and delicate. In 
' this case, except for the tedious trial methods, there are the Krivoshein 
method, 30 the Grimm31 method, the Erza method , a sin simple method of the 
writer, 32 the Leach diagram33 method and others. 
Also, the methods for computing backwater and dropdown curves for the 
7 =e channels of trapezoidal, parabolic, and other cross-sectional forms | 
can be divided in h two categories whether the effect of the change in kinetic | 
4 energy y owing the accelerated or delated flow of water is considered or not; 
Therefore, the methods of Bresse, Struckel, Kozeny, and others consider the 
effect of velocity- -head V 2/2 (in which, V is the mean velocity of flow and gs. 


the Schaffernak’s method, and others, do n not consider it. 
_ Tolkmitt, who was particularly interested in natural watercourses, assumed ‘ 

in his method a constant value of friction factor C throughout the whole 
range of depths. Tolkmitt chose as an “idealized” profile, a =e channel 
(Fig . 3). Let Tg be the top width and A, the cross sectional area. a2 


The raising z the width ' T, will 

4 and the wetted area: 


at On the other hand, the hydraulic radius: wes 
2/3 T (d+ z) _ 2 


— 


26 “Deter mination approchée de quendiont du remous dans les canaux 
laires ‘a faible pente,” Energie Electrice, December, 1951 ee ee 
“Procédé graphique rapide pout tracer la courbe de remous de 1’ecoulement per- 
dans les canaux prismatiques,” Energie Electrice, July, 1952, 
28 “Determination graphique de la ligne d’eau et calcul des remous,” La Houille 
Blapche, “Etude et trace des ecoulements permenents en canaux et riviéres,” Ed., Dunod, 
i? 30 * Méthode naturelle de calcul des courbes de remous, ” Rev, Universidad, National 
al _ 31 “Blackwater Slopes Above Dams,” Engineerin -News Record, Vol. 100, 1928 ext. 
Se 82 “Computation of the Backwater and Dropdown Curves in Natural Streams,” Tech- 2 
nical Records, Technical Chamber of Greece, No, 401-402, 1957. 
Methods of the Solution of Backwater Problems,” Engineering-News Records, 
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is the fall of the tees water per 


On the other the 


To 
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— In this case, the m 
— 


27 may take the form 


| 
_ The integration of Eq. 29 leads to Eq. 1. the arate to find the formula | 
or the dropdown curve, whose correct writing is: _ an 


(0) 


is s analogous. 
ot artificial channels or natural streams having other than parabolic form. In 
this case, the cross section of the channel must be deformated in oud ; 
parabolic profile of the same to the given channel top width (Fig. 5) and 
=! cross sectional area Ao. For instance, if the — a is of a regular 


"trapezoidal cross section (Fig. 5), 


3 
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Wins 


1. ‘channel, 98 ft wide and 4.6 ft in depth, has a 
slope of 1 in 2,000 ( = .0005). The discharge is 1,483 cfs. By constructinga 
e. weir, the original water level has been raised 2. 6 ft. Estimate (by the Tolk-— * 
: . mitt’s method): (1) ‘at which distance up: upstream from the | point of the obstruc- oe 
gz tion the water is raised 7 ft, and (2) what is the total —* of the backwater 
ash = 6,9 +2 ft and, d+z = 6,9 + 1,0 = 


1,38 38 and, = 1.14 

ren 


rom ‘the Tolkmitt’s table, for ‘ 1,38 f = (Xx) =1 238 ote 


= £(1,14) = 0. 


235 - 0.8 818) = 


(1,001) = - 0,507 
aid 
hz) = (1,235 + .507) 
) -0005 = 


for computation of backwater and dropdown curves, ‘is 

valuable and satisfactory. European engineers feel a special confidence in 

_ this method, By Goodrich’s paper, the Tolkmitt’s formulas and tables, slightly be 
by the addition of twenty” items in each table many 


G 
| 
4 | 
— 
— 
hen, f 


-~HOOD INLETS FOR CONDUIT 
BERNARD L. GOLDING, !! M. ASCE.— Mr. Blaisdell and the en engineers of 
Research Service are to. be congratulated for recognizing ill 
potentialities of, fully developing, the hood inket culvert by 
Karr and Clayton at Oregon State College, in 1954. 
7 In the writer’s opinion, the very simple hood inlet is the answer aie the 
_ sc of making culverts. on steep (supercritical) slopes flow full (prime). — 
= Bey could result in the saving of many dollars c on the cost of culverts of the 


ae nesheenge. It is the intention ofthe writer to discuss this paper as far as) 
its application to the design of culverts under highway fills Ri 
: ss the determination of a culvert size under a highway fill, the headwater 
_ depth upstream from the culvert entrance and the exit velocity of the flow from — 
the culvert are the two principal factors that must be considered. The head- 4 
- water depth upstream from the culvert entrance must be limited sothatex- 
tensive flooding of upstream | property does not occur. The exit velocity of - 3 
flow must not be so great as to cause excessive scour of the stream bedand = fj 
_banks which could conceivably cause a washout of the highway embankment. 7 
Generally, up to the present ti time, engineers inthe writer’ s Office have felt _ 
egg to assume that a culvert would not flow full (an entrance control a. 7 
orifice condition existed) in computing the headwater depth caused by culverts | 
on hydraulically steep Gupercritical) slopes, even though the normal depth of 


=headwater depth was” greater” than 1.5 times the height of the. culvert. This 
ae “design assumption was based on the many observations and reports of others | 


as well as our own field and laboratory 
4 _ Making the assumption that the culvert would not flow full, resulted in the | 
a e headwater depth being greater than that which would have occurred if the same _ 
oy pipe was to flow full, or for a particular allowable headwater depth, that is — 
. id the usual design situation. The full flow condition would result in a smaller cul- — 


vert size than the entrance control condition, An example of the possible ies 
7 ita in culvert size that would result from the use of the hood inlet to insure 


full pipe flow willbe shown subsequently, 
_-‘The use of a smaller culvert flowing full will not necessarily mean that the © 
exit velocity will be any greater than in the case of the larger pipe size with 
entrance control. Generally, the cross sectional areas occupied by the flow in 
each case will probably be the same. In the past, engineers in the writer’ . 
have made the reasonable assumption that culverts with entrance 
‘eae 


— 
* 
— 
— 
| 
— 
aa ‘amen and Bergendoff, New York, N.Y. 


and so forth, at outlets based on this assumption 

io computing the headwater depth for culverts acting as orifices (entrance © 
= _ control condition), nomographs of the Bureau of Public Roads12 have been used. | 
hs have been reproduced in two publications.13,14 
In order to compute the head H over the invert of the culvert entrance . 


t the full flow condition, the writer has s simplified th the e Eqs. : 4, 5, ane 6 as follows: 


In deriving Eqs. 19a to c, the following simplifications and substitutions were 


The term (sin- 

ng Spl has been substituted for Z, where So is the slope of the culvert in 
ft per foot and 1 is the length of th the e culvert (not including the e length of the — 

Manning’ s roughness coefficient “n” has been substituted for the Weis- 

In order to further | simplify the ‘computation ofthe head H over the culvert en- a 
trance, a nomograph is given in Fig. 16 that enables an immediate solution of a 

In the solution of Eq. 19a, the reasonable assumption c can also be 


- duit at a distance above the invert equal to 0.75 times the diameter of the cul- 
F vert (B = 0.75) that enables the complete solution of the head H over the invert | - 
2 Generally, the standard state highway department headwall is equal to the © 7 

- Headwall Type Anti-Vortex Wall shown in Fig. 8(e), so that no other special 

k. device to inhibit vortices seems necessary. However, the writer would appreci- 
ate ‘more discussion author concerning the effectiveness of this type of 


~ When a conduit flows full, itis 3 usually a a good idea for the designer to com- | 

pute the subatmospheric pressure at the entrance of the conduit to make sure 

- ‘it is not less than the vapor pressure or cavation that: will occur at the ¢ 


de 
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“he in whic “4 
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a 
— 
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July, 1953 (Preliminary), 
_ 13 “Highway Engineering Handbook,” K. B, Wood, ed., McGraw-Hill Book Co.Inc., 
14 “Design of Small Dams,” Bur, of Reclam,, Gov. Printing Off., Washington, D. 
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D= Diometer of pipe in feet 

one Monning s roughness coefficient 
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pipe | entrance, When the conduit entrance is of the re- entrant type, such as the 


hood inlet, the computation of the subatmospheric at the entrance should never 
2 be ignored. The actual computation of the pressure at the conduit entrance is = 
“quite simple « once the head H at the entrance is known. From Fig. 17, it he be 


he + hr = hsa + 


in which he= entrance loss in feet = Ke 2¢ ; ; hr = reduction in pressure head 


2 in which K;= pressure 


‘ caused by the ‘contraction at = “Kr 


entrance, Eq. 20 may be rewritten as follows: oor, 
(x, +x, +Ky) =h (H D) . 


Aga = 
we 


Be a but for the re-entrant type of entrance, the pressure reduction coefficient oll 


— 
— 
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approach “unity pressure a 25 ft at mean sea 
is considered allowable in hydraulic design. 


‘forth as follows: 


Design flow 

Slope of 


* a 42 in. pipe flowing full will cause a head H = 7.9 ft over the invert of the cul-| 
vert entrance. The exit velocity from the culvert will be 18.2 fps. Substituting 

in Eq. 22, the subatmospheric pressure at the crown of the culvert entrance . 

9.0 ft that is well within the minimum subatmospheric of 25 ft allowed at mean 
- sea level, if a pressure reduction coefficient Ky = 0.8 is used in design. _ 
Entrance Control Condition. — (Square edged entrance) Eq. 23 may be used 

in lieu of a nomograph for orifice flow — edged concrete — enkeenees 

under low heads H>1.0D,H<3.0D) | 

eis in Eq. 23, a 54 in. . pipe will cause a head of 7.2 ft over the invert 4 

of the culvert entrance. If the depth of flow is - assummed to be two-thirds of © 


the pipe diameter at the outlet, the exit velocity of the flow will be 20.0 fps. 
th the northeastern states, the cost differential between a 42 in. reinforced a 


_ concrete pipe and a 54 in. reinforced concrete pipe in place is about ten (10) 
dollars per ft, that would result ina savings of $2,000 in the preceding typical _ 


: culvert installation. If the reasonable assumption of three such culverts in — < 
each mile of roadway is made, a cost anving of $6,000 per mile or $600, 000 
for each 100 miles of roadway would result. 
The writer earnestly hopes that state highway departments will tak tak 
cognizance of the advantages of the hood inlet as a culvert entrance, = 
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1, 6 line lines" from bottom Varwich” ‘should be “ Varwick” 
45, 12 from the 


bottom, the ‘to” is omitted “formu la” 


45 line and p page _IR should 
, Table I in headings, R.10-% should be R. 10° 
6 above whole line should be. 
p. 66, Chintu Lai” ‘should be Lai, Chintu” 


— 
7 
— 
— 
— 
7 


- 


=, 


OF) 


~ 


NOVEMBER 1960 — 


* 


DIVISION 


OF far foe 


HYD DRAULICS Jf 


OF THE HYDRAULICS ? 
OF THE AMER AMERICAN SOCIETY CIETY OF CIVIL EN ‘ENGINEERS 


fe iim 
| 
= 
| 
| 
— 
— 
— 
| — 
| 


Agee 


1960-39 


HYDRAULICS DIVISION 
7 Procee 


of the American Society of f Civil | 


PURPOSE OF = HYDRAULICS DIVISION 
“The advancement and dissemination of scientific and engineering knowl- 
edge in all branches of hydraulicss hydrology, hydraulic engineering and peas = 
_ Water resources. In particular this shall embrace meteorology and hy- 


a drology as 5 the sciences dealing with the occurrence of water i in the i 


added, and new chairmen take over. The memberships of the various rom- ‘a 
a es _mittees of of the he Division will be somewhat as f follows for the e 1960-61 Society 


Fred A. Bertle 


_A.T.Ippen = Water Resources Coordinator 


| Note. .—No. 1960-39 is part 2 of the copyrighted Journal of the “Tiydraulics Division, 


_ Proceedings of the American Society of Civil Engineers, Vol my? No. HY 9, Ss 


Copyright 1960 by the American Society of Civil Engineers. 
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Hydraulic Structures 


Pratt, Chairman 


W. Howe Be Alvin J. Peterka 


Joseph N. Rex A. Elder 


Clarence F. Wicker, Chairman 

Frederick R. Brown 


ns 
Eugene P. Fortson, Sr 
Francis G. Christian Ray K. Linsley, Chairman 
‘Robert B, Banks Arvi O, Waananen- 


Norbert 


‘Carl E. ‘Kindsvater, new member 


Es A. Koelzer 

"Joseph 1 Perrey, term extended 1961. Charles 

Ground Water H arolog Lindsay P; Disney, term 


Raymond Boucher 
= M. I. Rorabaugh 


_ Karl Emil Hilgard Hydraulic Prize 
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The technical committees hi many ac active task New task fore 
are formed as the need arises, , and old ones are ‘dismissed 2 as bevel _ 


Executive Committee Meets Soo 
of the Executive Committee was held on forencons of Thurs- 
x day, Friday, and Saturday, August 18, 19, and 20, 1960. The a. 
a 7 Chairman, P Professor Arthur T. Ippen, presided, © Those’ in at attendance wer 


Maurice L. Dickinson, Vice- chairman 
Carl E. ‘Kindsvater 
Mr, Eugene P. Fortson, Jr 
Herbert S. Riesbol, Water Resources Coordinator 

Harold M. Martin, Secretary 
Don P. 


‘The following are from the minutes of this: 


@ Publications. 


ma measured; also, that the chairman of the committee on publications should | 
&é _ receive all reviews, and assess them before sending them on the “ape. of 


iy 
_ Research. Dr. Vito A. Chairman, | reviewed the minutes of the 


“ amples of the kind of research which the Committee feels should be undertaken “al 


4 


oe *e One important action favored by the Committee is to raise the annual dues ~ 
of the Society about $5 and to reserve these additional funds for increasing © , 


_ search is important in the eyes of those most intimately acquainted with = 


cent needs. The effect of this would be to increase the availability of funds © yb 4 


‘This and other conclusions and recommendations were endorsed by the 
Executive Committee with the request that Dr. Vanoni prepare the report of 
_ the Hydraulics Division Committee on ee for reproduction in — 


“for research from both public and private = i 


— 
k force on Salt Water Intrusion (in ene Membership is 
— im 
— 
| 
e been received this year and 
e 50 papers hav ublication in tl 
— Som en is required for publica papers bymore 
— | site rating of seven is require hat rating of paper -— 
38 declined. Compose Erecut ive Committee tha 
bo z= nal, It is the feeling of the Exe — 
ar 
Se 


review and endorse research proposals was discussed. was 
_ that it should for the technical aspects in assistance to the Society Research a F 
Committee. can also advise an applicant on preparing proposals. prt) 
‘Sessions’ ‘Prog . The Committee prepared Hydraulics Division Pro- 
rams for th the e following conventions held during | the year 1960: Bae = 


New 3 sessions--11 papers 
sessions-- 9 papers 
‘Seattle Conference 6 sessions--23 


n 


: “past year. Four of the committees were activated during the year while —- 
ee on Flood Control had been active and productive for many years. Ncad : 
_ The committees on hydrometeorology, ground water hydrology, and surface | 


co =e hydrology were activated in late January while the Committee on Water 
‘Resources Planning was authorized in late March. The Committee on G-ound Bs 
Water Hydrology met in Denver with the Water Resources Coordinator on He 4 
February 19- -20, ‘and the Committees. on Hydrometeorology : and Surface Water 
By, Hydrology, and the Chairman of the Committee on Flood Control met with the 4 
t Coordinator in Chicago, Illinois , on April 21 and 22. The Committee on Flood 
Control held a meeting in Seattle on August 18 during the Hydraulics Division Rs 
Conference. _ With the exception of the Committee on Water Resources Plan- 
a ning, which got off to a late start, all the committees have been very active 
; alt during the year in the planning and execution of their activities. As a result, 


the following recommendations were approved by the Executive ‘Committee: 


~The following six-point outline of committee functions was approved: 


_ a. To seek out oan encourage the presentation of ‘papers pertinent to all Ae 
phases of objectives and purpose. 
To seek out areas of needed research and to promote and 


—_—sResearch Commit 
— 
3 ; oposed for 196 
Way ted a brief report on activities 4 
— l, Water Resources Coordinator, presente — 
— 
=" — 
— 
— 
— 


ASCE » "Hydraulics Division 

ee. c. To review ¢ and appraise on a continuing basis the availability of a. 
data in meteorology and hydrology and related physical data, as needed 

for the long-range planning of water resource development. 
= _d. To review and appraise standards of data compilation, analysis, 
a publication and distribution in terms of technical adequacy and — 3 , 


f. To focus of the on water 
bring them to the attention of the Society, 


2. The new ‘aeestinaia on pane resources in the 1e organizational structure of 
should enable the Society to become the dominant professional or- 
Sanization dealing with water resources in the United States. This would 
— auite properly focus attention on civil engineers as the professional lead- i &§ 
in this important area of national development. 
_ 3. The ASCE policy and program for publication of papers should be care- 
_ fully reexamined as it | relates to water r resources to to be certain that ‘it is - ae 
designed to give best service to the e profession. 
4 The water resources committees will sponsor 2-1/2 sessions a the 
Phoenix Convention, April 10-14, 1961. This was a reduction from the 
: four sessions originally proposed by the water resources committees. i 


5. The committees will not sponsor sessions at the New York Convention % 
an 6. The water resources . cominittees will sponsor three sessions (9 to 12 Pant 
papers) at the conference to be held in Urbana, Illinois, 
in August 1961. 
1, The Task Force on Spillway Design Floods will bidlienens its assignment 
as a task force of the Committee on Surface Water Hydrology. Mr . Joseph 
8. The Committee on Ground Water Hydrology will establish a Task Force 
— on Salt Water Intrusion. 1. The Task Force will carry out such a such activities ines a 


a ge be aware of the s status of salt water intrusion in th in the United tee 


b. To recognize sources and causes of intrusion. 
. To analyze methods for control of intrusion. 
d. To keep abreast of an encourage research on babies ™ oe 


e. To encourage presentation « of of papers relating to intrusion at t society 


= . The five water r resources Committees will, if possible, hold a joint 

"meeting in Chicago during the same week that the Executive Committee — 


‘The Hj Hydraulics Division has been 
7 _ invited to have its 1961 conference (Tenth National Conference) at Champaign- | 


Urbana, Illinois, August 16, 17, 18, 1961. Host will be the Central Illinois 


ASCE, and the University ‘of Illinois. Ww. J. , Roberts, represent- 
ing the Illinois local Committee, made announcements for the 1961 conference 


— 
— 
— 
ia 
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— 
— 
= 
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‘ond 


The Eleventh was for Cali- 
fornia, in 1962, since the Sacramento Section issued an invitation 2 years ago. a “5 
“i ae Division Conference. University of Washington, Seattle, August _ 
gistration of over 200; in addition 50 
i _ The local committee arranged excellent housing facilities and the social 
program was attractive. Two busloads took advantage of the 
et wa } well attended, - the program was 
Tour. Engineers Joint Council is a lectureship financed 4 
byt the National Science Foundation. Dr. Thijsse of Delft, Holland, has been 
eee to determine whether he would be available in view of the limited | 
funds available. Dr. Ippen will firm | up an itinerary if Dr. Thijsse is willing» 
: bi Note by Retiring Secretary. After 6 years’ service on the Executive Com- . 
va mittee, 2 years as Secretary, Harold M. Martin is to retire at the end of this hs E 


to the Society and the pleasant relationship with so many in the profession ae 


oenix Meeting ‘ 


hydraulie engineer are being planned for this ‘meeting. Better to go 


Claude R. Hunter (Soil Conservation 
P. Simmons, Jr. (Bureau of Reclammation) 
Kenneth R. Wright (Wheeler and Wri 


: — meeting of the officers is studded with interesting speakers and wal 


‘April 10-14, ASCE, Convention 
August 16-18, 1961 pare: Hydraulics Division © Conference, Urbana, Mlinois 


"ASCE, , Detroit Convention 


— 
= 
— 
have similar activities which would be | 
rado Section. No doubt other always welcome news | 
| j of interest to members of the Hydra ee — 
Hydraulic Division Activities 
— 
— 
— 
— 
= — tember through May at 7: p. m. lz: 
— 
— 
— 
— 
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2, 50 Meeting Geological Society of Amer- 


Denver, Colorado 
July 2 2, = Congress International Committee on 
Large Dams, Rome, Italy 
3- 1, 1961 Ninth IAHR Convention, Belgrade, Yugoslavia 
Seventh Midwestern Conference on Fluid 


a Mechanics and Soil Mechanics, Michigan a 
ot! bo University, East Lansing, Michigan 


18- 21,1962 Fourth National Congress of Applied Me- 
chanics, University of California, Berke- 
31 Much credit for the success of the Newsletter the past year is due to you 


ibuted timely news items and to the promptness of the Executive a 
- Committee and other committee members in furnishing pertinent information 4 
ic _ on Hydraulics Division participation at ASCE Conventions. . Your editor for a 
= 1960 expresses deep appreciation for the many timely items furnished, for 
¥ eae the Newsletter worthwhile and lighten the workload of the editor. =e 


‘Use of the aulics Division Newsletter 
are urged to use the ‘Division Newsletter for inquiries, 


2c ng news, committee reports, surveys and other items of interest to a =a 


Deadline dates for contributions: 
November 20; ‘March 1961 issue--January 20. 


Please direct all future contributions to your new Editor. 
bas 


Wheatridge, Colorado 
i 1960 Directory is now available to members on request. ‘The Direc- an 
_ tory lists the entire membership of the Society, giving the membership grade, _ 
Position, and mailing address of each. In addition, there isa 


geographical listing of the members is also included. 2S 


| 
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| 
| 
: 
| 
| 
a 
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meetings is particularly requ 
iF 
4 
| Bi 
a 
+ 
ay 
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It goes without saying that the in 


5 alue to every I member, and every member can obtain this valuable inform 


tion. ‘To receive your -free copy of the Directory siraply fill out the — 


ee Prompt delivery depends on prompt return of the coupon, 


: 


_ ASCE members are entitled to receive, free of charge, the 1960 ASCE a 
rectory. To obtain the directory simply clip this coupon and mail to: Ameri- 
can Society of Civil Engineers, 33 West 39th Street, New York 18, N. Y. a. 
2 cea make the mailing label legible—correct delivery d depends on n you —" 


“Address 


= 


ASCE- _WPCF JOINT SEWER M. MANUAL ISSUED 
A major addition ‘to the ASCE series of Manuals of Engineering 
is now available in a new volume entitled “Design and Construction of Sanitary _ 
and Storm Sewers.” Identified as No. 37, this publication is the result of serv-  __ 


- eral years of joint effort by the Sanitary ‘Engineering Division of ASCE and the : 


Industrial Wastes 

ied The twelve-chapter sewer manual contains 283 pages, over 100 illustra- a 
ae tions, 24 tables, and more than 100 references, As the first extended collec- 

E a tion of information on the subject, it will make a valuable reference in an im-_ 


portant phase of wastewater technology. Individual subjects covered include 
quantity of sanitary sewage and storm: water, hydraulics of sewers, , design of 

sewer systems, appurtenances and special structures, materials for sewer ‘a 

construction, structural requirements, | plans and 


8 may order 
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Hydrau cs | Division 


er copy. The price to members of the Water P 


he same ‘upon ee to their organization, 


4 . copy(s) of Manual is nig remittance of 


— ol Federation cc 
llution Control Feder 
— --------- 
New York 18. N. — 


